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PREFACE 


The topics with which the two companion volumes 

Patterns and Configurations in Finite Spaces 

and 

The Mathematics of Experimental Design : incomplete 
block designs and Latin squares 

are concerned have a long history, but the main stimulus to new devel¬ 
opments has come from the interest of statisticians in the efficient 
design of experiments. The origin of some of the patterns can be traced 
to recreational mathematics, or to problems in pure mathematics, par¬ 
ticularly in number theory. 

The statistical aspect of most of the subjects mentioned has been 
dealt with in many excellent textbooks, together with their analysis 
which leads to inferences about the effectiveness of treatments or 
other choices, the outcome of which is subject to stochastic variation. 
No such analysis is described in the present books. They contain 
rather the combinatorial aspects of the construction of designs, with¬ 
out regard to their practical application or indeed other usefulness. 

In Patterns and Configurations in Finite Spaces, the first chapter 
contains the algebra which is needed in the subsequent pages; it is, of 
course, not a course in algebra, but only a review of certain parts of it. 

Chapters II and III contain the fundamentals of finite geometry, 
two-dimensional and more-dimensional geometry being treated in sepa¬ 
rate chapters. In Chapter IV we reach configurations. 

In The Mathematics of Experimental Design we start with a short 
review of algebraic facts, and then deal fully with some special cases 
of configurations treated in the first-mentioned book. Thus, balanced 
incomplete block designs are described in Chapter II, and partially 
balanced incomplete block designs in Chapters IV and V, the last 
chapter being concerned with a special case of designs introduced in 
the former, namely those of two associate classes. 

Chapter III is devoted to orthogonal arrays and to a special case 

of these, namely Latin squares and their orthogonal sets. It contains 

an account of the most recent results in this field, including the proof 

ol the falsity of a hypothesis by Euler, which was long considered to 
be plausible. 


v 


The bibliographies mention only those papers which are referred 
to in the text of the respective book. A complete bibliography of the 
subject would be too extensive. 

The exercises are meant to test the reader’s understanding, and to 
serve the teacher of these topics. Both the reader and the teacher will 
welcome the inclusion of solutions to the exercises. Examples are 
also found throughout both books. 


Birmingham, 1967 


S. VAJDA 
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CHAPTER I 
INTRODUCTION 

In this book we shall be concerned with the mathematical aspect of de¬ 
signs which have been developed for the purpose of analysing first agri¬ 
cultural, and then also other statistical experiments. The analysis of 
the results of such experiments is outside the scope of the booK. We 
shall concentrate on the description and construction of the designs. 

Much of the mathematics which forms the basis of the construction 
of designs possessing those features which make them useful for experi¬ 
mental layouts has been dealt with in the companion volume [129], We 
start here by stating briefly, without detailed proofs, those facts which 
we need in the following chapters, and refer the reader for a rigorous 
treatment to the presentation in that volume. 

Finite groups 

A finite group is a set of a finite number of distinct elements, say 
a i> ...» a s > with a rule of operation for combining any two of them in a 
given order. The operation satisfies the following rules: 

(1) If both a i and a j are elements of the group (not necessarily distinct) 
then their combination, written a ; ., and also called their product, is 
also an element of the group. 

(2) Associative law: if a it a jt and a k are elements of the group (not 
necessarily distinct), then ( a.a.)a k = a^a.a^. 

(3) The group contains a unique “identity” element, say /, such that 
for every element a. of the group we have a. / = la. = a.. 

«• tit 

(4) For every element a i of the group there exists a unique element, a ; - 
say, such that a i a. = a j a < = /. We write a. = aT 1 . Of course, / _1 = /. 

We do not necessarily assume that a. = a. a i for any pair of ele¬ 
ments of the group. A group in which this equality holds is called com¬ 
mutative, or Abelian. 

The number of elements in a group is called its order, and the small¬ 
est integer m such that a? 1 = I is called the order of a i . Every element 
has an order, and it is always a factor of the order of the group. 
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A group which consists of all the powers of one element is called 
cyclic. It is Abelian. 

Two groups are called isomorphic if each element of one can be 
made to correspond to an element of the other in such a way that the 
product of two elements of one group corresponds to the product of the 
corresponding elements of the other group. The relationship between 
two isomorphic groups is called an isomorphism, and if the two groups 
are, in fact, the same, then the isomorphism is called an automorphism. 

If the product of two elements corresponding in an automorphism 
are multiplied, they produce again an element of the group. If the pro¬ 
ducts of two such elements are all different, then the isomorphism is 
a “complete mapping” (see [64]). A group of order At + 2 cannot have 
a complete mapping. 

The permutations of m symbols form a group of order m\ = 1.2 ...m, 
where m is called the degree of any permutation group of m symbols. 

A permutation group is conveniently described by a sucession of cycles 
of non-overlapping sets of elements. Each cycle, say (a be ... d ), indi¬ 
cates that the permutation changes a into b, b into c, ... and finally 
d back into a. 

Examples 

The permutation which changes the order 12345 into 21543 may be 
written as (12)(35)(4), where the last “cycle” (4) indicates that 4 re¬ 
mains in its original position. This cycle could remain unrecorded, 
thus: (12)(35). 

The permutations /, (12) (34), (13) (24), (14) (23) form an Abelian 
(but not cyclic) group of order 4 and degree 4. 

The permutations /, (1234), (13) (24), (1432) form a cyclic, and 
hence Abelian group of order 4 and degree 4. 

The permutations I, (12)(34), (13)(24), (14)(23) 

(234), (132), (143), (124) 

(243), (142), (123), (134) 

form a group of order 12 and degree 4, which is not Abelian. For in¬ 
stance, 

(12)(34).(234) = (132), but (234).(12)(34) = (124). 

Galois fields 

* s a set of s (> 1) distinct elements which form a 

commutate group with regard to a combination which we call, and write 

1 Wn ’ the ldenti ty elements being denoted by 0, and the inverse 
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of a i by — a i . The elements are also subject to another combination, 
called multiplication. With regard to the latter, a ,0 = Oa. = 0 for all 
a, .. Apart from 0, the elements form a commutative group under multi¬ 
plication, and the identity element of this group is denoted by 1, and 
the inverse of a i by aT 1 . 

The elements of a Galois field satisfy also the “distributive law” 

a Mj + a *) = (S' + a k) a i = a i a j + a k . 

A Galois field exists for every s which is a power of a prime, and 
for no other s. All Galois fields of the same order are isomorphic as 
regards addition as well as multiplication, and we may therefore say 
that, in an abstract sense* there exists only one Galois field for every 
prime power. We denote it by GF(s). 

Example 

GF( 4). We denote the elements by 0,l,y,z. They satisfy the 
following tables of 

Addition Multiplication 



0 

1 

y 

z 


0 

1 

y 

Z 

0 

0 

1 

y 

z 

0 

0 

0 

0 

0 

1 

1 

0 

z 

y 

l 

0 

1 

y 

z 

y 

y 

z 

0 

l 

y 

0 

y 

z 

1 

z 

z 

y 

1 

0 

z 

0 

z 

1 

y 


Finite projective geometries 

A finite projective geometry of m dimensions consists of the 
ordered sets (x 0 , x v ..., x m ), called points, where the x. are elements 
of the GF(s ) and not all x i = 0. The space is denoted by PG(jn,s ) 

(for “projective geometry”). A PG(2,s ) is also referred to as a finite 
projective plane. (Other types of such planes also exist (cf. [ 129].) 

The x i are homogeneous co-ordinates of the point, i.e. for any 

k ¥= 0, the point ( kx Q kx m ) is considered to be the same as the 
point 

Any of the x i can have s different values, but they must not all of 
them be simultaneously zero. Moreover, multiplying all of them by any 
of the s elements of the GF(s), except by 0, produces the same point. 
It follows that the number of points in a PG Qn,s ) is 

(s m+1 - 1) /(s-1) = V(m,0;s), say. 

T be PG(m,s) contains /c-flats, (k = 0,1,..., m - 1), defined as the 
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set of all points whose co-ordinates satisfy m-k independent linear 
homogeneous equations with coefficients from the GF (s), not all simul¬ 
taneously 0 within the same equation. 

Alternatively, a /c-flat in a PG(m,s ) consists of all points with co¬ 
ordinates ( a 0 x 00 + ... + a k ..., a Q x Qm + ... + a k x km ) with all ele¬ 
ments from the GF(s). The a. are not simultaneously all 0, and the 
matrix 



has rank k+ 1. Each /c-flat contains V(/c,0;s) points. In particular, 
there are s+1 points on a 1-flat, also called a line. 

A 0-flat is a point, and an(m - 1 )-flat in a PG(m,s ) is called a 
hyperplane. 

It can be shown (see [129]) that the number of /c-flats in a 
PG ( m,s ) equals 


(s m+1 -l)(s m -l)...(s m -* +1 -l) 
(s* + 1 — 1) — 1) ... (s — 1) 


V(m,k;s), say. 


The number of f-flats which contain a given /c-flat (k 1 ) is 


V(m~k-l,m — t-l;s) = V(m~k — l,t-k — l;s). 

For instance, there are (s m -l)/(s-l) hyperplanes through a 
point, and (s m-1 -1)/(s-1) through a line. 



Example 

The PG (2,4) is a finite plane with 21 points and an equal number 
of lines. On each line there are 5 points, and through every point there, 
are 5 lines. Through any two points there is precisely one line, and 
any two lines meet in precisely one point. 

We can also define a conic in a finite projective plane, as the set 
of ail those points whose co-ordinates satisfy a homogeneous equation 
of the second order. If the quadratic form is not the product of two 
linear forms, then the conic is called non-degenerate, and consists of 
s + 1 points. As regards their enumeration, see [l29]. 


Example 

A conic'in the PG(2,4) is, for instance, that consisting of the 
five points (1, 1, z) (0, l,y) (1, l,y) (l,z,y) (l,z,z). Its equation 
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is 

*o + x f + y(* 0 *i + X 0 x 2 + x \) = 0. 

In a PG(2,2 e ) the tangents of a non-degenerate conic meet in the 
same point, the “nucleus” of the conic. In the example above, it is 
the point ( 0,1,1). 

Finite Euclidean spaces 

A Euclidean space of m dimensions, denoted by EG(m,s), is the 
set of all points of a PG(jn,s ) whose first co-ordinate x 0 is not 0. 

Because the co-ordinates of the points in a PG(/n,s ) are homogeneous, 
we can take that first co-ordinate to be 1, so that the points of a EG(m,s) 
are given by the ordered sets (x,,... , x m ), where the x i are elements 
of the GF(s). 

There are s m such points, and the number of fc-flats can be shown 
(see [l29]) to be 

V(m,k;s ) - V(m-l,/c;s) = V(m- 1, k- l;s) s m ~ k . 

For instance, there are s(s m -l)/(s-l) hyperplanes. 

A/c-flat in a EG(m,s) contains s k points. In particular, there are 
now s points on a line (a 1-flat). 

The number of i-flats which contain a given /c-flat is the same as 
that in a PG(m,s), namely V (m-k- 1, t-k- l;s). 

'd&L. In a EG.(2.s) there are conics with s + 1, with s, and with s-1 
points which we call, respectively, ellipses, parabolas, and hyperbolas. 

They number, respectively, s 3 (s-l) 2 /2, s 2 (s 2 -l), and (s 5 -s 3 )/2. 

Difference sets, and systems of difference sets 

A set of integers d v ...,d k is said to be a “difference set mod v.” ' 

if the k(k-l) differences d i -d j equal any non-zero value mod v the 
same number (say A ) of times. It follows at once that k(k-l) = A (v-1). 

Example 

v = 7. 1,2,4 form a difference set, with A= 1. 

v = 11. 1, 3,4,5, 9 form a difference set, with A = 2. 

We speak of a system of difference sets 

d yv ..., d^ k 
d mV ..., d mk 
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if the d.. are elements of a group of order s, and such that if we com- 

t J 

pute all the mk ( k- 1) values d^ (i = lj ..., ^ / 2 ) of any two 

elements in the same row, tnen each of the s — 1 non-unit elements 
appears the same number of times. If this number is A, then we have 
mk(k-l) = A(s-l). 

Multiplication, and reciprocity, are to be interpreted in the sense 
of the rule of combination of the group. If the group is that of addition, 
and m= 1, then we are back at the definition of a difference set. 

Example 

k = 3, m = 2,s = 7. 124 

3 5 6. 
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CHAPTER II 

BALANCED INCOMPLETE BLOCK DESIGNS 


In this chapter we deal with a design consisting of two types of objects 
with an incidence relation between them. We shall call the two types 
of objects varieties and blocks ; the latter contain the varieties. These 
expressions are taken from agricultural field trials, where these designs 
were first applied to statistical experiments. The number of varieties 
will, as a rule, be denoted by v, and that of the blocks by b. 

No variety will appear more than once in a block, and no two blocks 
will be identical. Also, to exclude trivial cases, we assume v > 1 and 
b > 1- Every block will contain the same number (say k) of varieties, 
and all varieties will appear in the same number (say r) of blocks. It 
follows that bk = vr. 


We assume, moreover, that any pair of varieties appears together 
in the same number (say A) of blocks. 

Any variety, x say, appears r times, and therefore together with 
r (k- 1) varieties in the same block, though they are, of course, not 
all different. Since every one of the v-1 other varieties appears A 
times together with x, we have r(fc-l)=A(v-l). 

Such a design is called a balanced incomplete block design 

(b.i.b.d.). It is a special case of complete, or balanced, configurations 

(see [129]). In [l3l], Yates called b.i.b.d’s “incomplete randomized 
blocks”. 


If b - v, and hence r - k, then the b.i.b.d. is called symmetric. 

If any two of the blocks have the same number (say /i) of varieties 
in common, then the b.i.b.d. is called of linked block type. We shall 
see that symmetry implies linkage, and vice versa. 

A b.i.b.d. is called resolvable if its blocks can be grouped into 
sets of equal sizes, so that the blocks of each set contain between 
them all the varieties once; if, moreover, any two blocks from different 

c! “;: the Sam , e T bet ° f Varieties in commo ". ‘hen the b.i.b.d. is 
called affme resolvable. 

We give first a few examples. 


? 


( 
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Finite spaces 

The points of a PG ( m,s ), 
and its t-flats as blocks, form 
design are as follows: 


or of a EG ( m,s ), considered as varieties, 
b.i.b.d’s. The parameters of such a 


PG(m,s ) 


EG (m,s ) 

V (m,0;s) 

V 

s m 

V (m,t',s) 

b 

V(m,t;s)-V(m-l,tis) = 

V(m-l,t-l;s) 

r 

Vs) 

V(t, 0; s) 

k 

s c 

V (m-2,1-2; s) if f *1, 

X 

V(m-2,t-2;s) if t* 1 

or 1, if f = 1 


or 1, if t = 1. 


As special cases, we mention t = m — 1: 


PG (m,s) 


EG (m,s) 

(s ffl+1 -l)/(s-l) 

V 

s m 

(s m+1 -l)/(s-l) 

b 

s(s m -l )/ (s-1) 

(s™-l)/(s-l) 

r 

(s m -l)/(s-l) 

(s“-l)/(s-l) 

k 

s m — 1 

(s m_1 -l)/(s-l) 

X 

(s m ‘Ll)/(s-l) 


These designs, derived from PG (jn,s), are symmetric. 
For special values of s we have, when t = 1: 


s = 

: 2 


s 

= 3 


PG(m, 2) 


EGim, 2) 

PG(m,3) 


EG(m,3) 

2 m+1_i 

V 

2 m 

&(3 m+ 1 -l) 

V 

3 m 

: m+ 1 -l)(2 m -l)/3 

b 

2 m " 1 ( 2 m _l) 

( 3771 + 1_ 1 ) ( 3 m _ !) /16 

b 

3771(3771 _ i) /6 

2 m -l 

r 

2 m — 1 

Vi (3 m — 1) 

r 

Vi(3 m -l) 

3 

/c 

2 

4 

k 

3 

1 

A 

1 

1 

X 

1 


For general s, m = 2, t = 1 we have: 
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PG(2,s) 


EG(2,s) 

1 + s + s 2 

V 

s 2 

1 + 5 + S 2 

b 

s(s + 1) 

1 + s 

r 

1 + s 

1 + S 

k 

S' 

1 

A 

1 


Conics 

From the enumeration of non-degenerate conics in a PG(2,s), 
given in [l29], it follows that if we take the points as varieties, and 
the non-degenerate conics as blocks, then we obtain a b.i.b.d. with 
the following parameters: 

v = s 2 + s + 1, b = s s - s 2 , r = s 2 (s 2 - 1), 

* = s + 1, A = s (5-1). 

If we take parabolas as blocks, and the points of a EG ( 2 ,s) as 
varieties, then we have a b.i.b.d. with 

v = s 2 , b = s z (s 2 -1), r = s (s 2 - 1), k = s, A = s (s -1), 
and if we use hyperbolas, we have 

v = s 2 , b = l / 2 s 3 (s 2 - 1), r = l / 2 (s + 1) s(s - l) 2 , 

k = s-1, A = H s(s- l)(s-2). 

More generally, using the results for non-degenerate quadrics 
given in [l29], we obtain b.i.b.d’s with 

v - (s 2t+ '- l)/(s- 1), b = s tl!+ " n (S 2i+1 - 1), 

i=1 

* = (^-l)/(s-l), r = 5* (fc+1) (s 2e - 1) iT (s 2<+1 -l) 

i= 1 

A = s tlt+ " (s 2 *-' _ 1) ff( s «+'_l) 

i =1 

from points as varieties and quadrics as blocks in a PG (2 t,s) where 
t ^ 2, or 

v = (s 2i -l)/(s-l), b =Vis t2 (s t + 1) n 1 (s 2 *+i -1), 

i=1 
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/ t = (s t -l)(s < - , + l)/(s-l), r = 54s*V"'+ 1) n (S 24t, -1), 

i= 1 

K = y 2 s‘ 2 (s , -'+ l)(s 2t_ ^+s t- ' - s‘- 2 - 1) ff(s*‘ +l - 1) 

i =1 

from ruled quadrics as blocks in a PG (2t-l,s), t > 2, and 

V = (s“- l)/(s - 1), b = fcs'V - 1) n' (s 2<+1 - 1), 

k = (s 4 + l)(s* _1 - 1)/(s - 1), r = tfs* 2 (s t-t - l)TlV i+1 - 1). 

fc = 1 

e-2 

If, in the last two cases, t = 2, then the term II s( 2i+t - 1) is re¬ 
placed by 1. * =1 

If we use hyperbolas, we have v = s 2 , b = s 3 (s - l)/2, r = 

(s + l)s(s - lf/2 and k = s - 1. 

A = 1 / 2 s e2 (s t ' 1 - lKs 24-2 - s 4_, + s t_2 - l)Ti 2 (s 2i+1 - 1) 

i— 1 

from unruled quadrics as blocks in a PG(2t — l,s), t > 2. 

Consider now a non-degenerate conic in a PG (2, 2*) = PG (2, s). Of 
the s 2 + s + l lines of the plane, ^ ^ ^ meet the conic in two points and 

do not pass through the nucleus, s + 1 are tangents, and pass through 
the nucleus, do not have any point in common with the conic. We 

call these lines “v-lines”. 

Of the s 2 + s + 1 points of the plane, one is the nucleus, s + 1 are 
points of the conic, s 2 - 1 points are neither on the conic, nor is one 
of them the nucleus. We call these points “b-points.” 

We take the b-points as blocks (not as varieties), and the v-lines 
as varieties. We have then a b.i.b.d. with 

v = CD- i ” s2 - 1 - *-*». 

Because any two v-lines meet in a single point, A = 1. 

This b.i.b.d. is resolvable (though, of course, not affine). This 
can be seen as follows. 

Take a line which is not a v-line. There are on it s - 1 different 
b-points. Each of these corresponds to a block with V 3S varieties, and 
two varieties in different blocks cannot be the same, because the line 

connecting two of the b-points is not a v-line. Hence all these blocks 
contain, between them, all varieties. 

“ we select any point on the conic, or the nucleus, and consider 
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the s + 1 lines through it, then there pass through such a point s + 1 
lines which are not v-lines, and thus we obtain s + 1 different sets of 

’ thltThe 63 ’ ° h ° f theSe S6tS C ° ntalns 811 valieties - We see 

also that the grouping into sets can be done ins + 2 different ways 

because there are s + 1 points of the conic from which we can start’ or 
we can start with the nucleus. ' 

In [9] R.C.Bose has shown that resolvable b.i.b.d’s exist for all 
parameters of the form v = s * + 1, b = s* (s* - s + 1), , = s 2, 

k = s + 1, X = 1, when s is a prime power. He takes as varieties the 
points of the curve + *f+< + . 0 in a PG( 2 >s «) and as blocks 

he lines which are not tangents of this curve. For s = 2, we obtain 
trio D.i.b.d. given on page 27. 

We shall now describe two methods of obtaining b.i.b.d’s from a 
given symmetric one, with parameters v, r, and A. 


Block section 

., r: deSign ’ CaIled “ residual ’- « obtained by omitting one 
block and all its varieties from all blocks. The new parameters are: 

v'=v-k, b'=v- 1, r'= r, k'= r - A, V= A. 


Block intersection 

The other new design, called derived, is obtained by omitting one 
block and retaining in the other blocks only those varieties which were 
contained in the omitted block. The new parameters are: 


= k, 


= v-1, 


" = k- 1, k" = x, A"= A - 1 


A method of construction, introduced in [ 4 ] (see also [ 5 ] and 
will be described by an example, taken from [ 22 ]. 

whichT eUeS T t den ° ted by numbers mod 7. n is an element 

blocks 4 Change When any number is added to it. Consider the 


(1,1 2,1 4,1), (3,1 1,2 5,2), (6,1 2,2 5,2) 
(5,1 4,2 6,2), (0,1 0,2 n, 0). 


d g ' res P ectl vely, 1, or 2, ..., or 6 to the first element of a pair and 
e ucing mod 7, we obtain a resolvable b.i.b.d. with v = 15 b = 35 
* = 3, r = 7, A = 1. 
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Incidence matrix 

We now introduce a way of representing a b.i.b.d., its so-called 
incidence matrix. This is a matrix of v rows and b columns. Its rows 
correspond to the varieties, and its columns to the blocks. We number 
the varieties and the blocks, and if the ith variety appears in the j th 
block, then we write into the ith row and the /th column a 1, otherwise 
we write there 0. 

It is easily seen that the sum of entries in any row is r , and that 
in any column is k. The inner product of two different rows is A, and 
if the b.i.b.d. is of linked block type, then the inner product of any two 
different columns is fi. 

We denote the incidence matrix by A. Its transpose, A', is the 
incidence matrix of the so-called “dual” design, and if we write 1 for 
0, and 0 for 1, then we obtain from A the incidence matrix of the “com¬ 
plementary” design. 

We quote from [129] that all rows and columns of the matrix product 
AA ' have the same total, namely rk, and that the determinant \AA ; | 
equals rk (r - A) v_1 . 

Because r > A, this is positive, i.e. AA 1 is of rank v. The v by b 
matrix A cannot have smaller rank, hence 

b > v. 

This proof is taken from [6]. R.A. Fisher has proved this inequality in 
[43] by a different method. 

As an example of the use of the concept of the incidence matrix we 
now prove that a b.i.b.d. is of linked block type if and only if it is 
symmetric. 

Denote the square matrix of order v whose elements are all 1 by J v 
and the identity matrix of order v by I y . If v = b, then we have 

AA'= (/c-A)/ v + A J y , and AJ V = J v A = kJ v = rJ v . 

Hence A'A = A~\AA')A = A~\k- A)/ v A + A AT% A 
= ( k-X)I v +\J y = AA '. 

This proves that a symmetric b.i.b.d. is of linked block type, and 

A = (jl = k(k - 1 )/(v - 1). 

Conversely, let a b.i.b.d. be of linked block type. Then 



AA' = 


k\i .. (i 

[xH •’ k 
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r\..\ 
A r.. A 

AA.. r. 


and A'A = 


I AA) - rk(r ~ A) v_1 and \A'A\ = [k + n{ b - 1)] (k - fz) 6 * 1 . 

Both these determinants are positive, and hence the rank of A A' is 
v, and that of A A is b. The rank of a product of two determinants 

cannot exceed that of either factor, and therefore we have b > v, and 
v > b, i.e. 

v = b, r = k, A = /x. 

In particular, let A = fi = 1. Then we have v=b = r 2 — r + 1, and 
denoting r - 1 by s, we have v=b = s z + s + 1, r = k = s + 1. 

Thus the blocks and varieties of a symmetric b.i.b.d. with A = 1 
satisfy the conditions for lines and points of a finite projective plane. 
We have already seen that the converse is also true: the two concepts, 
namely that of a symmetric b.i.b.d. with A = 1, and of a finite projec¬ 
tive plane, are equivalent. 

If the b.i.b.d. is resolvable, then we can derive b ^ v + r — 1, 
which is a stronger inequality than b > v. Taxe those columns of 
A which correspond to the blocks of one complete replication. The 
entries of each row add up to 1, and therefore not more than b-r+1 
columns can be independent. Hence b-r+1 ^ v, and the inequality 
follows. It was proved by a different method in [l29]. 

If the b.i.b.d. is affine resolvable, then we have r sets, oin = b/r 
blocks each. Call the number of common varieties in any two blocks 
from different sets /z. Take one block, say B 0k and consider all its 
varieties, and all their replications in the other sets. In the sets which 
do not contain B 0 , these varieties will appear k (r - 1) times, because 
every one is repeated r - 1 times in them. Also, since the n(r - 1) 
blocks in these sets each contain /x varieties in common with those of 
B 0 , we have k(r-l) = n(r-l)^z, i .e. p = k/n = kr/b. 

Then, consider the i^\ pairs of varieties in B Q . They appear 

A - 1 times in other blocks. In every one of these, n(r-l ) in number, 
there will be such pairs. Hence Q (A -1) = n(r-l). 

Replace n by b/r and /z by kr/b. Then we obtain, after some 
straightforward arithmetic, b = v + r - 1. 

From an orthogonal matrix of order At whose elements are 1 or -1 





14 


we can construct a b.i.b.d. with parameters v = b = 4t- 1, r = k = 2*~1, 
A = t _ 1 by a method indicated in [128]. 

We can change all signs in a row or in a column of an orthogonal 
matrix in such a way that all elements in the first row, and also in the 
first column, are 1. Orthogonality is not destroyed by this change. We 
then suppress the first row and the first column. As a consequence of 
the orthogonality of the original matrix each row will then contain 2r —1 
elements 1, and any two rows will have t — 1 pairs of 1 s in the same 
column. If we change — 1 into 0, then we have obtained an incidence 
matrix of a b.i.b.d. with the parameters as stated. (The procedure can 
be reversed, and we can construct a Hadamard matrix from a b.i.b.d. 
^ose parameters.) 

A design with such an incidence matrix is called a Hadamard con¬ 
figuration. An example is the PG(jn, 2), with hyperplanes as blocks, and 
points as varieties. To see this, write t for 2 m_1 . The complementary 
design, i.e. that whose incidence matrix is obtained by exchanging 0 
and 1, is also called a Hadamard configuration. Its parameters are 

v = b = 4t - 1, k - r - 2t, A - t . 


1 


Existence theorems 

The parameters of a b.i.b.d. satisfy the conditions bk = vr, r(/c — 1) 
= A(v- 1), and hence = ^ therefore necessary that 


r = A (v - 1)/(/c - 1) and b = A v(v - 1 )/k (k - 1) be integers. There 
are cases where it can be shown that these two conditions are also 
sufficient for a b.i.b.d. with such parameters to exist. We shall discuss 
a few such cases. 

When A = 1 and k = 2, then we obtain r = v - 1 and b = ) • In 

this case the blocks are all possible pairs of varieties. 

When A = 1 and k = 3, then 3 b = vr and 2r = v-1. These are the 
conditions for a Steiner triple system, and it is shown in [129] that 
they exist if and only if v = 6i + l or 6t + 3, where t is a non-negative 
integer (cf. [96] and [68]). 

When A = 2 and k = 3 (a case studied by R.C. Bose in [4]), then 
3b = vr and 2r = 2(v-1); hence b = r(r + 1 )/3 and v = r+1. There¬ 
fore either r = 3t, v= 3f+1, b = t(3t + l), or r = 3* + 2, v = 3(*+l), 
b = (f + l)(3f + 2). If t = 2u, then r=6wor 6u + 2. These designs can 
be constructed by combining two systems of Steiner triples with the 
same number of elements, either permuting their notation, or taking non- 


isomorphic systems. If t = 2u+l, then r = 6u + 3 or 6u + 5. These de¬ 
signs can also be constructed. 

In these cases we can express our results by saying that it is suf¬ 
ficient for the existence of a b.i.b.d. that the expressions for r and b 
at the beginning of this section should be integers. These same con¬ 
ditions are also sufficient for k = 3 or 4, with any value of A, and 
also for k = 5 when A = 1, 4, or 20, with the possible exception of 
v = 141 in the first of these cases. This was proved in [49]. That for 
k = 6 and A = 1 the integrality of the above expressions is not sufficient 
was shown in [l26]. 

A symmetric b.i.b.d. with v = 43, r = 7, A = 1 is impossible, because 
the equivalent PG (2,6) does not exist. It follows that a b.i.b.d. with 
v = 36, b = 42, r = 7, k = 6 cannot exist either, because it would be 
equivalent to a EG(2,6), which could be extended into a PG(2,6). 

For k = 5 and A = 2 the integrality of the expressions for r and 
b is not sufficient, as shown in [74]. On the other hand, considering 
special values of v, it is shown in [56] that one b.i.b.d. exists (unique 
apart from isomorphisms) for v = b = 11, r = k = 5, A = 2 and that three 
non-isomorphic b.i.b.d’s exist for v = b = 16, r = k = 6, A = 2. 

Designs with v = b = 37, r = k = 9, A = 2 are mentioned in [ 2 9] and 
they are studied in [59]. 

Nandi has shown in [75] that there exists only one design with 
v = 6, b = 10, r = 5, k = 3, A = 2. It appeared first in [ 131] and it is 
probably the very first b.i.b.d. to have been published. We quote it here: 

012 013 024 035 045 ,j^<^ 

125 134 145 234 235. ^ ' 


Nandi has also shown that there are three designs for v = 10, b = 15 

r = 6, k = 4, A = 2 and in [76] that for any v smaller than 15 there exist' 
only one symmetric b.i.b.d., if at all, while there are five for v = b = 15 

\ = 3 ' Fr ° m the latter ' four ^n-isomorphic b.i.b. d’s can be 
formed by block section, and four by block intersection. 

ceDt!,f ,t a11 "T derfVe existence theorems b V “«ing use of the con- 
cept of the incidence matrix. 

If the b.i.b.d. is symmetric, then, as we have seen, IAA'1 - Ml* 

= rHr - Therefore we have 1 1 ~ W 


(TV°T " J Iiv=bis even < the ” r- A = A - A must be 
is has also been pointed out in [112], [33], and [105].) 


a square. 
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Thus, for instance, b.i.b.d’s witn v - b - 22, k - r 7, A - 2 or 
with v = b = 46, /c = r=10, A = 2 are not possible. The former fact is 

also proved in [57]. 

For odd v, a necessary condition for a symmetric b.i.b.d. to exist 
is given by Chowla and Ryser in [33]. To prepare us for its proof we 
show first 

Theorem II.2 If v, k, and A are parameters of a symmetric b.i.b.d., 
and v is odd, then the equation x 2 = (k - A)y 2 + (-l) 2(v ’'Ax 2 has 
an integer solution where x, y, and z are not all 0. 


Proof: Let the incidence matrix be A = (a !y ) and consider the expres- 

V 

sions Lj = Z^a^Xi. By the definition of k and A we have 


2 l3j = (A-A) 2 xf + a(2 X.) 
1 v ' i=\ % \i=i V 



We refer now to Lagrange’s theorem which states that every positive 
integer is the sum of (not more than) four squares. 

Let k -A = af + a\ + a\ + a 2 . Then 

(k- A) = (x2 + x|+x|+x2) = (y 2 +yf+yf+y 2 ) 

where 

y, = ^1^ Xj H" ^3 ^3+ ^2 = ^2 ^2 ^4 ^3 ®3^4 

y 3 = - a 3 x, + a 4 x 2 + a, x 3 -a 2 x 4 y 4 = a 4 x, + a 3 x 2 - a 2 x 3 - a, x 4 


The matrix which transforms the x - s into the y - s is not 
singular. Its determinant is, in fact, (k - A) 2 . Thus the L. are 

If 

V 

linear rational functions of the y. , say L i = 2 h i . y.. 

Conversely, the. x i are linear rational functions of the y. . The 
equation (L) is therefore an identity in the y- and it holds whichever 
relations we establish between them. 

At this stage we have to distinguish between the two cases (i) 
-v = 1 (mod 4) and (ii) v s 3 (mod 4). 


* Simple arithmetic shows that if k — A is a square, say m^, then 

Lj = mxj +2(fc - m)x i /v 

satisfies this equation. The coefficients of the x i are rational, but to obtain 
an incidence matrix they would have to be 0 or 1. 
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Case (i) 

v - 1 is divisible by 4, so that we can write 


1 L 2 j 

j =1 3 


= (y , 2 + y 2 2 + y 3 2 + y 4 2 ) + ... + (>^_ 


4 +yl-3 + y 2 v -2+ yl-\) 


+ 


(*-A.)y v 2 + x(.i 


2 


where the y^. are linear rational functions of the x i ( and x y = y v ). 

If b n * 1, then let y, be such a rational function of y 2 ,y y that 
L, = y,; if = 1, such that L, = -y^ In either case L 2 = yf. 

We continue treating L 2 ,..., in the same way, thereby making y 2 
a rational function of y 3 ,... ,y v and so on, until L, = y,,... ,L , = 
y v _ 1 have cancelled out and we reach 

4 = (* - + X ( £ . 

Now both sides are rational functions of y, and we choose it to be 
different from 0. It follows that x 2 = (A - A)y 2 + Az 2 can be solved in 
rational, and hence in integral values, with y # 0. 

Case (ii) 

v + 1 is divisible by 4, and introducing y v+1 = x y+1 we can write 





* + yl+y* + y?>+-+(y *. 2 


+ y 2 1 

•V-1 


+ y v + y 


v-M 



Proceeding as we did in case (i), we obtain eventually 

<*-*>*?+, = y. 2 t i + • 

Hence x 2 = (A; - A)y 2 - Az 2 can be solved in integers, not all of 
them zero. 

This completes the proof of Theorem II. 2. 

Examples 

(1) Consider the case of a finite projective plane. Then v = 

ca + + (/C “ 1)2 iS ° dd f ° r any k ' The Di °P ha ntine equation is in 

(A n a’®* V S 1(m ° d 4) ' i>e ’ k ~ 1 B ° or 3 < mod 4 ). * 2 = 

~ L )y + z and can be solved when k - 1 = 4f by x = 2t + 2, y = 2, 
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z = 2t-2 and when k -1 = At + 3 by x = 2t + 2, y = 1, z = 2r + 1. 

In case (ii), i.e. when v = 3 (mod 4), i.e. k-1 = 1 or 2(mod 4) it 
is x 2 + z 2 = (k-l)y 2 . 

This can only be solved in integers if k-1 is the sum of two 
squares, which excludes the values k = 7,15, 22, 23, etc. It is not 
known whether any other values, not belonging to this series, are also 
excluded from leading to a finite projective plane. On the other hand, 
it is easily seen that the rule does not exclude any prime powers: 

2 n = 2 (mod 4) only if n = 1; any prime of the form At + 1 is the sum 
of two squares; ( a z + b z ) zt is a square, and 

(a 2 + h 2 ) 2t+1 =[a^ + b z ) t ] 2 + [b(a z + b 2 ) t f. 

(2) For the two Hadamard designs the resulting Diophantine equa¬ 
tions are x 2 = ty 2 -(f-l)z 2 and x 2 = ty 2 -tz 2 . The first is solved by 
x = y = z = 1, and the second by x = 0, y = z = 1. 

We prove now the theorem of Chowla and Ryser, namely 

Theorem 11,3 Let a symmetric b.i.b.d. with parameters v, k, A with odd 
v be given. If a prime p divides the square free part of k -A, then 

(-1)* <V_1> A is a quadratic residue of p. 

Proof: If the assumption of the theorem holds, then the Diophantine 
equation of Theorem II.2 can be written 

j X 2 = pty* 

where t is square free and relatively prime to p, and where y is the 
previous y multiplied by the square factor (if any) of k — A. We can 
choose x, y, and z so that they do not have a common factor. Then 
p does not divide z, because if it did, then it would also divide x, 
p 2 would divide x 2 , and hence pty 2 . But then p would have to divide 
y as well. 

Because p does not divide z, (-l)^ (v_1) A is a quadratic residue 

of p. 

Example 

v=b = 29, r=k = 8, A = 2. Here 3 divides 6, and (-1 ) i(v_1) .2=2 
is not a quadratic residue of 3. Therefore a b.i.b.d. with the given 
parameters does not exist (cf. also [55]). 
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Note that the theorem is trivial for p = 2, since then both 0 and 1 
are quadratic residues. 

By methods outside the scope of this book it is proved in [26] that 
if s = 1 or 2 (mod 4), and if the square-free part of s contains at least 
one prime factor of the form 4*+ 3, then no b.i.b.d. with v = b = l + s + s 2 , 
r = k = 1 + s, A = l(i.e. a finite projective plane P(2,s)), can exist. 

The following similar theorem refers to resolvable configurations. 


Theorem II.4 If in an affine resolvable b.i.b.d. with parameters v, b, r, 
k, A the number v of varieties is odd, then if r is odd, k must be a 
square, and if r is even, v must be. (In the latter case, since k z /v is 
an integer, this is also a square.) 


Proof (cf. [63]): Adjoin to the incidence matrix r - 1 rows to obtain a 
square matrix, A say. Let the ith new row have 1 in positions (i-l)n+ 1, 
(i-l)n+2, (i — l)n + n, and q i n a u other positions. Then the de¬ 
terminant \AA'\ equals 


r A ... All ... 1 
A r ... A 1 1 ... 1 

A A ... r 1 1 ... 1 

11 ... 1 n 0 ... 0 

... 1 0 n ... 0 


and can be computed by applying the following steps: 

« add the sum of the 2nd to vth row to the first (and remember that 
' + A(v-1) = rk)\ 

lhat?“ fir3t COh,mn £r ° m the 2nd> third .^ < a " d member 

'-1 columns; 1A UmeS 11,6 SUm ° f ‘ he 2 " d t0 V ‘ h C ° lumn £rom the las ‘ 
row SUbtIaC ‘ V/n = k ‘ imeS the SUm 0f £he last '- 1 from the first 


We obtain a determinant with zeros 


in the triangle above the diagonal, 



20 


and its value is 

kr'n'-Xrk-k(r-D) = KV' 1 - * r ' f+V '' 

This determinant must be a perfect square, and the theorem follows. 
An equivalent result is proved in [114], using results from [36], 


Structure 

We turn to the investigation of the structure of b.i.b.d s, in par¬ 
ticular as regards the dependence of sets of blocks on others. 

Order the blocKs of a b.i.b.d. with incidence matrix 4 in some way 
and consider only the first t of them. In the incidence matrix, the first 
t columns refer to these blocxs, and we call the incidence matrix con¬ 
sisting of these t blocks 4 0 . 


Consider the matrix 




say, where 0 is the t by b-t 


zero matrix. We have 



In [36] W.S. Connor has shown by simple arithmetic that the deter¬ 
minant equals /cr 1_fc (r-A) v_t_1 \C t | . Here C t is called the char¬ 

acteristic matrix of the chosen t blocks, and can be written 


XkJ t + r(r-\)I t - tAqAq. 


Because |4^^ | ^ 0 if the number of columns is not smaller than 
that of the rows of 4,, and because it is =0 otherwise, we see that 
\C t \ ^ 0 if b > t + v, \C t \ = 0 if b < t + v while, if b = t + v, |4,4j| 
is a perfect square, since all elements of 4^ are integers. 

When t = 1, then |C,| = A/c + r(r - A) - rA = (r-A)(r-fc) > 0, 
from which we obtain again the inequality r > k or b > v. 

When t = 2, then \C 2 \ = (r-A) 2 {r-k ) 2 - (U-n,/) 2 ’ >0. 

For a symmetric b.i.b.d. this results in n 12 = A, also a Known result. 
Another example is given by v = 16, b = 24, k = 6, r = 9, A = 3 which 
leads to 0 ^ ju < 4. 


Connor has studied in particular the b.i.b.d’s with v = 


( v —• luc u.i.D.a s witn v= i 

.2 J, r.k + 2, A = 2 (the case foi k = 3 is quoted on page 15), 

whichT 3 bl ° Ck - B " and Ca “ "* ‘ he " Umber of those othei blocks 
hich have i varieties in common with B,. Then 
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b -1 = V2k(k + 3), 


2 in. 

i=o * 


k(r-l) = k(k+l) 


It follows that 



J o («-l)(i-2) n< 


= 2 * (i-1) n i - 2 2 w. + 2 2 n . 
t=o 1 <=o z i=o 1 


= (& -1) - 2/c (k +1) + & (/c + 3) = 0. 

Because i as well as n i is non-negative, all terms (*-l)(i-2)/i i 
= 0, since otherwise their sum could not be zero. Hence n i = 0 for all 
i except for i = 1, or i = 2. For these two values we have, from n y + n z 
= Ytkik + 3) and n, +2« 2 = /c(* + l) that - 2fc, and « 2 = y 2 ^-l). Thus 

the matrix A A has k in all diagonal positions, and either 1 or 2 in all 
other positions. 

By a similar arithmetic Connor establishes, starting from two ini¬ 
tial blocks, that if these have t(= 1 or 2) varieties in common, then 
there exist 

(k-t)(k-t + l) + kt - y 2 k( k+3 ) + 1 blocks with two varieties from 
each of the initial blocks, 

2t(k-t) blocks with two varieties from one of them and one variety 
from the other, and 

k(2-t) - 2t(l-t) blocks with one variety from each of the two 
initial blocks. 


EXERCISES 


1) Construct a resolvable b.i.b.d. from the conic in the PG(2,4) given 
in the example on page 4 by the method explained on page 10. 

(2) Construct a design which is complementary to the PG (2,2). 

(3) Construct a Hadamard matrix from the PG(2,2 ). 

(4) Illustrate Theorem II.2 in the case of a PG(2, 2). 


p.w- 
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CHAPTER III 

LATIN SQUARES AND ORTHOGONAL ARRAYS 
Latin squares 

Consider n 2 cells, laid out in a square pattern, and containing n 
different symbols, each repeated n times, in such a way that the cells 
in each row, and also in each column, contain each symbol precisely 
once. Such an arrangement is called a Latin square (L.Sq.), because 
L. Euler, who studied such patterns in [ 40 ], used Latin letters for sym¬ 
bols. We call n the side of the L.Sq. 

An obvious way of constructing a L.Sq. is to write down the n 
symbols in some order and then to permute them cyclically to obtain 
other rows. A L.Sq. with side 5, constructed by moving each line two 
positions on to obtain the next, is attributed to the Norwegian Knut Vik 
in [ 44 ] • It is also easily seen that the elements of a group, written in 
the form of their multiplication table, form a L.Sq. 

Starting from a L.Sq. and permuting symbols, rows, or columns, or 
by reflections, we obtain in general new L.Sq’s. However, we do not 
always obtain, by these devices, all L.Sq’s of a given side. We shall 
not deal here with questions of enumeration, but mention a few perti¬ 
nent points. The symbols can be chosen in such a way that they appear 
in a prescribed order in the first row, and so that they appear in the 
first column in the same order. We say, then, that the L.Sq. is in Stand¬ 
ard Form, and n!. (n - 1)! L.Sq’s are represented by the same Standard 
Form. 

There exists just one L.Sq. in Standard Form for n = 2, one for 
n = 3, and four for n = 4, namely: 


0123 

0123 

0123 

0123 

1230 

1032 

1302 

1032 

2301 

2310 

2031 

2301 

3012 

3201 

3210 

3210 


For higher n the numbers increase rapidly. They are: 56 for n = 5 
(all found by Euler), 9408 for n = 6 (see [ 45 ]) and 16,942,080 for 
” = 7 (see [77] and [ 104 ]). 
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Further information can be found in [ 77 ] and in Chapter 4 of [ 48 ]. 
Interesting facts are also given in [ 9 s]. 


Latin rectangles 


A Latin rectangle of c columns and n>c rows is an array of n 
symbols such that each column contains all the symbols, and no row 
contains any symbol more than once. (This pattern is also called a 
Youden square.) If we take the rows as blocks, then we obtain a sym- ——~ 
metric configuration. The complementary design can also be written as 
a Youden square, as was proved in [129], and the two designs can be 
combined into a L.Sq. In other words, any Latin rectangle can be ex¬ 
tended into a L.Sq. (cf. [ 47 ]). For the number of ways in which this can 
be done, see [ 39 ]. 

We shall now show that we can also extend an r by c rectangle, 
containing n symbols in such a way that no symbol appears twice in 
any row or in any column, into a L.Sq. by adjoining n-c columns and 
n-r rows, provided the number N{i) of times that the symbol i appears 
in the rectangle is at least r+c-n, for all symbols. 


The necessity of the condition is obvious. If we want to extend the 
given rectangle into a Latin rectangle of n rows and c columns (as we 
must, at least implicitly) then we shall have, in the latter, c repetitions 
of every symbol. But in the additional n - r rows we cannot include any 
symbol more than n-r times, hence it must already appear c+r-n 
times in the original pattern. 

To prove the sufficiency, form the sets Sj , consisting of those 
elements which are not included in the yth column (y = l, ..., c) of the 
original pattern. Each of these sets consists oi n-r elements, and 
each symbol will, altogether, appear M(i) = c-N(i ) times. To write 
down these S,(y = 1 , .... c) as complementary to the c columns in the 
original pattern, we must be able to choose a system of c distinct ele¬ 
ments, one from each set, which we then use as a further row. Such a 
system exists by P. Hall’s theorem (see [129]), because any k sets of 
contain altogether k(n-r ) elements, none more often than n-r 
times, in view of M(i) = c-N(i)4c-(c + r-n) for all i. Continuing, 
we can thus construct n — r additional rows to complete an n by c 

Latin rectangle, and from this we can complete a Latin square, as we 
nave seen. 

This proof is contained, in a slightly different form, in [ 103 ]. 
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Example 

n = 6, c 

123 

= 3, r 

i 

= 4 

W) 

/v(0 

i 

246 

1 

2 

1 

4 

352 

2 

0 

3 

5 

465 

3 

1 

2 

6 

missing 511 

which 

can be written 

531 


634 614 


M(i) 

1 

1 

1 


/v(0 

2 

2 

2 


The Latin rectangle 






N 7 


* frt 


4^ 


v? v. 


123 

can be completed by adding 

456 

246 

(for instance) 

315 

352 


641 

465 


132 

531 


264 

614 


523 


As a consequence of the results just derived it is shown in [4l] 
that an incomplete L.Sq. of side rc(>4) with t different elements 
(i.e. a square of side n such that a subset of the n z places are occu¬ 
pied by numbers from among 1, 2, n, and no element occurs twice 
in the same row or column) can be embedded in a L.Sq. of side t, if 
t > 2n. 


Orthogonal sets of Latin squares 

In what follows, we shall be interested in pairs of so-called ortho¬ 
gonal Latin squares, i.e. L.Sq’s of equal sides and such that their 
superposition produces ordered pairs of elements, each of which occurs 
X onlyjjnce-in the same cell. For example, the two L.Sq’s of side 3, 




012 

012 



120 

201 



201 

120 


are orthogonal. Such a pair is also called a Graeco-Latin square, fol¬ 
lowing Euler who wrote the symbols of the first in Latin and those of 
the second in Greek letters, thus: 


aa 

b/3 

Cy 

by 

c a 

a/3 

c/3 

Uy 

ba 
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A set of L.Sq’s such that any two of them are orthogonal is called 
an orthogonal set. Some writers use such expressions as “mutually” 
or -pairwise orthogonal” L.Sq’s. Others, imagining all of them super 
imposed, speak of “hyper-Graeco-Latin” squares. 

There cannot be more than n -1 orthogonal L.Sq’s in any set. This 
is seen as follows. 

The symbols in each square can be permuted, independently, with¬ 
out destroying orthogonality. The first row in all the squares can there¬ 
fore be written as 0 1 . . . n — 1. 

The symbol 1 will in any pair correspond to itself in the top 
row. In the first column, this symbol cannot occupy the top position 
in any of the squares, because this position is already occupied by 0. 
Furthermore, it cannot occupy any other position in that column mote 


than once taking all the squares together, because it already corres¬ 
ponds to itself in the first row. Hence, not more than n - 1 squares 
can be constructed to form an orthogonal set. If a set contains precisely 
n-1 orthogonal L.Sq’s then it is called complete. (We shall see in 
Chapter V that the existence of n -3 pairwise orthogonal L.Sq’s of 
order n > 4 implies the existence of a complete set.) 

When n equals a prime power, s say, then a complete set can be 
constructed from a finite projective plane. 

Choose a particular line L. Through each of its s+1 points there 
pass s lines, apart from L. We call these a pencil of lines, and in 
each of the s+1 pencils we number the lines 0, 1, ..., s-1. Through 
each of the s 2 points outside L there passes one line from each of the 
s+1 pencils. 


Denote the points on L by P 0 , P,,..., P s . Select two of these, say 
P ° and p i > an d let the numbers of the lines through them refer, respec¬ 
tively, to rows and columns of a square arrangement. Enter in the cell 
of row i and column j the number of that line from P 2 which passes 
through the point defined by the line i from P 0 and line / from P,. 

Because any two lines meet in one single point, we obtain in this way 
a L.Sq. 


Construct other L.Sq’s by using P 3 instead of P 2 , and so on. It 
1S easi ly seen that-again because two lines meet in a single point- 
any two of the s-1 L.Sq’s are orthogonal. 

Conversely, if n — 1 pairwise orthogonal L.Sq’s exist, then a finite 
Projective plane with n+ 1 points on each line also exists. To prove 
ls » It is sufficient to show that a b.i.b.d. with v= n 2 , b = n(n+ 1), 

~ n +1,/c = n, X = 1 exists, i.e. a EG(2,n) which can be expanded 
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into a projective plane. 

For the purpose of this proof, number the cells of the L.Sq’s from 
1 to n in the first row, n+1 to 2n in the second, etc., until n 2 —n+ 1 
to n 2 in the last row. Then form blocks, by putting into the same block 
the numbers of those cells which have the same element in the 1st, 2nd, 
..., (n-l)th L.Sq. This gives n(n- 1) blocks of n varieties each. To 
these add the 2 n blocks which form the rows and the columns of the 
square with the numbers from 1 to n 2 . 

It is not known whether a complete set of orthogonal L.Sq’s exists 

when the side n is not a prime power. 

If we start from a PG(2,s) , then we can take as the first two pen¬ 
cils the lines x, = a v and x 2 = a 2 , where a, and a 2 take all the val¬ 
ues of GF(S), and number the rows by a^ and the columns by a 2 . In the 
cell (a,, a 2 ) enter the value a^x+a 2 , where x is a given non-zero ele¬ 
ment of GF(s). The various orthogonal L.Sq’s of the complete set are 
obtained by letting x have the various non-zero values. (This “Bose- 
Stevens” construction was given, independently, in [l25] and in [ 3 ].) 

If n is a prime, and x its primitive root, then 

0 1 ... n— 1 

x i x i + 1 ... x i + n — 1 

x i+ 1 x i+1 + 1 ... x i+1 + n- 1 

x i+n-2 x i+n-2 + l . , , x i+n -2 + „ _ j 

is a L.Sq., and two L.Sq's for different i are orthogonal. The proof is 
immediate. In this case the various L.Sq’s for the different values of 
i = 0, 1, ..., n -2 are obtained by cyclic permutation of the last n - 1 
rows. 

Example , 

Let s = 3. Then we have, for 

x = 1 x = 2 

0 12 0 12 

1 2 0 2 0 1 

2 0 1 12 0 

Conversely, if we number the cells 
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12 3 
4 5 6 


then we obtain a b.i.b.d. (which is 

resolvable) 


123 147 

168 

159 

456 258 

249 

267 

789 369 

357 

348 


In such a complete set any two L.Sq’s differ only through a permu¬ 
tation of their rows (this is called “property D 0 ” in [17]). [46] con¬ 
tains a complete set of L.Sq’s of side 9, with property D 0 , which does 
not arise from a PG(2, 9). 

A collection of n cells in a L.Sq. of side n, such that no row and 
no column contains two cells of the collection, and no two cells of it 
contain the same symbol, is called a transversal of the L.Sq. If the 
same collection of cells is a transversal in every L.Sq. of a set, then 
' we call it a transversal of the set. e 


For instance, if in an orthogonal set we single out one of the L.Sq’s 
and consider the cells which contain the same symbol, then these cells 
form a transversal of the set of the other L.Sq’s. 

We now show that a complete set of orthogonal L.Sq’s cannot have 
a transversal. Start with a cell in the first row. We may assume that 

the first row is the same in all L.Sq’s and that the symbol we start with 
is 0. 


To continue with the transversal, n- 1 cells of the second row are 
available. Take any one of these cells. It must not contain 0, and its 
contents must be different in the n - 1 orthogonal L.Sq’s. Also, it must 
not be the same as the symbol above it in the first row. Clearly, this 
is impossible. Thus there cannot even be two cells which satisfy the 
conditions for a transversal of a complete set of orthogonal L.Sq’s. 

By an application of this result we now prove 
Theorem lll.l If r- 1 L.Sq’s of side n form an orthogonal set, and 
if each L.Sq. contains a subsquare of side r such that all these sub¬ 
squares form a complete orthogonal set between them, then it is impos¬ 
sible to add an rth L.Sq. of side n orthogonal to all the original ones 
unless n> r‘. Moreover, if n is larger than r\ then it must be at least 
r +r (cf. [85]). 


Proof: 
be all 


Without loss of generality, the subsquares may be assumed to 
in the same position. If another L.Sq., orthogonal to all the r- 1 


7 
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given ones, can be added, then the set of the r - 1 L.Sq’s has n trans¬ 
versals, namely the collections of cells which contain all 0, or all l f 
etc. in the added L.Sq. These transversals do not overlap. As we have 
just seen, no two cells of a subsquare can belong to the same transver¬ 
sal, because the r—1 subsquares form a complete orthogonal set; on the 
other hand, every cell of a subsquare lies on some transversal. There¬ 
fore, the number of transversals must be at least as large as the number 
of cells in a subsquare, i.e. n > r 2 . 

If n exceeds r 2 , then there must exist a transversal with no cell 
in a subsquare. Consider the two rectangles within the L.Sq. and out¬ 
side the subsquare which have, respectively, r rows or r columns in 
common with the subsquare. Each contains r cells of a transversal, 
and therefore the remaining n — 2r cells of the transversal must be in 
the subsquare of side n — r, complementary to the original subsquare. 

The elements x u ...,x r of the subsquare must appear in those 
n —2r cells. But if a cell in one of the L.Sq’s contains x t , then the 
corresponding cell in any other L.Sq. cannot contain the same symbol, 
because the pair (x,, x,) is already accounted for within two orthogo¬ 
nal subsquares of side r. Therefore, r— 1 different cells of the trans¬ 
versal must be available to contain x, in the complementary sub-spaces 
of the several L.Sq’s. The same applies to the other x i , and hence 
«-2r>r(r—1), i.e. ra>r 2 +r. 

Two slightly more general theorems are also proved in [85]. 

Incidentally, Parker asserts that there do exist Graeco-Latin 
squares of side 10 without Graeco-Latin subsquares of side 3, and 
such squares were generated on an automatic computer (cf. [82]) . 

Graeco-Latin squares based on groups 

Every row of a L.Sq. contains the same symbols, so that the second, 
third, etc. rows are permutations of the symbols of the first. Denote the 
permutation which carries the symbols of the first row into those of the 
ith by Pi . If i ^ j , then P i Pf x does not. leave any symbol unchanged, 
because if it did, then that symbol would be in the same column after per¬ 
mutation by Pi and also by Pj . This property of 1 f or all pairs of 
different i and j is also sufficient for a L.Sq. to be generated. The P 
do not necessarily form a group, of course. 

Theorem III.2 Let the permutations P,P 1n generate a L.Sq. of 
side n, and let P 21 , P 2n generate another, orthogonal to the first. 
Then the permutations Ptf P 21 , ..., p+ P 2n will also generate a L.Sq. 
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Proof: Pu^ Z i transforms the ith row of the first into the ith row of 
the second, and since each pair occurs precisely once in a Graeco-Latin 
square, the set P H P 2i contains, for every pair x,y, a permutation 
which transforms * into y. There can only be one such transformation. 


The converse of Theorem III.2 is also true: if P H (i = 1, n) 
generate a L.Sq., and P^P 2i 0 = 1> •••> n) generate another, then so 
do the P 2i , and the square generated by the latter is orthogonal to the 
first, because otherwise P’. 1 P 2i would not generate a L.Sq. 

It follows by induction that the L.Sq’s generated by the sets P 
P ri form an orthogonal set if and only if there exist r(r—1) L.Sq’s 
generated by PJ>P ki (/, k = 1 ,... , r ; j ^k). ^ 1 ° 


If the rows in two orthogonal L.Sq’s are the same, though in dif¬ 
ferent order, and if the permutations of the rows form a group, then we 
say that the Graeco-Latin square is based on that group. 

For instance, any two orthogonal L.Sq’s constructed by the method 
of Bose and Stevens (page 26 ), are based on a group. Indeed, let P : 
be the permutation of a w x+a 2 into a u x+a 2 , and P 2 that of a 10 x+a 2 
into a 12 x+a 2 , then P, P 2 permutes the row of symbols a w x+a 2 into 
the row (a n +a 12 — a }0 )x + a 2 . Hence the permutations form a group. 

Also, two orthogonal L.Sq’s +a 2 and a^x 2 + a 2 contain the same 
rows m different order. Clearly the rows a u x } + a 2 and a 12 x 2 + a_ are 

the same, if ai1 x t = a 12 x 2 . 

If two L.Sq’s are based on a group, and if their permutations are, 
respectively, P. and Q A (i = 1, .... «), then P i " 1 Q i = R i} say, are also 
permutations of the same group, and they define a L.Sq. The mapping 
0 the P. on to the P i is a complete mapping. Thus, if a group has no 
complete mapping (e.g. a group of order At + 2), then a Graeco-Latin 
square of that side (may exist, but) cannot be based on that group. 

Conversely, if a complete mapping of a group does exist, then a 

the eC °" Latln SqUare based on that g rou P can be constructed. Represent 
oftifr° UP by 3 regular Permutation group P i , mapped on P i (elements 
k e same regular permutation group). Then, since such a group con- 

unch n ° 6 * ements ' a P art ^ rom the identity, which leaves any symbol 
the j 6d ' Pi ’ Rit and P i R i wil1 a11 generate L.Sq’s, and the first and 
in o ^ are ortho 6 onal - This follows from the converse of Theorem 
2 ’ becau se (/^- , )(P.P.) = P.. 
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Example 


Pi 

Pi 

Pi Pi 

I 

I 

I 

(12) (34) 

(13) (24) 

(14) (23) 

(13) (24) 

(14) (23) 

(12) (34) 

(14) (23) 

(12) (34) 

(13) (24) 


The L.Sq’s 


1234 

1234 

2143 

4321 

3412 

2143 

4321 

3412 


are orthogonal. A third L.Sq., orthogonal to these two, arises when we 
interchange the roles of R i and P i R i . This is here possible, because 


.P:R, = R. 
11 1 


The L.Sq. 1234 completes the set. 

3412 

4321 

2143 


Existence theorems 

We shall now mainly be interested in the existence of Graeco-Latin 
squares, i.e. of pairs of orthogonal L.Sq’s, and start with some prelimi¬ 
nary considerations. 

There are L.Sq’s which cannot possibly be one of an orthogonal 
pair. An example is given by the following- 

Theorem III.3 (see [66]) If in a L.Sq. of side 4t + 2 (t being a posi¬ 
tive integer) not more than t cells of the first 2t+ 1 rows and columns 
contain symbols from a given list of 2 1+ 1 symbols, then the L.Sq. can¬ 
not be part of a Graeco-Latin square. 

Proof : Divide the L.sq. into four squares, each of side 2 1 + 1, and 

denote them as follows: I II 

III IV 

If any symbol occurs in I m times, then it will occur 2t+ 1-m times 
in II as well as in III, and hence again m times in IV. Thus in I and 
IV together (briefly, in I/IV) each of the 4*+ 2 symbols appears an 
even number of times. 

Assume now that there exist two orthogonal L.Sq’s L and L' . 

Let the symbols in L be 0, 1, ..., At +1 and call the symbols 0, ..., 2 1 
A-marks, and the others D- marks. Assume that in I of L there are not 
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more than t cells with a B-mark in them. Then there will not be more 
than 2 1 of the A-marks outside I/IV and therefore at least 2t + 2 sym¬ 
bols of L occur, within I/IV, together with all the A-marks within 
I/IV. 

However, at most 2 1 symbols of L' occur in I/IV together with 
B-marks, and hence at least 2 remaining symbols of L' (out of those 
2t + 2) appear, in I/IV, only with /1-marks and therefore an odd number 
of times. This contradicts the fact, mentioned above, that in any L.Sq. 
every symbol appears an even number of times within I/IV. No orthogo¬ 
nal pair can therefore exist under the circumstances. 

It is clear that we may permute rows, columns, and symbols of L 
and obtain an equivalent theorem. Hence a L.Sq. of side 4?+ 2 which 
contains a L.Sq. of side 2t + 1 cannot be one of a pair of orthogonal 
squares. For instance, a multiplication table of a group of order 4t + 2 
forms a L.Sq. and contains one of side 2t+ 1, which is the multiplica¬ 
tion table of a subgroup. 

Also, if we construct a L.Sq. of side 6 from 012345 by cyclic per¬ 
mutations, then this L.Sq. contains L.Sq’s of side 3, consisting of 
alternate rows and columns. Tarry has, in fact, proved ([ 1261) that no 
Graeco-Latin square of side 6 exists and, as Mann [66] has pointed out, 
the above approach could simplify his proof. 

Ostrowski and Duren have proved ([so]) that the result in Theorem 
III.3 is best, in the sense that if more than t cells of I contain B-marks, 
a L.Sq. orthogonal to that given may exist. They exhibit the following 
examples of side 10: 


01234 

56789 

01923 

84657 

34012 

79865 

67895 

23104 

43120 

97658 

93746 

58210 

12407 

85396 

38254 

79061 

20375 

68941 

14507 

36982 

57698 

34120 

25619 

40873 

89756 

12034 

40138 

62795 

65981 

43207 

56480 

17329 

98563 

01472 

82071 

95436 

76849 

20513 

79362 

01548 


Three cells in the upper left-hand quarter of the first L.Sq. contain 
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R marks and yet the second L.Sq. is orthogonal to the first. 

Mann has also proved in [66] that if in a L.Sq. of side 4« + l the 
square formed by the first 21 rows and columns contains fewer than 
1/2 symbols out of a list of 2t, then this L.Sq. cannot be a member of 

an orthogonal pair. 


* y 


Orthogonal arrays . _ _ , 

For our further investigations we write sets of orthogonal L.Sq s 
in a more convenient way. If the side is n, then we write two rows of 

n 2 elements each, thus: 


0 0 ... 0 1 1 ... 1 • • • 

0 1 ... (n — 1) 0 1 ... (n- 1) . • ■ 0 1 -C"- 1 ) 

The same symbols will also be used for the elements of the L.Sq’s. 

Imagine the rows and the columns numbered from 0 to n- 1. We 
then write into the next row, underneath the symbols \j\, the entry of 
the first L.Sq. in the cell of the ith row and y th column. Each further 
row is allocated to a further L.Sq. of the set, in the same manner. 

Since the L.Sq’s are pairwise orthogonal, the array which we obtain 
has the property that in any two rows all possible combinations | j | 
appear precisely once. The array has n 2 columns and there are in it 
n different symbols. If it has c rows, corresponding to c-2 L.Sq’s, 
then we call it an orthogonal array of size n 2 , c constraints, level n, 
strength 2, and denote it by [n 2 , c,n, 2] . If the orthogonal set is com¬ 
plete, then c = n+ 1. Such an array is equivalent to a resolvable b.i.b.d. 
with v=-. n 2 , b = n(n+ 1), r = n+ 1, k = n, X = 1, which we construct as 
follows. 

Number the columns and the rows of the array. The number of a 
column will appear in the ith set, in that block which is indicated by 
the element of that column in the ith row. 


Example 

Two orthogonal 
orthogonal array 


Latin squares of side 3 are equivalent to the 

0 0 0 1 1 1 2 2 2 
012012012 
0 1 2 2 0 1 1 2 0 
012120201 


This leads to the resolvable b.i.b.d. on page 27. 

More general orthogonal arrays are defined as follows. 
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An orthogonal array of size N, c constraints, level n and strength 
t, written [/V, c, n, i], is a rectangular arrangement of n different sym¬ 
bols in N columns and c rows, and such that in any t rows the N 
columns contain all the n c possible ordered sets of n elements the 
same number of times. If we call this number the index of the array. A, 
then clearly N = \n c . (These definitions are due to C.R. Rao, see 
[90].) 

If we retain only c (^.f) rows from such an array, then we have 
a [N, c, n, f] with the same index. An array of strength t and index 
A is also of strength t' and index A, ( t> *'). If A is a power of 
n, then the array is called a hyper-Latin cube. When A = 1, and t = 2, 
then such a hyper-Latin cube is the same as an orthogonal set of L.Sq’s. 

Raghavarao defines in [87] a Latin cube of the first order, of side 
n, as a cube of n letters n 2 times repeated, each letter being exactly 
n times in each of its planes (there are three parallel sets of these, 
with n planes in each set). Two such cubes are orthogonal if every 
letter of one cube appears n times in the same position with every 
other letter in the other cube. A set of c pairwise orthogonal cubes is 
equivalent to a [n 3 , c, n, 2] . We shall see that the maximum number of 
cubes in a set of pairwise orthogonal ones is n 2 +n- 2 (see page 42). 

A set consisting of this number of cubes is called complete. The con¬ 
struction was studied in [60] and [ 24 ]. 

We give here an example of a complete set of 10 orthogonal Latin 
cubes of side 3 from [ 42 I. 


012 

012 

012 

012 

000 

000 

012 

021 

021 

021 

120 

201 

012 

012 

111 

111 

120 

102 

201 

120 

201 

120 

012 

012 

222 

222 

201 

210 

120 

201 

012 

012 

120 

201 

111 

222 

120 

102 

120 

201 

120 

201 

120 

201 

222 

000 

201 

210 

012 

012 

201 

120 

120 

201 

000 

111 

012 

012 

201 

120 

012 

012 

201 

120 

222 

111 

201 

210 

201 

120 

120 

201 

201 

120 

000 

222 

012 

021 

120 

201 

201 

120 

201 

120 

111 

000 

120 

102 

012 

012 


Construction 

In order to construct an orthogonal array, we may turn to a 
p Q(m, s), following C.R. Rao ([%]). After omitting a selected 
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/ 


Cm-11-flat “at infinity”, we obtain a EC(m, s). The (m l)-flat 

at infinity contains (.--»/(»-D different ( "" ”.'f“ f mm 
through each of them there pass S different ( m _ 1 )-flats apart from 
that at infinity. We call each of the (m-2Hlats at “taty a v* « 
of a parallel pencil of (m_l)-flats, and there will be one (m- 1)- 
flat from each of the (s"*-l)/(s-l) pencils through every point of 
the EO(m,s). Any two (m-D-HaU from different pencils intersect 

in an (m-2)-flat with s m ~ 2 points. 

To the (m-l)-flats of each pencil attach the elements of a 

GF(s), number the pencils in some way, and attach to each of the s m 
points of the EG(m, s) the ordered set of the numbers of the (w - 1)- 
flats passing through it. If we give a column to each point and a row to 
each pencil, then we have s m columns and (s m -l)/(s-l) rows. 
Write into the column of each point the numbers of its (m - l)-flats. 

In any two rows any combination of the elements of the GF(s) appears 
then s m ~ 2 times (this being the number of points in the (m-2)-flat 
in which the two pencils, represented by the rows, intersect), and we 

have thus obtained a [s m , s, 2 ] . It is equivalent to an affine 

resolvable b.i.b.d. with v = s m , b = s(s m — l)/(s — 1), k = s , 
r = (s m - l)/(s-l), A = (s m_1 - l)/(s-l). 

Another construction of arrays of strength 2 is based on Hadamard 
matrices. Let such a matrix of order 4 A be given. It can always be 
such that all elements in the first row are 1. Omitting this row, and 
replacing -1 by 0, we obtain an array [4A, 4A -1, 2, 2]. The proof 
is almost immediate. The resulting arrays are explicitly 


for all A < 25, except for A = 23. 

Also in [86], Plackett and Burman have constructed orthogonal 
arrays by cyclically permuting a row of s m — 1 elements (s m — l)/(s — 1) 
times and adding a column consisting entirely of the element 0. 

More general arrays of strength t can be constructed if we can find 
c pencils in a EG(m, s') such that any t (m-l)-flats from different 
pencils have s 1a ~ t points in common, i.e. such that they intersect in an 
(m-t)-flat. This means that the vertices of the pencils (they are again 
(m —2)-flats in a PG(m, s) in which the EG(m, s) is embedded) 
must be such that any t of them have an (m -1 - 1 )-flat, but no flat of 
higher dimension in common. We then obtain an [s m , c, s, f]. 

The f (m - 1 )-flats of a pencil have analytic representations 
a 0 + a i*i+ ••• + a m*m = 0 where the ai( 1 i = 1, ..., m ) are fixed and the 
value of a Q varies. Let a 0 be the value attached to the (m- l)-flat 
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whose equation contains it, out of a pencil defined by the other a i . 

Let c pencils be given by (a^ ,..., ajm ) (/ = 1,..., c ). Select c (m-1)- 
flats, one from each pencil, denoting them by a j0 . We want any t of 
them to intersect in a (m — O-flat, whatever the a j0 . Hence any sub¬ 
matrix of t rows of the matrix 



must have rank t. Each point (x 1f ..., x m ) in a EG(m,s ) receives 
from the jth pencil the mark , given by 


a j o + a j 1*1 + 


+ a jm X m 


0. 


Therefore the set (a 10 , ..., a t0 ) will be allotted to the solutions of 


a 


10 


+ a t1 x y + 


+ 

1771 m 


0 


a t0 + a i1 X \ + + a tm X m - 

i.e. altogether to s m ~ t points. 

It follows that if a matrix, as defined above, is given, then the s m 
columns 


formed from the s m sets (jc,, ..., x m ) will form an array [s m , c, s, t], 
(Cf. [11].) 

A matrix of the required property is that of the co-ordinates of c 
points in a PG(m,s), where no t points lie in a space of less than 
t— 1 dimensions. 

Example (from [ll]). 

The ten points of the quadric x* + x% = x 3 x 4 in a PG(3, 3) are 

(0010) (0001) (0111) (0122) (1011) 

(1022) (1112) (1121) (1212) (1221) 

No three of these points are on the same line. Therefore the co¬ 
ordinates can be used as the rows of a matrix, from which an array 
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r qi in 3 31 can be constructed. 

[ EstU; Seiden has given in [lod ™les for constructing arrays f 
strength 3 from suitable ones of strength 2, and also other™,se. We 

give here a simple example There will be A columns 

Take any three rows of a L4A, c, x, 

0 1 0 1 

each of which is either 0 or 1, A columns either 0 or 1, A e,ther 

0 1 1 

0 1 10 

1 or 0, and A. either 0 or 1. 

01 01 i u 

Proof: Let x ijk denote the number of columns / in three rows chosen. 

Then 


2 x,, 7 . = lx. jk = 2 


ijk 


For instance, * 0 oo + *oio = *oio + *iio * 


Hjk 


110 


= A. 

+ x m = A, so that 


= A . Similarly, * 0 oi+*iio - *mu .«■ — 

Thus if to a given array we add all columns obtained by changing in 

those'already given 0 into 1 and vice versa, then we obtain a 
[8A, c, 2, 3 ]» and adding a row consisting of 4 A times 0 and 4A 

times 1, we have a [ 8A, c+ 1, 2, 3] • 


010 +*101 ~ *100 +*011 


= A. 


Example oon 00111100 

[4, 3, 2, 2] 0101—■ [8, 4, 2, 3] 01011010 

0110 01101001 

00001111 


Another method of constructing orthogonal arrays will now be ex¬ 
plained. 

There are s b different polynomials 

a t _^x l ~^ + ... + a, x + a Q 

with coefficients from a GF(s). We call them )/. (*), j = 1, 2, s t . 
Denote the elements of the GF(s ) by x,, ..., . Let t^s. Then the 

matrix y. ( x i ), i for rows and j for columns, is an orthogonal array 

IV, s, s, t). 

For the proof, (from [ 27 ]), take any t rows. Without loss of gener¬ 
ality, we may assume that they are the first t. Since the array has s ( 
columns, it is sufficient to show that any two columns of the t rows 
are different. The identity of two columns would mean that 
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c-i 


t-\ 


t -1 


,Io a i x % ~ & a i x i> i■ e • ( a j — a /) x\ 

for i — 1) •••» t. But the Vandermonde determinant 

1 x , 


= 0 


XT ... x 


1 X. 


6-1 

1 


6-1 


is only zero if two of the x i are equal, and since this is not the case 
here, the only solution of the system of equations in the (a^ — ay) is 
a j~ a j = 0 for all y, so that the two polynomials would then be identi¬ 
cal. Thus the s* different polynomials lead to s* different columns. 

The [ s*, s, s, t] which we have obtained can be expanded into a 
[s 4 , s+ 1, s, t] by the addition of a row which contains, in each column, 
the leading coefficient (i.e. a t _ x ) of the polynomial to which that 
column corresponds. For the proof we need only investigate those pairs 
of columns in which the two added elements are equal. Then a t 1 = 

6-2 6-2 

a t-i anc ^ for S a j — .£ Oj'xi = 0 to hold, a Vandermonde deter- 

3 U 3 0 

minant of order t — 1 would have to be zero. We conclude, as before, 
that the array is still orthogonal. 

Example 

s= 4, t= 2. We obtain a [ 16, 5, 4, 2]. 

The elements of GF( 4) are 0,l,y,z=l + y. The multiplication 
table is 


0 


and the array is 

0 0 0 0 
1 

y 

z 

1 


y 

z 

1 

y 


o 

l 

y 

z 

1 

0 

z 

y 

l 


o 

y 

z 

1 


0 

z 

1 

X 


y 

y 

y 

y 

o 


y 

z 

0 

1 

1 


y 

o 

l 

z 

y 


y 

l 

z 

0 

z 


z 

z 

z 

z 

0 


z 

X 

1 

0 

1 


z 

1 

0 

y 

y 


z 

0 

y 

l 


The last line gives the leading coefficients of the respective poly- 
nomials, namely 
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0 x yx zx 1 x+1 yx+1 zx+1 y *+y yx+y zx +y 

z x+z yx+z z x+z. 


If s is a power of 2 and ( = 3, then yet another row can be added, 
by writing into each column the coefficient a t _ 2 of the corresponding 
polynomial. For three rows which do not contain either of the rows of 
a._. and a._ 2 , the argument is the same as before. It may be reduced 
to that argument if we consider the two rows and a third. However, if 


the relevant determinant is either 



1 X, 


i 

X, 2 

either 

or 


2 


1 x 2 


i 

x 2 


The first is not 


teristic 2, ( x, — x 2 ) 2 . 

We now describe a similar method, given in [ll], for the construe 
tion of a [As 2 , As, s, 2], where A and s are powers of the same 


prime, say A = p u , and s = p v . u+v-1 

The elements of a GF(p u + v ) can be represented as x j = .| 0 a. x\ 

with coefficients from GF(p). To each x^, let there correspond the 

V-1 

sum of the first v terms, i.e. y, = 2 a. x\ so that each y. corresponds 

3 i=0 


to p u = A different x .. 

Write down the multiplication table of the x. . Amongst the differ¬ 
ences of corresponding elements in any two rows each element of 
GF(p u+v ) appears once, and if we replace x^ by yj then each differ¬ 
ence appears A times. 

Now replace each element of the modified multiplication table 
(i.e. after replacing Xj by y^ ) by its row in the addition table of the 
y. . The result is a [As 2 , As, s, 2]. For the proof, which is quite 
simple, we refer the reader to [ll]. A further row can be added which 
consists .successively of As times the first, As times the second, 
..., As times the sth element of GF(p v )- 


Example 

p = 2, u = v = 1, A = 2, 
Consider 

x, 0 

and hence 


s = 2. We construct a [8, 5, 2, 2]. 
1 x 1+x 


y, 0 1 0 1 

The multiplication table of the Xj , i.e. 
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0 0 0 0 

0 1 x 1 + jc 

0 x 1 + x 1 

0 1+x 1 x 

is then “modified’ * into 

0 0 0 0 
0 10 1 
0 0 11 
0 111 

and the array becomes 

01010101 

01100110 

01011010 

01101001 

to which can be added 00001111. 

If u>v, then at least one more row can now be added. We refer to 
[ll] and merely mention that altogether c = 1 + A( s w+ '— l)/(s w —s w ~') 
rows can be obtained, where w is the largest integer contained in u/v. 
In our example, the two further rows of the array are 

00110011 

00111100. 


Maximal strength 

It is known (see [l29]) that m t , that is the largest number of points 
in a PG(m - 1, s),such that no t of them lie in a flat of dimension small¬ 
er than t — 1, satisfies the inequalities 

s ™ > . (j*) (s-lf when t = 2 h 

and 

* Jo C‘) (S - 1)4 + ("V X ) when * = 2h + 1. 

We shall now show that analogous inequalities hold for the largest 
number of rows in any orthogonal array, whether derived from a EG(m,s) 
or not. We have to replace s m by N. 

Suppose N, s, and t to be given. We want to find the largest pos¬ 
sible value of c, say /(/V, s, t). 
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Let a square matrix of order s be given, say («,,), i, I - 0, .... 

s _l, and consider the products F( I,, .... i.) - (^0°^) 

f S l a ■ y cj V These are sums of s° terms ... <* iojg ••• y C j B > 
each term being determined by the ordered set . i 0 , called an 

““insider also an orthogonal array [IV, c, s i] and retain in 
p U i) only those terms which correspond to assemblies given 

by the*columns of the array. We call the resulting expression G( , 
i ). Take two of these, say G( u„ .... u c ) an ( v,, .... v unc 

tions of the products y,. ...y^ . Two of them are orthogonal if 

S = = °’ 

where the summation extends over all assemblies from the given ortho- 

gona^array tr^y ^ ^ strength c and hence contains all possible s* 

assemblies A times, then 

S = A(a Ui ,a Vi , + ... + a u 0 c av c c ) • 

If the matrix ( a tj ) is orthogonal, then we see that 5=0, which 
ever {u,.u c } and \ v,.v c i were chosen, unless - v,, ..., 

C The case which interests us more is that where t<c. Therefore we 
assume that a 00 - ... - a 0s _, = 1. and also that both the sets M and 
{v\ contain 0 in i v=c-t common positions. Then the w a 
each term of S, which refer to those common positions, are equal o , 

and we have v t n 

S = A (a Ui ,a Vi , + ... + a u o c a v c c> ~ 

This remains true if w = c-i‘> c-t, because an array of strength I 
is also of strength ('(<»)■ We have thus shown . 

Theorem III.4 Any two expressions G{u v .... u 0 ) and G{ v v ..., c > 
where the sets I u I and 1 v I differ in not more than t arguments, while 
the other arguments are 0 in both expressions, are orthogonal. 

We turn now to the question of how to find a set of functions 
G(q, ..., i a ) which are, by virtue of Theorem III.4, pairwise orthogona • 

We distinguish two cases. 

Case (i): even t (= 2 h). 
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We take all those G(i,, ..., i o ) where not more than h arguments 
differ from 0. Their number is 

1 (all arguments equal to 0) 

+ c ( s — 1) (all but one equal to 0) 

+ ( 2 ) ( s ~ 1 ) 2 (all but two equal to 0) 

+ K) u-d*. 

Because the GO’,, i c ) are linear combinations of terms 
yy •••y C j c » total number of pairwise orthogonal GO,, i c ) can¬ 
not exceed the number of assemblies in the orthogonal array, so that 

■t/2 /r\ 

Case (ii): odd t (= 2k + 1). 

We take all those GO,, ..., i c ) where not more than h arguments 
differ from 0, as in (i), and to this set we add those GO,, ..., i c ) 
where t,, say, differs from 0, and so do h others of the arguments. 
There will still be, in any two G, not more than t arguments in which 
they differ, while the other arguments are 0 in either. The number of 

the additional functions G is ( C ^ X ) ( s-l) h +\ and hence 

N ^ Mo (?) + { C h X ) ( S ~1) M1 , where h = %(j-l). 

Example 

[4, 3, 2, 2] orthogonal matrix /1 1\ 

\1 -1/ 


v-oemcients 

of 

>oo >10 >20 

>00 >11 

>21 

>01 >10 >21 

>01 >11 >20 

G(0 

0 

0) 

0 

0 

& 

0 

0 

» 

0 

0 

0 

0 

0 

cf 

a oi 

O 

O 

O 

& 

O 

a oi 

a oi a oo 

G( 0 

0 

1 ) 

0 

0 

» 

0 

0 

0 

Q 

0 

0 

a 

0 

fl 11 

a oi a oo fl 11 

a 01 

a oi a io 

G( 0 

1 

0) 

a oo a io a oo 

0 

0 

a oi 

a 

0 

0 

0 

0 

a oi 

a i 1 a oo 

G (1 

0 

0) 

0 

0 

0 

& 

0 

0 

0 

G 

0 

a 0i 

a l 1 a 00 a 01 

a 11 

a oi a oo 


Compare G( 0 0 0) and G(0 0 1): 
a oo a io + a 01 a,i + a 01 a n + a Q0 a, 0 = 2(a Q0 a 10 + a 01 a t1 ) = 0. 
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Compare G(0 0 1) and <3(0 1 0): 


( a oo a io ) ( a io a oo 


) + ( 


+ ( a oi a io ) ( a n a oo ) 


) ( a ii a 01 ) + ( a oo a ii ) < a io a oi ) + 
( a oo a io + a oi a ii ^ = 


The inequality above gives 4 » 1 + 3(2-1) and we see that, in 
fact, equality holds in this case. However, this is not always so, and 
for special cases closer upper bounds have been established, as we 

= 2 we obtain f(N, s, 2) ^ (/V-l)/(s-l) from the above 
formula, and this inequality was also proved, independently, in |86J 
and in [l08]. There exist N and s for which the upper limit is reached 
We have seen, for instance, that arrays [s m , (s m - l)/(s- 1), s, 2J 
can be constructed from a PG(m, s). For a set of orthogonal Latin 
cubes of the first order of side s we have N = s 3 , and since 
f(s 3 , S , 2) = s 2 + s+1, a complete set consists of s 2 + s - 2 cubes as 

stated above (page 33). 

For t = 3, we obtain /(As 3 , s, 3) ^ (As 2 - l)/(s-l) + 1. For 
s = 2 this reduces to 4A, and we have seen earlier how to construct an 
[ 8A, c+ 1, 2, 3] from a [4A, c, 2, 2]. We find now that if the latter 
has the largest possible number of rows (c = 4A-l), then the former 

has also the largest possible number (c+ 1 = 4A). 

In [27] K.A. Bush has found expressions for f(s t , s, t) by build¬ 
ing up an array step by step. We merely quote here his results. He has 
established the following upper bounds for /(s*, s, t): 


s < t 
t + 1 

(equality, if s 
is a prime power) 


s 

> t 

s even 

s odd, and 


t ^ 3 

s+t- 1 

s+t — 2 


We have seen how to construct an array |V, s+ 1, s, t ] from poly¬ 
nomials in GF(s ) when s > t ; in that case s must be the power of a 
prime. Hence then /( s 2 , s, 2) = s+ 1, and /( s 3 , s, 3) = s+ 1 if s is 
odd. We have also seen that, if s is even, and t = 3, then /( s 4 , s, t ) 
> s+2. From the table above it follows that equality then holds in the 
last case. 

Bose and Bush, in [ll], have improved on /(As 3 , s, 3) < 

(As 2 - l)/(s- 1) + 1 when A—1 is not divisible by s —1. Call the 
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remainder b. Then the right-hand side can be reduced by the largest 
integer contained in t+ 1, where t is the (always existing) positive 
root of t 2 -t(2b — 2s + l)+ b(b-s+l) = 0. 

^ A. — 1= a(s — 1) and (s — l) 2 (s — 2) is not divisible by as + 2, 
then /(As , s, 3) ^ (As — l)/(s — 1) — 1. This is also shown in 

[ill- 

E. Seiden has proved in [ 107 ] and in [ios] that /( 81, 3, 3) = 10. 
The construction of the corresponding array was given earlier in this 
chapter (pages 35-36). 

The following formula follows from the structure of an orthogonal 
array [N, c, s, t]. Let n^ be the number of columns, not counting the 
first, in which there are i coincidences, i.e. i values which equal the 
value in the same row of the first column. Then 

Jo (0 n i - 00 Us'-' 1 - 1), (A = 0, 1.()• 


Proof : Take any h t ) rows of the array. This can be done in 

ways. Within h given rows, A s* h — 1 other columns will be identical 
with the first. Amongst these, any column with i coincidences will be 

counted (J^ times among the columns considered. The number of 
the latter is the right-hand side in the formula to be proved. (Note that 
\/iJ = 0 when i < h .) 

For h = 0 we obtain 2 = N — 1, which is evident from first 

principles 


MacNeish’s theorem 

In [28] K. A. Bush has shown how to construct an orthogonal 
array from given ones of equal strength. Let an [N v c,, s y , f] and an 
[/V 2 , c 2 , s 2 , t] be given, and denote the smaller of c, and c 2 by c. 

Let the first array be (a i; ) and the second (b^). We take as elements 
of the new array the s t s 2 pairs ( a ef , b gh ) and construct the array: 


^ a n )•.. (a Wi , b u ) (a,,, b 12 )... (a w ^, b 12 )... (a 





(a .„ b c , )...(a Wi , i>„ 2 )...(a oW) , b c2 )...(a o ^, b ). 

itwitVT“ ly See ” th3t thiS U a ” atray lN ' N *’ C ’ s iV Combining 
th other arrays, we can obtain a [II/V*, c, Es*, t], where c = mi 


min c A 
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f ♦ .^tinned in \l ]. From it we obtain 
The result was first mentioned in l/j 

_ r + m n a i n ar Then an orthogonal array 

Theorem III.5 Let N = P, -P r • , , 

[iV 2 , c, N, 2] exists, where c is the smallest of p. 

. _ r rt 2o t + 1. d°*. 2 , i.e. complete 

This is so, because arrays yp i > P i J 

sets of orthogonal L.Sq's do exist whenever the side is a prime power. 


(Cf 'u“im 2.5 was proved by MacNeish in [62]. He conjectured 
that he had, in fact, proved that c-2, as defined above was the 
largest possible number of pairwise orthogonal L.Sq s. His proo o 
the latter proposition was shown in [61] to be erroneous, though the 
proposition itself might still have been correct. (See, however, the 

next section.) 


Orthogonal Latin squares 

We shall now concentrate on orthogonal arrays of strength 2 and 
index 1, i.e. on sets of pairwise orthogonal L.Sq’s. To begin with, we 
consider the consequences of Theorem HI.5 and of the conjecture of 
MacNeish. 

If N is odd, then Theorem III.5 guarantees the existence of at 
least two orthogonal L.Sq’s (a Graeco-Latin square), because then the 
smallest prime factor cannot be smaller than 3. If N is even and con¬ 
tains the factor 4, then this is also true. But if N is of the form 
2(2 1+ 1) = 4r+2, then no Graeco-Latin square would exist if min(p^-l) 
were indeed the largest possible number in a set. Euler conjectured 
that no Graeco-Latin square of side At + 2 existed, and this is obvi¬ 
ously true for t = 0, and was proved in [l26] for t = 1 by a method 
of exhaustive enumeration. However, we shall see that it is incorrect 

for all other values of t. 

We shall first describe the breakthrough achieved by Parker in [8lJ- 

Let an orthogonal array [ k 2 , c, k, 2] be given. If a b.i.b.d., D 
say, with parameters v, b, k, r and with A = 1 exists, then an orthogo¬ 
nal array [ v 2 , c- 1, v, 2] can be constructed as follows. 

The [k 2 , c, k, 2] is equivalent to a set of c-2 orthogonal L.Sq’s 
of side k. We can write the elements, and hence the array, in such a 
way that in the latter the first k columns and the first 2 rows are as 
shown in the table on the following page. 

Omitting the first row, and the first k columns, we obtain an array, 
A:I , say, such that in any two rows all k(k — 1) ordered combinations 
without repetition of elements appear precisely once. (C.R. Rao denotes 
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0 o ... o 

o 1 ... (*-i) 

0 1 ... (fc — 1 ) 

0 1 ... (k-1) 

such an array by [ k(k — 1), c — 1, k, 2] : I; see [ 92 ], where a different 
type of modified array is also defined.) This array can be divided into 
(/c — 1) subarrays such that in every row of each subarray every symbol 
appears precisely once. 

We take one block of D, and replace in A :/ the element x (a 
number) by that in the xth position of that block of D. We do this for 
all blocks and place the resulting arrays next to each other. Thus we 
obtain a matrix of bk(k-l) columns and c —1 rows. There appear in 
it altogether v symbols, and in any two rows any ordered pair 
(y ^ z ) appears precisely once, because y and z appear once together 
in the same block, say in its u,th and u 2 th position, and thus they 

1J 

were inserted in the final design where 1 appeared in A : / when 

u 2 

that block was considered. Since A = 1, and therefore bk(k — 1) = 
vr(k — 1) = v(v — 1), we have obtained a [ v(v— 1), c — 1, v, 2] :/. 

We call it B :/. 

By adding a matrix of c — 1 rows, all of them listing the v varie¬ 
ties of the b.i.b.d., we obtain finally a [v 2 , c— 1, v, 2]. 

We can add one more row to this array, obtaining a [v 2 , c, v, 2], 
if B:I can be subdivided into v-1 subarrays of v columns each, such 
that every subarray contains all v varieties precisely once in each row. 
This is the case if the b.i.b.d. D was resolvable or symmetric, assum¬ 
ing that in the latter case it was written as a Latin rectangle (Youden 
square). In these cases we can construct, from the [ v(v — 1), c 1, v, 2] 
a [v 2 , c, v, 2] in the same way as the [A 2 , c, k, 2] given at the 
beginning of this section can be reconstructed from its truncated part, 
the A: I. We have proved 

theorem III.6 If a b.i.b.d. with parameters v, b, k, r, and with A = 1 
exists and if an orthogonal array [k z , c, k, 2}, i.e. a set of c - 2 ortho¬ 
gonal L.Sq’s of side k also exists, then a [v 2 , c-1, v, 2], i.e. a set of 
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c-3 orthogonal L.Sq’s of side v can be constructed. 

If the b.i.b.d, is symmetric, or resolvable, then a [ v , c, v, 2], i, e 
a set of c-2 orthogonal L.Sq’s of side v (as many as there were of 
side k ) can be constructed. 

We are, in particular, interested in those cases where c — 2 is 
larger than p a - 1, when p° is the highest prime power factor contained 

in v. 

Let s be a prime power, v = 1+s + s 2 , k = s+ 1. Then the b.i.b.d. 
whose existence is postulated in Theorem III.6 exists as a PG(2,s), 
and is symmetric. If s+1 is also a prime power, then a [k , k+ 1, k, 2] 
exists and hence also a [v 2 , k+ 1, v, 2], i.e. k— 1 orthogonal L.Sq’s 
of side v. 

Both s and s+1 are prime powers when s is a Mersenne prime 
(2 P — 1, where p is a prime), or when s + 1 is a Fermat prime (2 2 +1). 
If s > 3, then in either case s = 1 (mod 3), and hence l + s + s 2 = 3 
(mod 9). MacNeish’s conjecture leads to two orthogonal L.Sq’s of 
side 1+s+s 2 , while Theorem III.6 gives s. 

If i +s + 5 2 is itself a prime power, then a set of as many as s + s 
orthogonal L.Sq’s exists, which is more than s. The smallest value 
which gives something new is s = 4, when 1+s + s 2 = 21, not a prime 
power. A set of 4 orthogonal L.Sq’s of side 21 is exhibited in [si]. 

It is the first counter-example to the conjecture of MacNeish. 

We have seen (page 10) that a resolvable b.i.b.d. with v = > 

k = y 2 s, A = 1 exists if s = 2 m . Let m = 2 n, then v=2 Zn_1 (2 n +l)(2 7l -l), 
k = 2 2n_1 , so that 2 Zn_1 -l orthogonal L.Sq’s of side v exist, while 
MacNeish’s conjecture gives only 2 n -2 ([2l]). This is true for all int¬ 
eger values of n. 

We have also seen (page 11) that a resolvable b.i.b.d. with v = 
s 3 + 1, k= s+1, A = 1 exists if s is a prime power. Hence there exist 
as many orthogonal L.Sq's with side s 3 + 1 as there are for side s+ 1. 

Example (from [2ll) 

s = 31, k = 32 = 2 5 , v = 29792 = 19 x 32 x 49. There are at 
least 31 orthogonal L.Sq’s of side v, which is more than 18 (= 19 — 1)- 

No counter-example to Euler’s conjecture has been produced in this 
way. To prove that this conjecture is false, further investigation is 
required. (The counter-example on p.31 was constructed only later.) 

For a design which is balanced, i.e. in which each pair of varieties 
appears precisely once within the same block, and which contains blocks 
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of sizes k y , k m out of v varieties altogether, we write (v; ky ,..., k^. 

It has been generalised, but we shall not deal with the more general 
cases. 

The falsity of Euler’s conjecture 

If there exist arrays [k 2 , c if k it 2] for each i, then an array 
[ v 2 , c — 1, v, 2] can be obtained, where c is the smallest of the c 4 . 

For the purpose of its construction, we treat the [k 2 , k it 2] as we 
treated the arrays of k 2 columns before, put the resulting arrays next 
to each other, and retain only c — 1 rows. It is easily seen that we thus 
obtain a [ v( v -1), c - 1, v, 2] :/ which we can expand into a 
[v 2 , c — 1, v, 2]. 

It is shown in [ 23 ] that we can, in fact, obtain c — 1 rows where c 
is the smallest of Cy +1, ..., c r +1 , c r+1 , c m , if the blocks of sizes 

k y , ..., k T are such that no two of them, of any of these sizes, have any 
variety in common. In this case we use, instead of a [k i (k i — 1), q — 1, 
k it 2] :/ (i = 1, ..., r) the orthogonal array [k 2 , c v k it 2], and the 
matrix of c — 1 rows which we added in the earlier case will now only 
contain those elements, if any, which do not already appear in the blocks 
of sizes ky, ..., k r . 

Example 

(v;ky,k 2 ) = (6, 3, 2). The blocks are given by the columns: 

0 1 2 4 3 5 0 

1 2 3 5 4 1 2 

3 4 5 0 

Any two blocks of size 2 are such that they have no variety in common. 

There exist arrays [9, 4, 3, 2] and [4, 3, 2, 2], and they give 
012012 0011 

120201 and 0101 respectively 
201120 0110 

From the latter arrays we obtain 


013013 

124124 

235235 

450450 

3344 

5551 

0022 

130301 

241412 

352523 

504045 

3434 

5111 

0202 

301130 

412241 

523352 

045504 

3443 

5155 

0220 


No further columns have to be added, because the blocks of size 
two already contain all six variables between them. (This example 
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merely illustrates the procedure. A [ 36, 3, 6, 2] is equivalent to one 
single L.Sq. and, in fact, no pair orthogonal to it easts.) 

Yet one more row can be added if the design from which we start 
is separable (a concept introduced in [ 21 ]), i-e- lf “S' set of blocks of 
equal size can be divided into subsets such that each variety occurs in 
a subset the same number of times. If this number is unity, then we 
have a generalisation of a resolvable b.i.b.d.; if it is larger than one, 
then that of a symmetric b.i.b.d. 

The proof that a further row can be added is analogous to that in 
the earlier case. 

Balanced designs of the type we have been considering can be 
obtained by adding or by omitting varieties. First, take a case of addi¬ 
tion. 

Consider the resolvable b.i.b.d. with v = 15, b = 35, r = 7, k = 3, 

A = 1: 

123 246 437 365 514 

4 12 8 3 9 13 6 14 10 7 11 12 2 10 11 

7 9 14 5 14 11 1 11 8 2 8 13 6 12 9 

6 11 13 7 8 10 5 13 12 1 10 9 3 14 8 

5 10 15 1 12 15 2 9 15 4 14 15 7 13 15 

6 7 1 7 5 2 

5 8 9 1 13 14 

4 13 10 3 10 12 

2 12 14 4 9 11 

3 11 15 6 8 15 

Add to each block of the i th set of 5 blocks (each set contains all 
varieties) a new variety x i , and add a further block (*,, x 2 , • ••> x 7 )- 
This produces a (22; 7, 4). There exist arrays [ 16, 5, 4, 2] and 
[49, 8, 7, 2]. Hence an orthogonal array [22 2 , 4, 22, 2] can be found. 

This is the first known counter-example to Euler’s conjecture. The 
Graeco-Latin square is given in [ 20 ] and also in [ 21 ] - 

The existence of an infinite series of Graeco-Latin squares of 
2 sides 36m+ 22 is proved in [ 21 J, as follows. 

It is known (see Chapter X of [ 98 ]) that resolvable b.i.b.d’s exist 
with parameters v = 6t+ 3, b = ( 2t + 1) (3 1 + 1), r = 3t + 1, k = 3, A = 1- 
If t = 4m + 2, then we have v = 24m + 15, b = (8m + 5) (12m + 7), r = 

12m + 7, k = 3, A = 1. 
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To each block of the ith replication, i.e. the ith set of 8m + 5 blocks, 
add a new variety, and combine all the new varieties in an additional 
block. This produces a (36m+ 22; 12m + 7, 4). 

There exist arrays [16, 5, 4, 2] and [(12m + 7) 2 , 6, 12m+ 7, 2] 
because the smallest possible prime factor of 12m+ 7 is 5. Hence at 
least two orthogonal L.Sq’s of side 36m + 22 exist. 

Yet another counter-example to Euler’s conjecture can be exhibited 
by starting from a resolvable b.i.b.d. with v= s m , b= s m ~\s m - l)/(s- 1), 
r - (s m — l)/(s — 1), k = s, A = 1, where s is a prime power. These 
b.i.b.d’s can be constructed from a EG(m,s) by taking lines as blocks. 
Add another variable to the blocks of the first (and no other) replication, 
to obtain a ( s m + 1; s +1, s). If s as well as s + 1 are prime powers, then 
an [(s m + l) 2 , s, s m + 1, 2] can be found. 

If, in particular, s is 2 e — 1 and a prime number larger than 3, and 
m = 2n , then s m + 1 = (2 s — l) 2n + 1 = 2(mod 4), so that once more at least 
two orthogonal L.Sq’s exist of a side of the form At+ 2. 

We now give an example of the construction of a balanced design by 
omitting varieties. We start from the b.i.b.d. which is equivalent to a 
PG(2, s ). Let s = 16, and omit six treatments from all the blocks where 
they appear. This gives a (s 2 + s+ 1 -6; s+ 1, s, s- 5), i.e. a (267; 17, 
16, 11). Therefore a set of at least 9 orthogonal L.Sq’s of side 267 
exists, while MacNeish’s conjecture gives only two, since 267 = 3x89. 

We mention one more scheme, and because we shall make use of it 
later, we express it as a theorem. 

Theorem 111.7 Denote the largest number of L.Sq’s of size n in an 
orthogonal set by N(n). Then A/((c-l)m + x) is at least as large as the 
smallest of 

N(c)-l, N(c-1)-1, N(x), and A/(m), 

provided x < m. 

Proof : Let an [m 2 , c, m, 2] be given and omit the first row. The 
remaining array can be divided into m subarrays such that in every row 
of each subarray each symbol (let them be 0, 1, ..., m - 1) occurs pre¬ 
cisely once. Replace the symbol i in the /th row by a new variety 
*(i-i )m+i ‘ Take co l umn s as blocks of a design. This is a configura¬ 
tion with v = (c — 1 )m, b = m 2 , k = c — 1, r = m, but not balanced. (In 
fact, it is a partially balanced incomplete block design of group-divisible 
fype; for these, see Chapter V.) To make it balanced, we add m 
further blocks, each consisting of 




so 


Viiffl+o’ x u-v™+ m ~' ' 


We illustrate this procedure with 
row of a [9, 4,3, 2] we obtain 

m = 3, 

£7=4 

012 

012 

012 

012 

120 

201 

012 

201 

120 

and from it the design 

036 

048 

057 

147 

156 

138 

258 

237 

246 

To this we add the three blocks 



012 

345 

678. 


= 1 , m. 

After omitting the first 


To obtain the balanced design we want, we add a new variety to 
each of the first x (<m) replications and combine all these new varie¬ 
ties in a further block. The result is a (( c-l)m + x ; c,c- 1, x, m). 
Notice that no two of the added m + 1 blocks contain common varieties. 

From this the Theorem follows. 


Example 

m = U, c = 9, x = 7. 

N(9S)>mrn(N(9)-l,m-hNO),N(m - nun(7,6, 6, 10) - 6. 

In [83] E.T. Parker gave the following example of a Graeco-Latin 
square of side 10: Q4m98365 0786935412 

8152739406 6178094523 

9826374510 5027819634 

5983047621 9613782045 

7698415032 3902478156 

6709852143 8491357260 

3071986254 7859246301 

1234560789 4560123789 

2345601897 1234560978 

4560123978 2345601897 

The last three rows and columns contain two orthogonal L.Sq s of 
side 3, and because 10 > 3 2 , but not ^ 3 2 + 3, no third L.Sq. exists 
which is orthogonal to the two squares given. 
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The same Graeco-Latin square is reproduced, with the rows in a 
different order, in [23]- The latter paper contains a somewhat simpler 
argument for its construction than Parker’s original demonstration, and 
we reproduce it here. 

To begin with, form a [(3m + l) 2 , 4, 3m +1, 2] for odd m, as 
follows. We define the matrix A 0 as 


0 

0 

... 0 

1 

2 

... m 

2m 2m — 

1 ... m+1 

*1 

X 2 

•” x m 

1 

2 

... m 

0 

0 

... 0 

x \ x z 

••• x m 

2m 

2m-1 

... m+1 

2m 

2m — 1 

... m+1 

*1 

X 2 

- X m 

0 0 

... 0 

1 

2 

... m 

*1 

*2 

••• X m 

2m 

2m — 1 

... m + 1 

1 2 

... m 

0 

0 

... 0 

where the 

Xi are 

indeterminate symbols. 






Form Aj by adding / (0, 1, ..., 2m) to the elements of A 0 , except 
to the x- which remain unaltered. Reduce mod 2m+ 1. In the matrix 


(i4 0 , A y , ..., A Zm ) any ordered pair of two distinct residuals of 2m + 1, 
or any ordered pair consisting of one of the residuals and one of the x if 
occurs once in any two-rowed submatrix. To obtain an orthogonal array, 
we add the matrix consisting of four identical rows 0 1 ... 2m, and 

also an orthogonal array [m 2 , 4, m, 2], which exists because m is odd. 

' * 

In the latter array we use the x i as elements. 

The number of all the columns we shall then have is 4m (2m + 1) + 
(2m +1) + m 2 = (3m+l) 2 , and they form an orthogonal array. When 
m = An + 3, then this array can be written [(12n + 10) 2 , 4, 12n + 10, 2]. 

We have thus obtained an infinite series of Graeco-Latin squares of 
sides At + 2, and 10 is one of these numbers. 

The complete theorem on Graeco-Latin squares 

We turn now, finally, to the proof that Graeco-Latin squares exist 
for all sides except 2 and 6. The proof starts with listing various 
cases, and we give only the gist of the argument. 

[20] and [ 23 ] contain lists of lower bounds (not necessarily the 
best) of possible numbers of orthogonal L.Sq’s in a set which have been 
proved for various values of n, using methods outlined above. These 
numbers are, in all cases, at least 2, and cover all sides larger than 
6, up to 726 inclusive. 

All numbers of the form At+ 2 which are larger than 726 can be 
written 4a+144b+10, where 0 4 ; a 4 . 35 and b > 5. In Theorem III.7, 
take c = 5, m = 36b, and x = Aa+ 10. Then m ^ 180 and x £ 150, so 
that x < m. It follows that 
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N (4x36b + 4« + 10) > min[/V(5)-l,IV(4)-l,W(4a + 10),iV(36i)]. 

Now N(5)-l = 3, JV(4)-1 = 2 , W(4a + 10) » 2, and N(36fc) > 3, 
because the smallest prime factor of 366 is 4, in view of 6 > S. 
Hence N(4a+1446 + 10) 2, which proves 


Theorem III.8 : Graeco-Latin squares exist for all sides except for 
side 2 and for side 6. 


In [84] Parker mentions that he has found many Graeco-Latin 
squares of side 10 by trials on a computer, and he shows one “totally 
unlike the examples obtained in 1959 by purely mathematical methods". 
But he did not find any set of three orthogonal L.Sq’s of that side. 

On the other hand, N.S. Mendelsohn and others have found four 
sets of five orthogonal L.Sq’s of side 12, but no larger set of this side 
(see [67])- 

In [ 32 ] the authors have proved by methods of number theory that 
the maximal number of pairwise orthogonal L.Sq’s of side n tends, with 
increasing n, to infinity. There exists a number n Q such that N(n) 
n 1/9 ' 

exceeds - for all n > n 0 , but they remark that 1/91 is far from 


being the best exponent. 

To end this chapter, we mention as a matter of interest that Rojas 
and White discuss in [ 97 ] designs of a.b varieties in a.b columns and 
b rows, such that each variety occurs in every row, and once in every 
group of a columns. They call these designs “modified L.Sq’s", and 
they are the same as Yates’s “semi-Latin squares" (see [l3o])- We 
give here an example: 


1 2 3 
4 5 6 
7 8 9 


4 5 6 
7 8 9 
1 2 3 


7 8 9 
1 2 3 
4 5 6. 


EXERCISES 


(1) Construct two orthogonal Latin squares using the following points 
and lines of a PG (2,3): 

Points: A (100) B (010) C (001) a (211) b (121) c (112) 

Q (120) r (102) S (012) X (Oil) Y (101) Z (110) 
P dll) 
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Lines: 


*0 + X 2 
acRB 

= o, 

x 0 + 2x 2 
bPYB 

= o, 

x 2 =0 

QZAB 

(through B) 

X 0 + *1 
abQC 

= o. 

x 0 + 2x, 
cPZC 

= o, 

X, =0 

RYAC 

(through C) 

x t + x 2 

AbcS (1) 

= o, 

x, + 2 x 2 
aYZS (0) 

= 0, 

XQ + x^Xg = 0 

PQRS (2) 

(through S) 

x 0 + jc, + 2x 2 
QYcX (2) 

= 0, 

x 0 + 2x, + x 2 
RZbX (1) 

= 0, 

x, + 2 x 2 = 0 

APaX (0) 

(through X) 


(2) Construct a complete set of orthogonal Latin squares of side 4. 

(3) Construct a complete set of orthogonal Latin squares of side 5. 

(4) Extend the following rectangles into Latin squares of side 5: 



012 


412 

(i) 

340 

(ii> 

341 


123 


123 
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CHAPTER IV 

PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS 


Definition 

In this chapter we deal with a generalisation of b.i.b.d s. Let v 

varieties be given. We define an association scheme as follows: any 

two varieties are either first, or second..., or mth associates, and 
this relationship is symmetric. Each variety has different ith asso¬ 
ciates. Given any two varieties that are ith associates, there are 
p i k = p*. varieties which are ;th associates of the one and kth asso¬ 
ciates of the other; these values do not depend on which pair of ith asso¬ 
ciates we started from. (The cases i = j or i = k are not excluded 
from this definition.) 

We simplify our formulae by taking each variety to be the 0-th 
associate of itself, and only of itself. Hence n 0 = 1, pj£ = fij when 
j = k, and = 0 otherwise, while P^- 1 when i = k, and = 0 otherwise. 

We define a partially balanced incomplete block design (p.b.i.b.d.) 
to be a b.i.b.d., except that not all pairs of varieties appear equally 
often within the same block. Instead, any two varieties which are ith 
associates appear together in A^ blocks. The A^ need not all be dif¬ 
ferent, but if they are all equal (or if m = 1), then the design reduces 
to a b.i.b.d. We also use A 0 as an alternative to r, the number of 


repetitions of any variety. 

These designs were introduced in [ 16]. For a broadening and re¬ 
phrasing of the original definitions see [73] and [19]. 

The p i and n A define the association scheme, but they do not 
define the design by themselves. The parameters v, b, r, k, X i and 


are called of the first kind, and the p\ parameters of the second kind. 
There are altogether 2m + 4 of the former (not counting n 0 ), and 


Vi(m + l) 2 (m- 1 - 2 ) of the latter. 


We have n 0 + n, + ... + n m = v and n 0 A 0 + ... + n m \ m = rk, 
because whenever a variety appears (which is r times), it appears 
together with k - 1 others in the same block, A i times with each of 
its ith associates, and n 0 A 0 = r. 
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Further relations between the parameters can be derived as fol¬ 
lows:— 

Take any two varieties, say v, and v 2 , which are ith associates. 
The /th (*#/) associates of v, number rij, and any of these are either a 
first, or second, ..., or mth associate of v 2 . Hence p*. + ... + p i . = n . 

This formula holds also for i = 0, when it reduces to p?. =™h jt or 

for / = 0, when it means p\ = n a = 1. 33 

0i U 

We shall now prove that m pj k = n, p( k = n k p$ . Consider a variety, 
say v,. It has n* different ith associates, say x u x 2 ,..., Xn. and nj 
different jth associates, say y,, y 2 , ..., y . 

Every one of the x’s is related to p* k of the y’s which are its fcth 
associates. We thus obtain mpf k “connections”. But the relationship 
of being fcth associates is symmetric. Therefore, the same connections 
are obtained if we start from the y’s and relate each of these to the 
x s which are its /cth associates. There are then ^jP^ k connections, 
and their equality establishes the formulae to be proved. 

These relationships establish constraints on the choice of the 
parameter values. Altogether, given the n it only m(m 2 - l)/6 of the 
Pjk are independent. (In this count we exclude 0-th associates, since 
their parameters are all fixed in advance.) 

Proof: To begin with, consider the matrix (pj k ) for some fixed i, 
say for i = 1. Since the matrix is symmetric, it is defined by l m + 1 ) 

values; since the row and column totals are fixed, only ( m ) of these 
are independent. ' ^ ' 

If the p\ k are known, then the relation = n 2 p? k determines 
al p ik = Ply Hence in the matrix (pf k ) only the second to mth 
columns and rows remain free to be determined, and these form a mat¬ 
rix in which ( 2 J values are independent, since row and column 
totals are again fixed . 

Continuing, we find {"'f] independent values among the p.*, 
and so on, until ( m = 1 value remains among the p™~>. The 


values of the p» are then implicitly fixed and we have in all "f 2 ( m ~ ‘= 

ITU^ 6 inde P endent Values of the parameters of the second kind, 
v he last formula is easily verified by induction.) 

This completes the proof. 

n .„, F ° r m = 2, we have m(m2 ~ 1)76 ■ and ‘his will be used in the 
next chapter. 
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motors satisfy the necessary conditions which 
Even when the param p.b.i.b.d. with such para- 

we have established, it is ^ 13 ! ,-4. *-4, A,- 2, 

raet ers exists. For instance, let 

A 2 =0,n, = 6,n 2 = 6 - / 3 3 \ 

Gjl> - ( 33 )' W = M' 

For such a design to exist, the. m us. be two blocks to contain 

two given varieties which are^ associates of both 
The other two elements in either block must oe 

The other two ^ such elements exist , and therefore the 

X and y (A 2 )• £ if we allowed this to happen, it 

two blocks would be identical. Even 11 

would be impossible in this case, because b is ocm, H 

the number of blocks if each appears twice. (This example is taken 

from [35].) 


Construction 

We can construct p.b.i.b.d’s by selecting blocks from suitable 
b i b d’s. In Chapter II we considered b.i.b.d’s whose blocks were 
hyperplanes in parallel pencils in a EG(m,s). This time we take only 
the m pencils * - c* 0-0, 1, ... , m-1), where the q are the ele¬ 
ments of GF(s). The points are varieties, and we define as ith asso- 
ciates any two points with m - 1 equal co-ordinates. 

We have v = s m , b = m.s, k = s m ~', r = m. The ith associates of 
the point , xj = x , say, are found by selecting m-i co-ordin¬ 
ates-this can be done in ways - and giving to each of the remain¬ 

ing co-ordinates any value different from that in x -which can be done 
in (s-l)* different ways. Therefore Two varieties 

which are ith associates and have, therefore, m - i equal co-ordin¬ 
ates, appear together in all blocks defined by any one of those co¬ 
ordinates. Hence h = m- i. 

Now to the parameters of the second kind. If two varieties, say 
x and y, are ith associates, then they have m-i co-ordinates with 
equal values, and i co-ordinates with unequal ones. A variety which 
is a yth associate of x and a /cth associate of y will have m - j co¬ 
ordinates in common with x and m - k in common with y. Of the 
m - j former ones, there will be u (any value from 0 to m - i) taken 
from those which are common to x and y, and m - j - u from the re¬ 
maining i co-ordinates of x, while m - k - u will be among the 
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m - (w - i) - (m - j - i) = i + j + u - m co-ordinates of y which do 
not have common values in x and y, and which have not been chosen from 
among those that are not common. The number of all these varieties is 
(m-i\ l i \(i+j + u-m\ 

[ u }\m-j-u)\ m-k-u ) 

Each of the remaining m — i — u co-ordinates out of those common to 
x and y can have any of s — 1 values, and each of the other remaining 
co-ordinates, numbering -u) - (m-k-u) - i + j + k-2m-2u = a, 

say, can take any of s-2 values (not those of either x or y). Thus 
we obtain 




i \ li-m + ] + u\\ 
m-]-ul \ m-k-u 


(s - 2)°. 


When s = 2, then only that term in this sum will differ from 0 in which 
a = 0, i.e. where u = l / 2 (i + j+ k-2m). 

When m = 2, then A, = 1, A 2 = 0, n, = 2(s—1), n 2 = (s-1) , 

( n _ (s-2 s — 1 \ ( _ ( 2 2(s-2)\ 

KP > k) \s-l (s-l)(s-2)J' ^2(s-2) (s-2) 2 j‘ 


This is a Latin square type scheme L 2 (see Chapter V). 

When m = 3, s = 2, then we can imagine the design to be derived 
from a cube with vertices 1,2,...,8. Let each block contain vertices 
of the same face. Two vertices are first associates if they lie on the 
same edge, second associates if they are diagonally opposite on a 
face, and third if they lie at the two ends of a diagonal of the cube. 
The design is then 1234 2367 5678 1256 1458 3478, where v = 8, 
J>=6, k=4, r=3, n 1 = n 2 = 3, n 3 =l, A 1 = 2, A 2 =l, A 3 =0. The parameters 
of the second kind are 



0 2 0 \ 

2 ° 1 , (Pi) 
oio/ 3k 




0 3 0\ 

3 0 0 . 

,0 0 0 / 


These designs are special cases of “quasi-factorial” designs with 
S\ xs 2 x,,.xs n varieties (>:, xj, x { = 1,..., s*, and such that 
(*i i ..., x n ) and (y,,. .., y n ) appear together in A c blocks with 

I 7T 

Q = 1 if xi = y it and c< = 0 otherwise. The above designs are ob¬ 
tained by suitably collapsing the association scheme. 

Other constructions of p.b.i.b.d’s from finite geometries will be 
given in the next chapter. 
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[ 123] contains examples of the construction of p.b.i.b.d s from 
difference sets, when v is the power of a prime. Bose and Nair have 
already considered such constructions in [16] for m = 2 and arbitrary 
values of v. The following example illustrates the procedure. 

Take 1,3, 9 (mod 13). The differences 2,5,6,7,8,11 appear once 
each, and 1, 3,4,9,10,12 do not appear at all. The blocks (1+*, 3 + *, 


9 + 0, t-0 


<■£> 


- (a 3) • <*> ■ 


12 form a p.b.i.b.d. with n 1 = 



Another example starts with 1, 2, 3 (mod 5). 


ti 2 — 6, A, — 1, A 2 — 0, 
We obtain a p.b.i.b.d. 


whose association scheme can be illustrated by a pentagram, where 
any two vertices connected by a line are second associates, and any 
other pair are first associates. 



In [9l] C.R. Rao mentions the following “circular lattice designs”: take 
n concentric circles and draw n diameters. Let the intersections rep¬ 
resent 2n 2 varieties, and let the circles and diameters represent 2 n 
blocks, with k = 2n varieties in each. Every variety will be twice re¬ 
peated. The association scheme is as follows: two points on diamet¬ 
rically opposed ends of the same diameter are first associates, two on 
the same circle, or on the same diameter, but not on both, are second 
..ssociates, and all other pairs are third associates. Then n, = 1, 
n 2 =4(n-l), rz 3 =»2(n-l) 2 ; A, = 2, A 2 = 1, A 3 = 0; 


0 0 


(p') =[04(n-l) 0 

3 ' 0 0 2(n-l) 2 


01 0 

<p£) = ( 1 2(n- 2) 2fo-D 

0 2(n-1) 2(n-l)(«-2) 


0 0 1 

(Pi) « I 0 4 4(n-2) 

3 V 1 4(n-2) 2(n-2) 2 


These designs are resolvable into two replications, one containing 
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the n circles as blocks, and the other the n diameters. 

In[93], P.V. Rao refers to [ill] where three associate classes are 
treated and where it is shown that certain "partially balanced” matrices 
can be used to find p.b.i.b.d’s with 2m + 1 associate classes. 

Various association schemes 

B. Harshbarger has developed designs (in[50] and [5l]) some of 
which are p.b.i.b.d’s. They were systematically studied by K.R. Nair 
in [72]. 

The “simple rectangular lattice” of v = n(n- 1) varieties may be 
described by allocating to each variety an ordered pair ( x,y ), x^y, 
x,y = 1, , n. We then let the first associates of ( x,y ) be the 

varieties which have either x or y in common with it (hence 
= 2(n-2)), second associates those which have no mark in common 
(hence n 2 = (rc-2)(n-3)), third those with one mark in common, but 
not in the same position (n 3 = 2(n-2)), and fourth the variety ( y,x ) 
(hence n 4 = 1). The blocks are (x,l). . .(*,«) for all x, and 
(l,y) •. • (n,y) for all y. Then b=2n, k = n- 1, r= 2, A, = 1, A 2 = A 3 
= A 4 = 0; 

(n -3 n- 3 1 0\ 

\ I n -3 (n -3)(n-4) n -3 0 | 

n- 3 n- 3 1 

0 10/ 

2(n - 4) 2 0\ 

(n-4)(n-5) 2(n-4) 1 \ 

2(n-4) 2 0 

1 0 0/ 

n -3 n -3 1 \ 

-3)(n-4) n -3 0 \ 

n-3 1 0 

0 00 / 

0 2{n - 2) 0\ 

(n-2)(n-3) 0 0 \ 

0 0 0 / 

0 0 0 / 

These designs are also mentioned in [ 70 ]. 

When n = 3, then no second associates exist. Renumbering the 
third and fourth association to be second and third, we obtain a scheme 


= 


c*S> - 


(* 4 ) - 


= 
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With three classes, and the parameters of the second kind are then 


0 1 o 

(pi ) = (10 1 

^ \0 1 o. 


oj!> - 



10 2 0 ' 
Cpp =200 
1 \0 0 0 , 


A “triple rectangular lattice” is constructed as follows: 

Takean[n(n-l),3,n, 2] : / (definedin Chapter III) and number 
its columns from 1 to n(n-V). The 3n blocks of the design are formed 
by all those numbers of columns which contain the same entry in the 
first, or in the second, or in the third row of the orthogonal array. Then 
k= n —1, and r= 3. 

Nair has shown that there are only two p.b.i.b.d’s among these 
lattices, namely those for n= 3 or n = 4. 


n = 3 


v= 6 , b= 9, k= 2, 
A,= 1, A 2 =0 

rci = 3, n z =2 

<*> - (" 0 ) - 


■=3 


3 0 

,0 1 


n = 4 




v=12, b= 12, fc=3, r=3 
A, = 1, A 2 =A 3 =0 

n 1 = : 6 , n 2 = 3, n 3 =2 

2 l\ /4 0 2\ /3 3 O' 
0 1. 0 2 0, 3 0 0 


n n 


n n 1 


The association schemes can easily be recovered from the para¬ 
meters. 

Other p.b.i.b.d’s with r = 3 are studied in [71]. 

Resolvable* p.b.i.b.d’s are introduced in [ 18] as follows. 

Consider the incidence matrix of a symmetric b.i.b.d. with v 
varieties and blocks, and k varieties in each block. Replace every 
cell with 1 by p different varieties, which are different in each cell, 
and every cell with 0 by q different varieties, again different in each 
cell. We have then v(kp + vq - kq) varieties, kp + vq - kq in each 
row of the matrix, and in each column. We take the rows, and the 
columns, as blocks and thus obtain 2 v blocks, and each variety is 
twice repeated. 

A rectangular association scheme is defined for v = s.t varieties 
as a rectangle with s rows and t columns, with first associates in the 
same row, second in the same column, and all other pairs being third 
associates. Hence n, = t-1, n 2 = s- 1, n 2 = (s- l)(t-1), and 


* Defined as for b.i.b.d’s. 
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- 



&£> = 


0 0 (-1 

0 S-2 0 

t -1 0 (s-2)(f-l) 



/ 0 1 f-2 \ 

I 1 0 s-2 J . 

\f-2 s-2 (s-2)(t-2)J 


Conversely, these parameters define the association scheme 
uniquely. Indeed, let x, and x 2 be first associates. Because n, = f-1 
and p^ = t-2 =71,-1, the s.t values can be grouped into s rows each 
of t elements, so that those in the same row, and only those, are first 
associates. This property is symmetric, and in the present case also 
transitive. 

x, has s — 1 second associates, and because p|=0 there will be 
precisely one second associate in each row of which x, is not a mem¬ 
ber. We can therefore write pairs of second associates into the same 
column, since p| = s-2 = n 2 - 1. 

A design constructed by taking as blocks the union of a row and a 
column of the rectangular association scheme has parameters v=b=s.t, 
r=k=s + t- 1, A, = r, A 2 = s, A 3 = 2. (This is a quasi-factorial design 
with n = 2.) 

A different design results when the variety common to the row and 
the column united in a block is omitted. 

In certain cases we can make an association scheme collapse into 
one with less classes. For instance, if in a rectangular association 
scheme s = t, and hence also A, = A 2 , and if we then let first associates 
be any two varieties either in the same row or in the same column, 
while any other pair are second associates, we obtain the L 2 association 
scheme (see next chapter). 

On the other hand, if we let two varieties in the same column be 
first, and any other two second associates, then we obtain a so-called 
“group-divisible” association scheme, which we study in detail in the 
next chapter. 

We obtain yet another association scheme if we denote 
v ~ fy x..,xr m varieties by Vv , , i , where ij = 1, 2, Tj . Let first 

associates be two varieties with all but the last subscripts equal, and 
generally ( m - i')th associates those with the first i subscripts equal, 
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but the (i+ l)th different. Then n f = r m - 1, n 2 -r m (r m - 1 1 ),..., 

^ r m-i • • • r 2( r i “I)* , . r i 

This scheme was introduced in [100] and discussed m [102]. In 

[88] it is proved that the values of the parameters of the second kind 

define the association scheme uniquely. These designs are special 

cases of the intra- and inter-group b.i.b.d’s defined m 1129]. We 

do not discuss this scheme any further here. In the next chapter we 

shall deal with a special case, namely that of group-divisible designs, 


Association matrices 

In [15] Bose and Mesner use another method to describe an asso¬ 
ciation scheme. They define it by Association matrices Bi, of order 
v and symmetric: 



where b a ti = 1 if the varieties s and t are ith associates, and = 0 
otherwise. The total of the entries of each column, and also of each 
row of B it is ri{. Clearly, 

(4.1) B 0 = I y (i.e. the identity matrix of order v). Furthermore, b 8 ti 
is unity for one and only one value of i when s and t are 
given, so that 

771 

(4.2) .2 Bi =J V (i.e. the matrix of order v whose elements are all 

t = 0 

1). A linear combination of the Bi can only be 0 if all the 
coefficients are 0. 

It follows from the definition of the b £ that 

m 

(4.3) pj k =^ o b° t Pj k if c and r are ith associates, because only 

one term in the sum of the right-hand side is different from 0. 
Consider the matrix product Bj B k . It contains in its rth 

V 

row and cth column the sum 2 bb° . A term in this sum is 

t =i V & 

1 if and only if t and r are y'th, and t and c are &th asso¬ 
ciates. Otherwise it isO. Therefore, if r and c are ith asso¬ 
ciates, then 


< 4 - 4) = 


Comparing this result with (4.3), we have 
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m m 

(4.5) Bj l B k .i ) B t tt jll - ( 2 B t p‘ kJ = 

This formula has appeared in [127]. We have not assumed 
that j and k must be different, and therefore any power of B* 
can also be expressed as a linear combination of the B*. We 
have 


= Bil $B„ = 2 ,B„ = 2 2 ^fi?„B t , 


and also 

«Wfl* = - 2 IpJp^B,. 

Matrix multiplication is associative, and the B i are linearly 
independent, as we have pointed out. It follows that 

m m 

(4 ' 6) J 0 P*pL = J 0 PijPL for 0 4 i,i,k,t 4 m. 

For t= 0, in particular, this means that n k p£ = n iP v jk . Given a 
p.b.i.b.d. with incidence matrix A, we have 

m 

(4.7) = 2 \Bi. 

t=0 

Conversely, let symmetric matrices of order v with ele¬ 
ments 0 and 1 be given, satisfying (4.1), (4.2) and (4.5) with 
some constants p? fc . Then these constants are parameters of an 
association scheme. 

We have to show that 

4 - p«; 

fb) Pq/c — 1 if t—k and = 0 otherwise; 

(c) t>? 4 = 0 for k*j, and pj> = Uj , the total of every row 
of B, ; 

(d) ? £ - 

W = "iPi- 

Proofs: 


(a) B k B, - B k B] = = ilp^Bi)' _ 2pi Bj - 

- Bj B k , and hence (i. 


(b) B 4 _ B„Bj, « 2 p' t B t , hence (b). 

V t 

(e) .1 = 0 for j^k because of (4.2), and equals the 

row total of the tth row of Bj when , = fc. Because of 
U5), Pjj IS the total of every row of Bj , say n jt while 

= 0 for j*/c. From (4.2) it follows that 2n =v. 
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(d) 

(e) 


(? Bj)B k = J,B k = and aIso 

i(BjB k ) = 11 P^B t = l(lpOB t . 

Comparison of coefficients establishes (d). 

Consider the matrix products MjM k and M k Mj, where 


U i " 


Pq; 


mj 



Because of (a), the two matrix products are equal, and from 
this fact we deduce (4.6) and, in particular, (e). 

■% 

Self-dual designs 

Let a symmetric p.b.i.b.d. be given with incidence matrix A and 
with association matrices = 0,..The dual design will have 
the same parameters v=b and r=k. In [54] A.J. Hoffman has derived 
sufficient (though not necessary) conditions for the dual to have also 
the parameters p* k and A^. 

V/e establish conditions which ensure that symmetric matrices, say 
Q, with elements 0 and 1 exist, which satisfy (4.1), (4.2), (4.5) with 
the same parameters as the B., and satisfy also, instead of (4.7), the 
relation (4.7a) A'A = A^Q. 

We prove that the matrices Q = A~ 1 A satisfy these conditions, 

provided the following conditions hold: 


Condition I A is not singular (i.e. A~ ] exists). 

Condition II The greatest common divisor of all determinants of 
order m contained m the (m + 2) by m matrix 


is unity. 

Proof of symmetry 
AA'B t 

Hence AA! B t 

C t 



= IA iBiB t , and B t AA' = 

I i 

= B t AA' , and B t A = AA'B t A'~\ Therefore, 
= A-'B t A = A-'AA'B t A‘-' = C t ‘. 
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Proof of (4.1) C 0 = A~ y B 0 A = A~'l y A = l y . 

Proof of (4.2) IQ = 2 = A~'J V A. 

We have to show that this equals J v , and this is a consequence of 
J y A = AJ y , which will now be proved. 

We have AJ V = rJ v . Also A 1 J v = A~'AA'j v = = cJ v , 

where c is a constant, equal to the sum of the entries in any one 
column of A. Since the total of all entries of A is rv, it follows that 

c=r. But A J y = rJ v as well as J y are symmetric. Hence J y A = rJ = 

AJ V - 

Proof of (4.5) 


c i c * = (A-'BjJtXA-'BkA) = XA-'p‘ k B t A = }p‘ k C t . 

Proof of (4.7 a) A 1 A = A~' AA‘A = A"' 1 kiB t A =2X i C j . 

t % 

The Q will be association matrices if their elements are 0 or 1. 
Multiply (4.7) on the left by A’ 1 and on the right by B f A. Then 

A’B jA = A'lX^A = SSA.piC, (, = 0,1.»). 


Imagine these matrix equations written down explicitly for all ele¬ 
ments of the Q. Reference to a theorem in [ 122] proves that these 
elements are all integers, by virtue of Condition II. We show now that 
these elements are 0 or 1. The sum of the squares of all elements in 
C i is equal to the trace (i.e. the sum of the diagonal elements) of 

c i c i = C { . By virtue of (4.5) this is also the trace of Id { C Th* 

£ ^ £ * 1 AlC 

trace of C t equals that of B t , i.e. 0 for r*0, and v for ( = 0. Hence 
the sum of the squares of the elements of C, is rqv. 

The sum of all elements of C, is also n 4 v, because all row (and 
column) totals are n t . But the sum of integers can only be equal to 
the sum of their squares if all elements are either 0 or 1. 

For m= 1 this theorem means that the dual of a symmetric b.i.b.d. 

is again a b.i.b.d. with the same parameters. This is a result which 

je haye already obtained. (In this case the conditions are always 
satisfied, provided r>\.) y 

ab ,e X t h o e ; K eth : d r f bl0Ck Secti0n and block intersection are applic- 

C.R Nafr'has d S “ 7“’ Pt ° Vi<led they “ e ° f linked block type, 
new'des ““ cases whete ‘he association scheme of the 

esign equals that of the original one (see [69]). 



66 


The dual of a b.i.b.d. with incidence matrix A is the.b.i.h.d- 
Whose incidence matrix is A'. The latter design is also said to have 

been obtained from the former by inversion. 

If the parameters of the original b.i.b.d. are denoted by v,b k, 
and r, with their usual meaning, then the dual design has b varieties 
V blocks of r varieties each, and each variety is repeated k times. 
Any pair of varieties will appear in the new design as often as there 
were common varieties in the blocks to which the two varieties corres- 

pond by their duality. . 

We shall now quote cases in which configurations turn, by inver¬ 
sion, into p.b.i.b.d’s. 


Case(l) ([101]) / (v-2) t 

A b.i.b.d. with v varieties, b=Cj c ), r= (J_i) » ^ = U-2/ ^ an 
“unreduced”b.i.b.d.), turns into a p.b.i.b.d. unless v = k+ 1, or *-l. 

Let the varieties be 1,2, ...,v, and take the block B 
There are blocks with no element in common with B. There 

afe blocks with only the variety 1 in common with B, and 

k ^ blocks with precisely one element in common with B. 

Similarly, there are blocks with precisely i elements in 

common with B. Of course, 2, (jt_f) = (k) = as ^ mus * 
be. (We have included B itself in this count.) It will be noticed that 


( = 0 when 1 > k. or when i <2k-v. 

Upon inversion there will be (J) {k-f) varieties which a PP ear 
i times in the same block with any other variety. If two varieties 
which appear k-i times together are ith associates, then A i = k-i 

wise we again have a b.i.b.d.) These numbers are independent of 
which variety we started from, and so are the p* which we are now 
going to determine. (We use t as the second subscript, because k is 
used for the number of varieties within a block of the original b.i.b.d.) 

Consider, in the original b.i.b.d., two blocks with k-i varieties 
in common. Let these be a 1f ..., a k ~ it and let the remaining varieties 
in the first block be b 1f ..., and in the second block cy ,..., c». 
Denote those varieties which are not within the sets just mentioned 


-k 
i 


j . (It is here that v-k > 1 is necessary; other- 
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by d ,,. ,dy-k-i . We want to find a block of k varieties with t from 

among the a’s or the b* s, and / from among the a's or the c’s. 

If s varieties are taken from among the a’s, k — t — s from among 
the b’s, and k-j-s from among the c’s, then there remain t+j + s — k 
to be taken from among the d’s. Given s, this can be done in x 

(k-t-s) (fc—;'-s) (f+/> s ~-/c) different ways, and since s may 
be 0, or 1,..., or k (if we take too high a value, then the corres¬ 
ponding binomial coefficient will be zero), we have to sum the above 
product of four binomial coefficients over all those values of s to 
obtain p x jt . 

Case (2) ([lOl] and [113]) 

If A=l, then the inversion of a b.i.b.d. produces a p.b.i.b.d., 
unless the original b.i.b.d. is symmetric, in which case its dual is also 
a b.i.b.d. 

Take a block B and one variety in it, say x. There are r-1 more 
blocks in which x appears, and this is true whichever of the varieties 
in B we take. However, since X=l, any other variety will reappear in 
r— 1 blocks different from those in which x reappears. Hence, if after 
inversion we call first associates those varieties which appear in the 
same block, and all other pairs second associates, then n, = fc(r-l), 
n 2 = b-k(r- 1)-1, A 1= l, A 2 =0. 

For a symmetric b.i.b.d. we have fc(r-l) = r(/c-l) = v-1, so 
that n\ = v-1, i.e. all varieties are first associates. 

To find the parameters of the second kind, we determine first p*. 
From the original blocks let us take two with no variety in common, say 
the blocks a = (a\, ..., a k ) and b = (b\ ,..., b^). There is precisely 
one block which contains a* as well as bj, and this block has precisely 
one element in common with block a, and also with block b. The sub¬ 
scripts i and j take, independently, k different values. Therefore 
p, t = k , and from this we obtain, through the relations which connect 
the parameters, 



r-2+(/c-l) 2 (r-k)(k-l) \ 

(r-/c)(fc-l) b + r-k(2r-Jc)-l) 



( * 2 

\ k(r-k- 1) 


k(r-k-V) \ 
b+ k 2 -2k(r- 1) -2/ 
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, ^ use of the relations b = vr/k and 

In this derivation we have made use 

r(k- 1) + 1 - v ‘ . special case of (1). In [lOl] P.M.R oy 

,£ fc=2, then we h n Tib“s by inverting b.i.b.d’s. In[U3]S.S. * 

has constructed more p.b. . ■ y are dua , s of b.i.b.d’s with \ = 1 
Shrikhande has studied p.b.i.b.d swlucn 
obtained from systems of difference sets. 


Case (3) ([lOl] and [113]). . 

The inversion of a b.i.b.d. with A=2 which was obtained from 

symmetric b.i.b.d. by block section leads to a p.b.i.b.d. with para- 
. fk\ - ou A. = 2. A, = 1, and 


a 


meters ri\ = 



For the proof, we refer the reader to the original sources. 

Case (4) ([113]) / a / r \ 

The dual of a b.i.b.d. with v = [ 2 J , b = ( 2 j , = r " 2 ’ A = 2 

is a p.b.i.b.d. with n 1 = 2(r—2), ri 2 = ( 2 ) ’ ^ — 


/ r-2 r-3 \ / 4 2(r-4) \ 

(p i ) = (r-3 ( r 2 3 ))’ (P ^ = (2(r-4) ( r 2 4 ) ) 

The proof depends on a theorem in [58]. We refer the reader to the 
source given above. S.S. Shrikhande has shown in [119] that the 
p.b.i.b.d. thus obtained is always of triangular type (see Chapter V). 


Case (5) 

Consider an affine resolvable b.i.b.d. with v=n 2 m, b = nr, k=nm. 
After inversion we have a p.b.i.b.d. with n, = b-n, n 2 = n-1, A, = m, 
A 2 = 0, 



(b-2n rc-l\ 

U-l 0 / ' 



Case (6) ([ 116]) 

If from the affine resolvable configuration of case (5) we omit all 
blocks that contain a particular variety, then we obtain a symmetric 
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configuration with n 2 m - 1 = (n-l)r varieties and blocks. The blocks 
contain nm = r —A varieties each, and each variety is repeated ntn 
times. The dual of this design is a p.b.i.b.d. with A, = 0, X 2 = m, 
n \ = n-2, n 2 = n m-n. A further examination of duality and partial 
balance is contained in [94]. 

Examples for (5) and (6) 

The design 


123 145 167 

18 9 

468 269 258 

247 

579 378 349 

3 5 6 

is affine resolvable. Inverting, we have 


1 4 7 10 2 4 9 11 

3 6 7 11 

1 5 8 11 3 4 8 12 

2 6 8 10 

1 6 9 12 2 5 7 12 

3 5 9 11 

", = 12-3 =9, n 2 = 3-1 = 2, A, = 1, 

O 

II 

Ol 

°&> " (2 o)- 04) 

(9 0' 

\0 1 


Omit all blocks with 9, say; we obtain 


123 167 145 247 

468 258 378 356 

Inverting, we have, with simplified notation, 

1 3 5 2 3 7 4 5 7 1 6 7 3 6 8 2 4 6 1 4 8 2 5 8 



It should be noted that in all the cases (1) to (6) the parameters 
y,b,r, and k refer to the original, and not to the inverted design. 


EXERCISES 

(1) Construct a design from the simple rectangular lattice for n = 4. ? f 
etermine the parameters of the p.b.i.b.d. derived from the 
difference set 1,2,3 (mod 5)(see page 58). 
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(3) Find the association matrix of the scheme of the following design, 
and determine A A' and all BiBj. 



Invert 


3478, 1234, 2367, 5678, 1256, 1458. 


A 0 = 

= 3, 

^ = 2, 

A 2 =l, A 3 -0, 

riy =n 2 

— 3, Wg 

= 1. 


12 

23 

123 

145 

167 

189 

(i) 

13 

24 

(ii) 468 

269 

258 

247 


14 

34 

579 

378 

349 

356 


(5) Construct a p.b.i.b.d. by inverting the design obtained by block 
section from 


12345 1569E 34689 

12678 236XE 3578 E 

1379X 2479E 4567X 

148XE 2589X 

v=b= 11, r=k = 5, A=2. 

(Example from [54].) 

(6) Find the two triple rectangular lattices which are p.b.i.b.d’s (see 
page 60) and determine their association schemes. How would 
you describe them as configurations? Invert the first p.b.i.b.d. 
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CHAPTER V 

PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS 
WITH TWO ASSOCIATE CLASSES 

Introduction 

To begin with, we repeat for convenience the formulae for the 
parameters of a p.b.i.b.d. as they apply to designs with two associate 
classes: 

vb = rk, «! + n 2 = v - 1, n, +A 2 n 2 = r(k— 1), 

P 11 + P t2 = "l ~ 1 > Pf\ + P 12 = n V P\2 + P& = ^ P,| + P 22 = n 2~ 1 

n ' p 22 = n 2Pf 2 'hPia = n 2 P 2 

When the parameters and n 2 are known, then only one of the 
Pj k remains independent. 

The definition of the demands that they should be independent 
of which pair of ith associates we start from. In [12] Bose and 
Clatworthy have pointed out that for designs with two associates it is 
sufficient to demand this independence only for p^ and p 2 ; that of 
the others will then follow. If x and y are any two varieties which 
are first associates, then let p k (x,y) be the number of jth associates 
of x and /cth associates of y. We have then 

p\ (x,y) + pl 2 (x,y) = n, - 1, 

P\ y (x,y) + p^(x,y) = n, - 1, 

and P' 2 ,(x,y) + p 22 (x,y) = n 2 . 

It follows that if p^ ( x,y ) is independent of x and y, then this is 

also true of all the other p* , and the same follows for p*(x,y) in a 
similar way. 3k 

On the other hand, it is possible that pf k (x,y ) is independent of 

the pair x and y, but p\(x,y) is not. The following example is given 
in [ 12 ]: 


01 12 23 34 45 56 60 

03 14 25 36 40 51 62 
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This proves the statement, if first associates are varieties which 
appear in the same block. In this example pJ/O.l) = 1, for instance, bllt 

P ' ,( °Adesign with equal block sizes (k) and an equal number of teplj. 
cations (r) of all varieties is called a (6*, v r ) configuration. If if is 
symmetric, i.e. b=v, and r=k, then we denote it by (v,) (see [129]). 

Let varieties in the same block be first, and other pairs second 
associates. Then A, = 1 and A 2 = 0. 

Examples A (7 3 ) is the PG( 2,2), a b.i.b.d. The (8 3 ) 

124 167 235 457 
138 278 346 568 

is a p.b.i.b.d. with parameters v=b = 8, r=/c = 3, n y = 6, n 2 =l and 

= (l o)' ( p P = (o o)- 

There exist three different (9 3 ), with v=b-9, r-k-i, n^-b, 
ri 2 = 2 and the following values of 


4 1 

and (pf k ) 

for the configuration 

5 1 

123 145 189 248 269 

1 1 

1 0 

356 379 467 578 

3 2 

5 1 

123 145 168 257 289 

2 0 

1 0 

356 349 478 679 

3 2 

6 0 

123 148 159 247 269 

2 0 

0 1 

357 368 456 789 


(This is the configuration illustrating 
the theorem of Pappus) 


The Desargues configuration illustrates a (10 3 ) which is a 
p. b.i.b.d. with v=b= 10, r=k = 3, n^ = 6, n 2 =3, (p} k ) = (| l) > 

( p P = (2 0) ' 

Other p.b.i.b.d’s can be obtained from the points and faces of 


polyhedra. 

A design with either A, = 0 or A 2 = 0 is called simple in [ 13]. A 
thorough examination of simple designs with k > r > 2 is contained in 
[12]. See also [8] for k> r= 2. An enumeration of p.b.i.b.d's with 
two associate classes and r=3 is given in [99], 

If in a projective geometry we take all flats of a given dimension 
as blocks, and all the points on it as varieties, then we obtain a 
b.i.b.d., as we have seen. If we omit suitably chosen blocks or 
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varieties, and call first associates those varieties which appear to¬ 
gether in one of the remaining blocks, while any other pair of points 
are second associates, then we obtain simple p.b.i.b.d’s. We shall 
give a few examples of this procedure. 

Omit one point, say P, of a PG(m,s) and all f-flats through it, 
but no further points in the latter. There remain v = s(s m - l)/(s-1) 
varieties. 

Originally we had V(m,t;s) different t-flats. Of these we have 
omitted all those through P, i.e. V(m- 1, t-lys), so that there 
remain 

b = F(w-1, t-1; s)(s m+ i-5 e+1 )/(s t+1 - 1) 

(s m - 1) ... (s m “ t+ i_ 1) 

= s t+1 —T-• -. 

(s f -l) ... (s-1) (s t+1 - 1) 

Every one of these blocks contains k = (s fc + 1 - l)/(s-1) varieties. 

It follows that the number of times each variety appears is r = bk/v = 

( s m - st) ( s m ~ 1 - 1) ... ( s TO " t+1 -1) 

- 1) ' (s t_1 -l) ... (s-1) * 

Two varieties will appear in the same t-flat unless their line 
passes through P: then and then only will all the *-flats have been 
omitted in which these appeared together in the original PG{m,s). 

Thus n 2 = s-1 (all points on the line through P and through a chosen 
point, excluding these two) and hence n, = s 2 (s m_1 -l)/(s-1). j s 

the number of f-flats through a straight line, less those through this 
line and P, i.e. 

A, = V(m-2,t-2,s) - 7(m-3,f-3;s) 

" ”(■$*-1-1) * (s fc -2-1) ...(s-1) * 

i We turn now to the parameters of the second kind. We have 
^ 22 ~ 0 , because a variety’s second associates are all in the line 
trough its corresponding point and P, and only two points which 
a re on a line with P can have a common second associate. From this 
e °ther values follow: 
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• • . + 
5-1 




s m + ... + S z 
0 



We notice that these latter parameters are independent of t. If, in 
particular, t = m- 1, then the expressions reduce to v= s(s m -l)/(s-l), 
b = s m , k = (s m -l)/(s-1), r = s™' 1 , = s 2 (s m ~ 1 -l)/(s-l), n 2 = s- 1, 

A, = s m-2 A 2 =0. If m = 2, f=l, then v=s 2 + s, ft=5 2 , &=s+l, r=s, 

«i= s* n 2 =5-l, A, = l, Ag = 0. 

More generally, we could start from any b.i.b.d. with parameters 
v,b,r,k, A= 1 and omit one variety and all blocks that contain it. We 
obtain then a p.b.i.b.d. with v- 1 varieties, b-r blocks, r-1 repli¬ 
cations and A, = 1, A 2 =0. 

An analogous construction can start from a EG(m,s ). 

We omit one point of the EG(m,s ) and all the f-flats through it. 
There remain v=s m -l points. Originally we had V{m-l,t-l;s)s m ~ t 
different t-flats. Of these we have omitted all those through a chosen 
point, V(m- 1, t- 1 ;s) in number, so that there remain b= V(m- 1, t- l;s)x 
f-flats. All of them contain k = s t points, hence r= bk/v = 
V(m-2,t-l;s)s t . By the same argument as before, we now obtain 
n 2 = s- 2, and hence n^ = s m ~s. Also we have A, = r( A - l)//j, = 
V(m-3,t-2;s)s t ~\ unless m = 2, when A, = l. Again, P 22 =0, so that 






(When s = 2, then only first associates exist, and we have a b.i.b.d.) 

In particular, when t = m- 1, then v = b = s m - 1, r=k= s m ~\ A = s m ~ 2 . 

A p.b.i.b.d. where p’ z , p 2 and p 2 2 are not all zero can be derived 
from the correspondence between a point (y Q ,... ,y m ) and the hyper¬ 
plane 2(2 a tj yj)Xj, =0 with skew-symmetric matrix (a^ ). Let m be 

odd, and let the points of a PG(m,s ) be varieties, and let blocks con¬ 
tain the varieties of those points on a line in a hyperplane passing 
through the point corresponding to the latter. There are then 
v = (s m+1 -l)/(s-l) varieties. On every line there are k = s+1 points, 
and through each point there are r = (s m_, -l)/(s-1) lines of the hyper¬ 
plane; this is the number of replications of the variety corresponding to 
that point. Hence the number of blocks is 
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b = vr/k = [(^ +1 -l)( S --i-l)]/[( s _l)2( s+ i)]. 


Let first associates be varieties in the same block, and let any 
other pair be second associates. We then have a p.b.i.b.d. with 

rii—s(s — l)/(s—1) (i.e. s times the number of lines through a point 
in the hyperplane), n 2 = (s^-l)/( S -l) - 1 = s™, 

A, = l, and A 2 = 0. 

The parameters of the second kind are determined as follows: two 
points which are first associates lie each on the hyperplane correspond¬ 
ing to one of them, and on that corresponding to the other. Hence those 
other varieties which are first associates of both of the former ones lie 
in the intersection of the two hyperplanes. Excluding the two points 
themselves, we have thus p, 1 , = (s^-l)/(s-l) - 2, and hence 


ts m ~ 1 -2s+l) 





s m- i(s-1) 







This construction is a fairly obvious extension of that given by 
Clatworthy, in [34] for m = 3. For m = 3, s = 2 the p.b.i.b.d. is men¬ 
tioned in [19], It is of triangular type, see later, page 84. 

A further generalisation, taking flats of higher dimension instead 
of the lines used here for blocks, is contained in [l]. The values of 
Pjh and Pf k ar e again independent of the dimension of the flats. 

We mention one more construction of a p.b.i.b.d. from a PG(m s) 
Choose a point P, and t of the (s»- l)/(s- 1) lines through it. Let 
the other s points on these lines correspond to varieties, so that 

there are v=s.t varieties altogether. Let the b s™" 1 _ m 

s -1 “ 5 

hyperplanes not passing through P correspond to blocks. Then each 
point that corresponds to a variety will be contained in r = - sm ~ 1 
s m "i-l S_1 

S — 1 =s hyperplanes representing blocks, and each block will 

contain k=t varieties. We call “first associates” points on one of the 
t chosen lines, and any other pair “second associates”, so that 
*i = s-l and n 2 = s(t-l). Also A 1 = 0 and A 2 =s m ~ 2 . We find 




9 



0 s-l\ 
s-1 s(t-2)j ' 
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Another method of obtaining families of simple p.b.i.b.d’s was 
described in [109] by E. Seiden. 

We have seen that in a PG(2,2 n ) there exist sets of 2 n +2 points 
each (a conic and its nucleus), no three of which are collinear. There 

are /2 n + 2\ = (2 n-1 + l)(2 n +l) lines through two points of such a 

\ 2 r . 2 ^ ^ . 

set, and no lines with precisely one, hence there are (2 + 2 +1) - 

(2 n-1 + l)(2 n +l) = 2 2n_, -2 n “ l lines with no point of the set on any of 
them. 

We take as varieties the 2 2n —1 points outside the chosen set and 
as blocks the lines containing two points of the set (each point lies on 
2 n ’+1 of these lines). This gives us one family of p.b.i.b.d’s. In a 
second design we take as blocks the lines containing no point of the 
set (each point lies on 2 n_1 of these). 

Let two varieties be first associates if they appear in the same 
block, and second associates if they do not. Then Aj=l and A 2 =0. 
The other parameters are found to be 


1st design 


2nd design 

2 Zn -1 

V 

2 2n -l 

t_1 + l)(2 n +1) 

b 

2 n ~\2 n -l) 

2 n_1 +1 

r 

2 n_1 

2 n -l 

k 

2 n +l 

2 2n-1 - 2 

n y 

2 2 n-1 


2 2n-z -3 

2 2n-2 


2 2n " Z 
r\2n~2 


22n-2 _ j 2 2n ' 2 _ 1 
22n~2 _ i 2 2n ~ 2 


n 2 

oi> 
«) 


2 2n-1 _2 

22n“2 2 2n “ 2 - 1 

»2n-2 2 ^ “2 j 

22n~2 2 ^ti“2 

22n-2 ^ 


In [95] p.b.i.b.d’s are constructed from the geometry of quadrics 
in spaces of higher dimension. 

Shrikhande and Singh have pointed out in [ 121] that if in a p.b.i.b.d. 
with two associate classes p' n = p 2 , then a symmetric b.i.b.d. can be 
constructed by forming v blocks from the first associates of each 
variety. We obtain a b.i.b.d. of linked block type with /c = n, = r. 

Example 

We start from a b.i.b.d. with A=l, r=2/c+l, and hence v=/c(2/c-l), 
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b - 4k 2 - 1, and take its dual (see case (2) on page 67, obtaining a 
p.b.i.b.d. with v = 4k -1, n 1 = 2fc 2 , p ^ k \ p 2 _ fc a if we then form a 
symmetric b.i.b.d. by the method just described, we obtain a Hadamard 
configuration (see page 14) with v=4k 2 -l, r=2k z , X=k 2 . 

In [10] R.C. Bose has applied graph-theoretic methods to the study 
of association schemes of p.b.i.b.d’s with two associate classes. 


CLASSIFICATION 

In [19] the authors have classified p.b.i.b.d’s with two associate 
classes according to their association schemes. We shall now consider 
this topic. 


(A) Group-divisible designs 

We start with the group-divisible (GD) scheme and shall call a 
p.b.i.b.d. with such a scheme a GD-design. 

Let v = m.n and write the varieties into m columns (called groups) 
and n rows. Let two variaiK&"in the same column be first, and any 
other pair second associates. Then n 1 = n-l, n 2 = n(m-l). Also 
p’ 2 =0, because the relation of being first associates is transitive (as 
in the rectangular scheme, of which this is a modification, already 
mentioned in Chapter IV). From this, or just as easily from first 
principles, we obtain 



° ) 
n(m-l)J 1 



"- 1 V 

n(m-2) ) 


We have already met this design, derived in a different way, on page 
75. Conversely, if a p.b.i.b.d. with two associate classes has p 1 =0, 
then it must be a GD-design. We have then = n y - 1, so that a^first 
associate of any variety is also a first associate of all its first asso¬ 
ciates. Hence the v varieties can be divided into groups, each of 

71 1 + 1 = ”, sa y, varieties, and those in the same group are first 
associates. 

If in a design p 2 =0, we again have a GD-design, but have to call 
all varieties in the same group second associates. 


A fairly obvious way of constructing GD-designs is to start from a 
b-i.b.d. with parameters v,b,r,k,\ and to replace each variety by n 
new ones. We have then altogether v.n new varieties, and if any two 
these which have replaced the same original variety are called first, 
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and any other pair second associates, then we have a OD-design with 

A.=r, A 2 = A, and m=v (see Exercise V.5). 

If in a GD-design A 2 =0, then it consists of disconnected groups 

of blocks, i.e. of m isomorphic b.i.b.d's each of n varletles ' 

Families of symmetric GD-designs are given in [35] as follows: 

(1) h=v=2n blocks (a,,.. .,n„,*i) and ( b w- ’ b < 
a = 1n) with the association scheme 


'n» 


a, ft, 

0£ ftg 


a n ftn 

Then m = 2, r=fc=n + l, n, = n-l, ^2 = ”* = ^ = 2 ' 

(2) b=v^3n blocks (a,,..., a^, ft*), (ft,,..., ft n » c t)» ( c v • • •» c n» a i) 
(i‘ = 1,..., n) with the association scheme 

o, ft, c, 


ftn ^n 

Then m = 3, r = k = n+ 1, n, = n— 1, n 2 = 2n, A, =n, A 2 = 1* 

Symmetric GD-designs can also be constructed from a set 

of integers mod v=m.n such that among the k(k- 1) dif¬ 
ferences every integer congruent to multiples of m appears A, times 
and all others A z times. (As regards a necessary condition for the 
existence of such sets when n = 3 (mod 4) see [ 129], page 22.) The 
ith column of the association scheme is then (i, m + i, 2m + i ,..., 

(n -1) m +1), and the blocks of the design are (d^ + t, ..., d k + t) 
t=0__ v— 1. 

In [124] D.A. Sprott constructs a GD-design as follows. Consider 
a b.i.b.d. with v = ft = s 4 + s 2 +l, r=k = s 2 + 1, A = l, represented, for 
instance, by the lines of a PG(2,s 2 )- We omit the s z +s+ 1 blocks 
given by the lines of a PG(2,s) contained in the PG(2;s 2 ),and the 
s 2 +5 + 1 varieties corresponding to the points of the PG(2,s). We then 
obtain a GD-design with v = ft = s*-s, k = r=s 2 , m = s 2 +s+ 1, and 
n=s 2 -s. 
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Example 

(s = 2). Association scheme 



1 

3 

5 

7 9 11 

13 



2 

4 

6 

8 10 12 

14 


12, v = b 

•% 

rH 

II 

r=k 

= 4, 

A 1= 0, A 2 

= 1. 


2 

3 

5 

9 

1 

4 

6 

10 

4 

5 

7 

11 

3 

6 

8 

12 

6 

7 

9 

13 

5 

8 

10 

14 

1 

8 

9 

11 

2 

7 

10 

12 

3 

10 

11 

13 

4 

9 

12 

14 

1 

5 

12 

13 

2 

6 

11 

14 

1 

3 

7 

14 

2 

4 

8 

13 


In [l9] we find this example, and another non-isomorphic GD-design 


with the same association scheme, namely v 

= 14, b= 8, 

V 

II 

0 

■r 

11 

to 

• 







1 

3 

5 

7 

9 

11 

13 

2 

4 

6 

8 

9 

11 

13 

1 

3 

6 

8 

10 

12 

13 

2 

4 

5 

7 

10 

12 

13 

1 

4 

5 

8 

9 

12 

14 

2 

3 

6 

7 

9 

12 

14 

1 

4 

6 

7 

10 

11 

14 

2 

3 

5 

8 

10 

11 

14 


We have seen in Chapter IV that the dual of a b.i.b.d. with \=1 
is always a p.b.i.b.d. (see case (2) on page 67). If we start, in 
particular, with s 2 varieties, s 2 + s blocks, s varieties in each block, 
and s+l replications of each variety (lines in a EG(2, s)),then inver¬ 
sion produces a p.b.i.b.d. with v=s 2 +s, b=s 2 , r=s, k=s+l, A,= l, 
A 2 = 0, n, = s 2 , n 2 = s- 1, 






This is a GD-design with m = s+1, n=s where the role of the first 
a nd the second associates have been reversed. It is, incidentally, of 
linked block type. 
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The design obtained from an affine resolvable b.i.b.d. by inversion 
(see Chapter IV, case (5), page 68) is a GD-design if we exchange 
the names of first and second associates. That obtained in case (6) is 
a symmetric GD-design. 

If in an orthogonal array [An 2 , m,n, 2] we replace an element x in 
the ith row by x+(i-l)n, and take the columns as blocks, and the 
elements of the same row as belonging to the same group of the asso¬ 
ciation scheme, then we have a GD-design with v=mn, b = A 2 n 2 , 
r=A 2 n, k-m, Ai = 0, A 2 = A. 

If a completely orthogonal set of Latin cubes of the first order, of 
side s, exists, then a GD-design with v=s 3 , m=s, n=s 2 , b = s (s+1), 
r = s(s+ 1), k = ns 2 , A 1= s, A 2 =s+1, also exists. For the proof, write 
the varieties of the ith group (i = l,..., s) as a square of side s. We have 
then s 2 varieties in 5 squares. Superimpose the s 2 +s -2 cubes on 
this arrangement. Picking up the treatments coinciding with the several 
elements of the 1st, 2nd,...»cube and writing them in different blocks, 
we get s(s 2 + s-2) blocks, each of s 2 varieties. Add 2s further 
blocks, formed by writing into the ;th block the varieties in the ;th row 
of all the squares, and into the (s + /c)thblock the varieties in all the 
&th columns. The design obtained is of GD type, with the parameters 
as quoted above. 

Consider now the incidence matrix of a GD-design, denoted by A. 
Then 



and F = A 2 4 1 . Both E and F are n by n matrices. 

To evaluate the determinant |AA'|, add the 2nd,3rd,.... , mnth row 
to the first. The latter row has then every element equal to rk. We take 
this factor out and subtract the first row A 2 times from all other rows. 
Then we have 



1 1 ... 1 


CM 

1 

• 

• 

• 

1 

S< 

1 

V- 

II 

Ai-A 2 r-A 2 ... A 1 — A 2 

• 

1 

IT 

1 

• 

• 

• 

1 

V 

IS) 


Ai-- A 2 A t -A 2 .., r—A 2 


A,-A 2 A^ — Ag ••• r-A 2 


The column sums of the second matrix are all (r- A 2 ) + (n-1)(A, - A 2 ), 
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and taking this factor out, we are left with the first matrix. Thus: 

1 1 ... 1 
A,-A z r-A 2 ... A,-A. 


\AA'\ = rk 


- ^2 A-i - A 2 

k m(n-1) 


r-Ar 


[(r-A 2 ) +(n-l)(A 1 -A 2 )] m " 1 


This equals r«r-A,) (rk-vX^T', as can be seen by subtracting 
the first column from the others in the matrix, and by noting that 
r^ 2+ (n-l;(A, A 2 ) ,+ *,„,_(! * r + A ^ 

A 2 (l + n 1 +n 2 ) = rk~v\ 2 - 

Types of group-divisible designs 

^classified according to the sign of the factors of 
1 (see [14]) tit is, of course, always positive. r-A, cannot be 

n T ™' IT ” the " We Cal1 the desi 8" gpgulvj In fact, any 
p.b.i.b.d. with two associate classes and r = X, must be a GD-design 

because in this case all first associates of a variety appear in a blocL 
if and only if that variety appears in the same block. Hence the re- 
ation of being first associate is transitive, p^ - 0 (and p* = 0), and 
the design must be a GD-design (see [l4]). 

In a singular GD-design all varieties that are first associates of 
one another can be telescoped into a single variety: we then obtain a 
b.i.b.d. If the parameters of the original p.b.i.b.d. were v=m.n, b, 

k > and V then the resulting b.i.b.d. will have v' = m varieties 
r = r replications, b=b blocks of k‘ = k/n varieties each, and any ’ 
pair o varieties will appear A'= A 2 times in the same block. It follows 
trom b >v' that in a singular GD-design we have b > m. 

Furthermore, in such a design we must have rk~v\ 2 > 0, because 
the left-hand side equals ( r’k ' - v'A')n, and this cannot be negative. 

1 r> A,, then again rk - vA 2 ^.0. We can see this as follows. 
Consider the submatrix A v consisting of the first 2 n rows of A, 
l e ‘ of those referring to the varieties of the first two groups. Then 

= |p 3 = [r + (n-l)A 1 + nA 2 ](r-A 1 ) 2l "-’ l (rk-vA 2 ), 

But I'VV must be non-negative, because 4, has at least as many 
lumns as it has rows* (see footnote on following page). 

Since the expression in square brackets is positive, rk - v\ 2 ^ 0 
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follows. A GD-design with r - A, > 0, rk - vA 2 > 0 (and hence \AA'\ > o) 
is called ^regular ])and one with r - A, > 0 and r/: - vA^O^Cand hence 
|AA'| = 0 asTfor a singular design) is calleddem^re^wS^ 

For regular designs.it follows, as for b.i.b.d s, that b ^ v, since 
the rank of AA!, namely v, cannot be higher than that of A. If a regu¬ 
lar GD-design is resolvable, then we can again derive b>v + r- 1, as 
on page 13. 

Take those columns of A which correspond to the blocks of one 
complete replication. The entries of each row add up to 1, and there¬ 
fore not more than b~r +1 columns can be independent. Hence 
b-r+ l^v, or b^v+r-1. 

Structure 

In [37] Connor chooses t < b blocks, and calls their incidence 
matrix A 0 . He then evaluates the determinant 

= (rk - vA 2 ) m_ * _1 1 C t | 

and shows that, for a regular GD-design, | C t | > 0 if b> t + v, | C t | = 0 
if b < t + v, while, if b = t+v, then 

r-»-« (r-A 1 ) v_t ~ m (r 2 -vA 2 ) m " 4_1 1 C t | 

is a perfect square. When A, = A 2 , then this case reduces to that of a 
b.i.b.d., considered in Chapter II (page 20). To obtain the same formulae 
as in that case, note that rk-v\= r- A, and that the matrix C t defined 
in the present chapter reduces to k(r-X) times that in Chapter II. 

Hence the two determinants |C e | are in the ratio of k\r-Xf . 

We have seen in the last chapter (see also [54]) that if the incidence 
matrix of a symmetric p.b.i.b.d. is not singular, and if the greatest 
common divisor of certain determinants there defined is unity, then the 

* Remember that 


AA! A 0 
K It 


11 


+ d\2 + 0, 


13 


a,,a 21 + a, 2 a 22 + 


a n d\z 

2 

a 11 a 13 

2 

a 12 a 13 


+ 


+ 

a 2 1 a 22 


a 21 a 23 


a 22 a 23 


a 11 a 21 + a 12 a 22 + a 13 a 23 

° 21 + < 4 + 4 


^ 0, and similarly 


for more rows and columns. 
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dual of the design is a P-b.i.b.d. with the same parameters. For a GD- 
design the determinants in question are those of the sub-matrices of 



The determinant of the first two rows is A 2 -A,, and that of the 
first and third rows is (A 2 - A,) n, + A, - r. In a symmetric design, 
r 2-r = A 1 n 1 +A 2 n 2 , so that the latter determinant reduces to (A, -A 2 )- 
(r 2 -vA 2 ). This is relatively prime to A 2 -A,, if r*-vA 2 is. The same 

result for GD-designs was, in a different form, proved by Connor in 
[37]. 


Shrikhande has shown in [115] that if a is an odd factor of m, and 
if a symmetric regular GD-design is constructed from a difference set 
such that multiples of m appear A, times as differences, and the other 
non-zero values A 2 times, then x 2 = (k 2 _ V A 2 ) y 2 + (_ l ) (a-1)/2 az 2 has a 

solution in integers which are not all zero. If c is an odd factor of 
n = v/m, and c and m are relatively prime, then a symmetric GD-design 
with k > A t can be constructed from a difference set as described 
above only if * 2 = (fc-A^y 2 + (_ l) (c 11/2 cz 2 has Such a so i ution ( C om- 
pare similar results for b.i.b.d’s in Chapter II.) 

In a semi-regular GD-design we have rk- vA 2 = (r- + 

(n-l)(A 1 -A 2 ) =0, i.e. r- At = n(A 2 -A 1 ), and hence A 2 >A!. 

We now have | A A' | = 0. If we omit from AA' the 2nth, 3nth ... 
mnth rows and columns (but not the nth), we obtain a product of rk 
and m matrices, but apart from the first of these they are now only of 
order n- 1. Proceeding in a similar manner to that applied before, we 
find that the determinant of the remaining matrix equals 


rk(r - A ,) n -\rk - vA 2 - A, + A 2 f'’ (r - A0 lm " fl(n “® 

= rA(r-A i r" 2m+, (A 2 -A 1 ) m_1 

which is not zero. Hence the rank of AA' is at least v-(m-l), and 
because it cannot exceed b, we have for semi-regular GD-designs 

b >v + r - m . 

If 

11 ej varieties of the first group of a GD association scheme appear 

* n the / th block of the design (/ = 1,..., b), then 2 e,- = nr = bk/m 

j= 1 



84 


(since v=mn). Further, because every pair of varieties trom the first 
group appears together in A, blocks, we have J = Ai^ j . 

From this, and from r-A, = n(A 2 - A,) (see above), we have in a semi¬ 
regular GD-design 2 ef = A^n- l)n+ nr = n 2 A 2 , and hence 

1 ( e , -2 eJbf = n%-bk 2 /m 2 = %(vX 2 -bk 2 /nm) 

3 =1 3 3 3 

- ^(vA 2 -r/c) =°. 

But a sum of squares cannot be zero unless all its terms vanish, 
and therefore e ; - = k/m for all j. The ratio k/m must, of course, be 
an integer. This applies not only to the first group, but to all of them, 
i.e. all blocks of a semi-regular GD-design contain the same number of 

varieties from all groups. 

(B) Triangular association scheme 

If the numbeT of varieties is of the form (2 ), then a so-called 
“triangular association scheme” can be constructed. We denote the 
varieties by (i,j = 1 , ..., n, i < j) , and write the scheme 

• #12 #13 • * ■ #1 n 

#12 • #23 • • • #2 n 


#1n a 2n a 3n ‘ • 

First associates are pairs of varieties in the same columns, and 
all other pairs are second associates. Fvery variety appears in two 
columns and therefore has n 1 =2(n-2) first associates. Hence n 2 = 

If two varieties are first associates, then they appear once 

in the same column, so that pj = n- 2. From this we obtain 

f n-2 n-3 \ / 4 2(n-4) \ 

^ = (n -3 { n f)) ; ^ = U„-4) (»j<) } 

If n = 3, then we have a b.i.b.d. 

We give an example of a p.b.i.b.d. with such an association scheme, 
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where w = 6, A, = 0, A^ 1. 


1 

10 

15 


2 

7 

15 


3 

8 

12 

1 

11 

14 


2 

8 

14 


3 

6 

15 

1 

12 

13 


2 

9 

13 


3 

9 

11 



4 

6 

14 


5 

6 

13 





4 

7 

12 


5 

7 

11 





4 

9 

10 


5 

8 

10 




This is the design mentioned in [19] and referred to on page 75. 

Conditions for the existence of such designs are given in [78] and 
[79]. (The latter deals also with other types of designs.) Again, the 
association scheme does not define the design. In [19] four designs 
are given for the same triangular scheme, namely that for v = 10 

W.S. Connor has shown that if in a scheme with two associate 
classes n, - 2(n-2) holds, and if the parameters of the second kind 
have the values given above, then the association scheme is triangular, 
provided n > 9. Shnkhande has shown the same for n 4 6 (see [38], 

11?] and [ 30 D’ 3X1(3 Hoffman has proved in [52] that it is also true for 
n = 7 but not for n = 8. He has also described, in [53], how to enumer¬ 
ate all counter-examples for n = 8. A complete list of these is con¬ 
tained in [31], 

If in a P-b.i.b.d. with a triangular association scheme the equality 
2{rk-v\) =n(r-\) holds, which is the same as 2 rk-nr- 
n K n- 2) A, = 0, then 2k is divisible by n, and every block contains 
precisely 2k/n varieties from each of the n rows of the association 
scheme. This is proved in [89] as follows: 

Let e tj be the number of varieties in block , from row i of the 
scheme. Then 


t e i • 
j=i 13 


r(n-l), I (M = A, 
\ 2 ] 1 


71 - f 
2 


= ("-1) [r + A 1 (n_2)]. 


We have e- = 


^ e n/b = r(n-l)/b = 2k/n, 


b 

= (n- 1) [r + \(n-2) ] - Abk Z /n z = 


so that 


and 
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„ tzl [nr + \,n(n-2)] -4vM/« 2 = + *,n(n-2)-2,k] 

n = 0 . 

Hence e# = e* = 2 k/n for all i . This proves the statement. 

For the special case b = (” 2 *)' = r = n-2, A, = 1, A 2 = 2 

this was proved, by a different method, in [llO]. In that case 2 k/n = 2. 
A computation similar to that for the GD-design shows that for a 

triangular design 

' \AA'\ = rJk(r- 2 )(A 1 +A 2 )" n(n ' 3> [r + (n-4)A 1 -(n-3)A 2 ] n V 


(C) Latin-square type designs 

Consider an orthogonal array [n 2 ,c,n, 2], equivalent to c -2 
orthogonal Latin squares of side n when c > 2. Label the columns 
from 1 to n 2 and let these labels be the varieties of a design. Two 
varieties are first associates if their columns contain, in any row, the 
same element; otherwise they are second associates. 


Example [25, 4, 5, 2 ] 
123456789 10 

0000011111 
0123401234 
0 123412340 
0 12342340 1 


11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 
222223333344444 
012340123401234 

2340 1340 1240 123 

4012312340340 12 


Second associates of 1 are 7,8, 13,15,17,19, 24, 25, and all other 
varieties are its first associates. 

An association scheme thus defined is said to be of Latin-square 
type” and is denoted by L c . L 2 , in particular, can be described equiva¬ 
lently as an arrangement of n 2 varieties in a square of side n, where 
first associates are any two varieties either in the same row, or in the 
same column. This scheme has already been mentioned in Chapter IV 
(page 61) as a modification of the rectangular association scheme. 

To determine n u take any column and its topmost element. There 
will be n—1 other columns with the same top entry. This is also true 
of the 2nd, 3rd,...,cth entry in the column. Since A,= 1, there are n, = 
(n-l)c columns that correspond to first associates. Hence n 2 = 

(n 2 - 1) - (n - 1 )c = (n - l)(n + 1 - c). 

To determine the parameters of the second kind, we first find 

Two columns of the array correspond to first associates if in any of 
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the c rows they contain the same element. Without loss of generality 
we may assume that we choose two varieties for which this is so in 


the first row. 

We must now find all those other columns which correspond to 
first associates of both varieties in the pair chosen. All those further 
columns, n-2 in number, which also have the same element in the 
first row, are first associates of both. 


Take any of the ( 2 ) Psirs of rows from the array, ignoring the 
first row. Assume that the columns of the original pair read, in such a 

pair of rows, ^ j ^ j . A column with (^j in the same two rows 
will correspond to a first associate of the variety with by virtue 
of a u and to a first associate of that of Q'j by virtue of b 2 . Hence 
this column must be counted when enumerating , as must also that 


with 


a 


\ / 

Accounting for all ( c ~ 1 ) pairs of rows apart from the first, we 
thus have 2 ^ c 2 J more columns, so that altogether p' u = 2 ( c *) 
n-2 = n+c(c- 3). From this value we obtain 


+ 



/ n + c(c-3) 
y(c-l)(rc-c + 1) 


(c-l)(n-c+1) 
(n-c)(n-c +1) 


( p2 ) . ( C< - C ~V c(n-c) \ 

,k \o(n-c) (n-cf + (c-2)/ 

It should be noted that reversing the roles of first and second asso¬ 
ciates changes these matrices from those for c into those for n-c + 1. 

We know that if n is the power of a prime, then c can be as large 
as n+ 1. Then n 2 = 0, i.e. all varieties are first associates. 

For c = 1, the parameters of the second kind are 




f 0 "-M 

\n -1 n(n-2)J 


This is a GD-design, and so is, by the foregoing remark, that of 
c ~ n. 
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For c = 2, we have 


/n -2 "-1 \ 

^ ~ \n-l ( n-l)(n-2)J 



/ 2 2(n-2)\ 

\2(n-2) (n-2) z ) 


From this we can obtain the parameters for c n 1. 

Raghavarao has proved, in [89], that if rk-vky = n(r- ,), then 
every block contains k/n varieties from each of the n rows or columns 
of the association scheme. The proof is analogous to that of the simi¬ 
lar theorem for triangular association schemes. 

For a Latin square design, similar calculations to those which we 

have carried out for GD-designs lead to 


\AA'\ = rk[r-(c-n)( A 1 -A 2 )-A 2 ] [r - c(A 1 -A 2 ) - Ag] 

Shrikhande has shown in [ 118] that when n^, 2, but not 4, and 
n, = 2n-2, pj, = n- 2, p 2 =2, then the association scheme must be 
L 2 . If n = 4, then a corresponding L 2 exists, but so does another 
association scheme with those parameters. D.M. Mesner has given 
corresponding results for L c , c > 2 (see [118] and, for large n, a 
thesis mentioned in [38]). 

If we construct a design by putting into the same block those 
varieties which have the same element in the first row (there will be 
n such blocks), then thoje which have the same element in the second 
row, etc., we obtain a p.b.i.b.d. with v=n 2 , b = cn, r=c, k. = n, A, = 1, 
A 2 =0. For c= n+ 1, we have a b.i.b.d., equivalent to a Euclidean 
plane. For c = n- 1, we have, after reversing the roles of first and 
second associates, the parameters of an L 2 scheme, because n-(n-l)+l 
= 2, and provided that n is not 4, we must, by Shrikhande’s theorem, 
be able to write the varieties into a square scheme in such a way that 
any two in the same row or in the same column are first associates. ' 

If we then construct 2 n new blocks additional to the original 
scheme by writing into the same block those varieties which are in the 
same row, or in the same column, of the square scheme, and then in¬ 
corporate these blocks into the orthogonal array [n 2 , n- 1, n,2] we ob¬ 
tain an [ n 2 , n+ 1, n, 2]. We have thus proved, following Shrikhande (see 
[l20]), that the existence of a set of n- 3 orthogonal Latin squares 
of side n > 4 implies the existence of a complete set of n-1 orthog¬ 
onal Latin squares of side n, and thus that of a PG(2,n). 
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Example 

In the [25, 4, 5, 2 ] given above, we reverse the roles of first and 
second associates. The association scheme can then be written 


1 

25 

19 

13 

7 

8 

2 

21 

20 

14 

15 

9 

3 

22 

16 

17 

11 

10 

4 

23 

24 

18 

12 

6 

5 


(Note, incidentally that this is a magic square: all rows and 
columns add up to 65.) 

We can then add the following two rows to the array: 


0 


1 2 3 4 4 0 1 2 3 3 4 0 1 2 2 3 4 0 1 1 2 3 4 0 
1 2 3 4 3 4 0 1 2 1 2 3 4 0 4 0 1 2 3 2 3 4 0 1 


side^ 6 6 roWS a ' e 6,UiValent t0 four orthogonal Latin squares of 

The following generalisation of this theorem was later proved 
by Brack in [25], If n > *72 + *’ + *V3*/2, x=n- ( -2, then a set 
of t parrwrse orthogonal Latin squares of side „ can always be 

omp eted. Also, rf n > (n-t-Zf, end if the set can be completed 
then the completion is unique. ’ 

, In [69] C.R, Nair has given necessary and sufficient conditions 

Pr ° Vided the 1 ~ * 


(1) 

( 2 ) 

( 3 ) 


Derive designs from (i) a tetrahedron, (ii) an octahedron 

^It/lines'throuS it” 110 " 8 ' ^ 3 ^ 


4(1>1 ' Z ), B (D.l,y), C(l,l,y), 

(nee Exercise 1 of Chapter II). 
Construct p.b.i.b.d’s by E. Seiden’ 


D (W), E(l,z,z) 
s methods (page 76). 
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(4) From a EG(2,3), construct a b.i.b.d. with v = 9, b= 12, fc = 3, 
r = 4, and invert it. 

(5) Determine the parameters of the following group-divisible 
design: 
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SOLUTIONS TO EXERCISES 

Chapter u 

(1) The 15 b-points are 

d(100) a(001) e(101) /(10y) *(10z) 

/c(lyl) fc(110) n(lzO) /(lyO) b(010) 

c(01z) * (111) /(lyy) o(lzl) m(lyz) 

and the 6 v-lines 

(1) x y = 0 (2) x 2 = 0 (3) x 0 + zxy =0 

(4) x 0 + x 2 = 0 (5) x 0 + yxy + zx 2 =0 (6) x 0 + x 1 + yx 2 = 0 

This gives the following b.i.b.d. : 

12 13 14 15 16 

34 26 25 23 24 

56 45 36 46 35 

(2) 3567 4671 5712 6123 7234 1345 2456 

(3) 11111111 

1 1 1-1 1 -1 -1 -1 
1-1 1 1-1 1 - 1-1 
1 - 1-1 1 1-1 1-1 
1 -1 -1 -1 1 1-1 1 
1 1 -1 -1 -1 1 1-1 

1-1 1 -1 -1 -1 1 1 

7 1 1-1 1 -1 -1 -1 1 

(4) This is a b.i.b.d. with v= 7, k = 3, A = 1, k- A= 2. The Dio- 
phantine equation x 2 = 2y 2 —z 2 has the solution x = y = z = 1. 

Chapter III 

(1) a b Q 0 12 0 12 

c P Z leads to 1 2 0 2 0 1 
R Y A 201120 
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SOLUTIONS TO EXERCISES 


Chapter III (continued) 
( 2 ) 


uX + v 



V = 


0 


X = x 
1 x 


x + l 



V = 

u =0 
1 
X 

x+l 


X = x+ 1 
0 1 x 


x+l 


0 1 

x + l x 
1 0 

X x+l 


X 

0 

X + 1 

1 


x+l 

1 

X 

0 


or, in other symbols 


0123 0123 0123 
1032 2301 3201 
2301 3210 1032 

3210 1032 2310 


(3) 


01234 

01234 

01234 

01234 

12340 

23401 

34012 

40123 

23401 

40123 

12340 

34012 

34012 

12340 

40123 

23401 

40123 

34012 

23401 

12340 


(4) 


(i) 


012 

34 

340 

12 

123 

40 

234 

01 

401 

23 


(ii) Clearly impossible, because 
the fifth symbol is missing 
in every row. 



SOLUTIONS TO EXERCISES 
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Chapter IV 
(1) v = 12 


• 

12 

13 

14 


>1 

B 

C 

21 

• 

23 

24 

D 


E 

F 

31 

32 

• 

34 

G 



1 

41 

42 

43 

• 

J 

A. 

L 



The blocks are the rows and the columns of this scheme. 


Association scheme 



First 

Associates: 

Second 

Third 

Fourth 

4 

BCHK 

IL 

EFGJ 

D 

B 

ACEL 

FG 

HIDJ 

G 

C 

ABF I 

EH 

DGKL 

J 

D 

EFGJ 

IL 

BCHK 

A 

E 

DFBL 

CJ 

AKGI 

H 

F 

DEC I 

BG 

AHJL 

K 

0 

HIDJ 

FK 

ABEL 

B 

H 

GIAK 

CJ 

DFBL 

E 

I 

GHCF 

AD 

CFJK 

L 

J 

KLDG 

EH 

ABDG 

C 

K 

J L AH 

BG 

CIDE 

F 

L 

JKBE 

AD 

CFJK 

I 



/ 2020 \ 

(P**> - f 0001 ] 

l 2020 j 
\ 0100 / 
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SOLUTIONS TO EXERCISES 


Chapter IV ( 1 ) ( continued) 
For instance: 4 

pfs 


4 A ... BCHK 
D...BCHK 

2 A...EGFJ 
I ... CK FJ 


( 2 ) 


< 

II 

cr 

II 

a\ 

r - k = 

: 3 , Aj - 2 , A 2 

= 1, n, = n 

(pi*) 

II 

O 

J) • (pf*) 

- (1 

( 3 ) B Q is the identity matrix 



B, 


B* 

% 

01011000 


00100101 

00000010 

10100100 


00011010 

00000001 

01010010 


10000101 

00001000 

10100001 


01001010 

00000100 

10000101 


01010010 

00100000 

01001010 


10100001 

00010000 

00100101 


01011000 

10000000 

00011010 


10100100 

01000000 



b,b 3 

b 2 b 3 

02022030 


00100101 

01011000 

20200203 


00011010 

10100100 

02023020 


10000101 

01010010 

20200302 


01001010 

10100001 

20300202 


01010010 

10000101 

02032020 


10100001 

01001010 

30200202 


01011000 

00100101 

03022020 


10100100 

00011010 


32122101 

23211210 

12320121 

2 \ i B i = AA' = 21231012 
1 * 21013212 

12102321 

01211232 

10122123 


SOLUTIONS TO EXERCISES 
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Chapter IV (continued) 

( 4 ) (i) abc ade bdf cef 


Association scheme 



First Second 


First 

Second 

a 

bcde f 

d 

abef 

c 

b 

acdf e 

e 

acdf 

b 

c 

abef d 

i 

bcde 

a 


•* 

co 

II 

II 

VO 

II 

> 

r = 2, 

n\ = 4 , n 2 

= i, 


A, = 1, A 2 = 0, (p.y 

= (2 1 
\ 1 0 

/*N 

S" 

II 

O 

O O 


(ii) adgj bdik cfgk 

aehk cdhl bfhj 
afil begl ceij 

v = 12 , b = 9 , k = 4 , r = 3 

n \ = 9 ' n 2 = 2 ' = 1, A.2 = 0 

Association scheme 


First Second First Second 


?| 

be g'l 

hi 


* defghijkl ac h 

" abcdefjkl gi 

C J 

ab i . 

gh 


ef i) 

kl 

e 

' abcghijkl df k 

> abedefghi jl 

/. 

de l j 

jk 


(5) 


c p/ k ) 

- & 

2 

0 

) 

,)■ 

- ( 9 
l 0 

a 

678 

d 

69 E 

£ 

89 X 

b 

79 X 

e 

6 XE 

h 

689 

c 

8 XE 

f 

79 E 

i 

78 E 



i 

67 X 
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SOLUTIONS TO EXERCISES 


Chapter /V( 5 ) ( continued ) 

Inverted: adehj bdfgh abfij bcegj acghi cdefi 
v = 10, b = 6 , r = 3 , k = 5, n y = 3 , n 2 = 6 A, = 2, A 2 = 1 


Association scheme 



First 


First 

a 

hij 

f 

bdi 

b 

fgj 

g 

bch 

c 

egi 

h 

adg 

d 

efh 

i 

acf 

e 

cdj 

i 

abe 



(pp = 

(2 

2 

4 , 

) (p J* } = (23) 

(6) 

First: 14 24 

34 


This is a (6 3 , 9 2 ) 


15 25 

35 



16 26 

36 




Association scheme (1st associates): 




4 5 

6 

5? 12 3 



6j 

: 123 

147 456 

258 789 369 

«, = 4 , 

n 2 = 4 , 

II 

ro 

II 

O 

II 

r*i 

(2 2)- 

(p * )= (2 ?)■ 


Second: 


2 

3 

4 
1 
4 
3 


7 12 

8 10 

6 11 
8 11 
7 9 

5 12 


5 

6 
8 
6 
5 
7 


10 

12 

9 

9 

11 

10 


This is a 

( 12 ,). 





SOLUTIONS TO EXERCISES 
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Chapter IV ( 6 ) ( continued ) 


( 1 ) 



Association scheme 



First 

Second 

Third 

1 

6 7 8 10 11 12 

2 3 4 

5 9 

2 

5 7 8 9 11 12 

13 4 

6 10 

3 

5 6 8 9 10 12 

1 2 4 

7 11 

4 

5 6 7 9 10 11 

12 3 

8 12 

5 

2 3 4 10 11 12 

6 7 8 

1 9 

6 

1 3 4 9 11 12 

5 7 8 

2 10 

7 

1 2 4 9 10 12 

5 6 8 

3 11 

8 

1 2 3 9 10 11 

5 6 7 

4 12 

9 

2 3 4 6 7 8 

10 11 12 

1 5 

10 

1 3 4 5 7 8 

9 11 12 

2 6 

11 

1 2 4 5 6 8 

9 10 12 

3 7 

12 

1 2 3 5 6 7 

9 10 11 

4 8 

?r V 

(i) 

123 124 134 234 

v = b = 4, k = r = 

(A b.i.b.d.) 

3, A = 2. 


(ii) 

134 135 146 156 

235 234 256 

246 


v - 6, b - 8, k - 3, r = 4, A 1 = 2, A 2 = 0, n 1 = 4, n 

(p * ) = (l o) ’ (p ^ } = (o o) ■ 


= 1 


(2) Omit (0,0). The remaining points are 

“(0,1) 6(0,2) c(l,0) d(l,l) e(l,2) /(2,0) *(2,1) 6(2,2) 
The remaining lines (blocKs) are 

= l...c<ie x 2 = l...adg x t + x 2 = 1 ... ach x^+2x 2 = 

1 ••• beg X] = 2 ... fgh x-£ — 2 ... beh x^ + x 2 = 2 ... bdf 
x, + 2x 2 = 2 ... ae/ 

Association scheme: 

The following pairs are second associates: 
a ~b, c-f, d-h, e-g. 

The points A, 3 ,C,D, and E form a conic, and their nucleus is 
F(0 ,l f l) (cf. page 4). Let the remaining 15 points be denoted 


(3) 
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SOLUTIONS TO EXERCISES 


Chapter V (3) ( continued ) 

by l.c. letters as in Exercise II. 1. Then we have the following 
15 lines, each of which contains two of A, B, ..., F and three other 


points: 

*o + y* f + *2 = 0 ABeln *o = 0 BFabc 

x 0 + x y =0 ACahi x 0 + zx 2 = 0 CDbfl 

x^ + yx 2 =0 ADcdk x 0 + zx, + x 2 = 0 CEcej 

*o + y * 2 = 0 AEbgm Xq + yx, + yx 2 = 0 CFgkn 

Xq + zx, + zx 2 = 0 AFfjo X 0 + yx 1 =0 DEano 

x, + zx 2 =0 BCdmo x 0 + x, + x 2 = 0 DFehm 

x Q + zx, + yx 2 = 0 BDgij x, + x 2 =0 EFd*7 

x 0 + x, + zx 2 = 0 BEfhk 

This is the first design. 

Second design: 


6 points with 5 points on each line: 


II 

o 

adefg 

O 

II 

CM 

X 

bdhjn 

x 0 + x 2 = 0 

beiko 

x 0 + zx, =0 

ajklm 

x 0 + x, + yx 2 = 0 

cghlo 

Xq + yx, + zx 2 = 0 

cfimn 


(4) The b.i.b.d. is 

123 145 167 189 

468 269 258 247 
579 378 349 356 

Inverting it, we obtain 


1 7 8 9 3 6 8 11 1 4 6 10 1 5 11 12 
25 8 10 345 9 247 11 26 9 12 
3 7 10 12 



This is a group-divisible p.b.i.b.d., with the groups as follows 
(after exchanging first and second associates): 

1 2 3 4 8 12 5 6 7 9 10 11 . 


SOLUTIONS TO EXERCISES 


99 


Chapter V (continued) 

(5) m = 3, n = 2, v = 6, b = 3, r = 2, k = 4 

ra, = 1, m 2 = 4, = 2, A. 2 = 




100 


bibliography 


[ 1 ] 

[ 2 ] 

[3] 

[43 

[5] 

[ 6 ] 
[7] 
[S] 

[9] 

[ 10 ] 


Berman, G., A three parameter family of partially balanced incom¬ 
plete block designs with two associate classes. Proc. Amer. 
Math. Soc., 6 (1955), 490-493 

Bhattacharya, K.N., A new solution in symmetrical balanced in¬ 
complete block designs (v=b= 31, r=k=10, A = 3). Sankhya, 

7 (1946), 423-424 

Bose, R.C., On the application of the properties of Galois fields 
to the construction of hyper-Graeco-Latin squares. Sankhya, 

3 (1938), 323-338 

_ _ On the construction of balanced incomplete block 

designs. Ann. Eugenics, 9 (1939), 353-399 

_ } On some new series of balanced incomplete block 

designs. Bull. Calcutta Math. Soc., 34 (1942), 17—31 

_, A note on Fisher’s inequality for balanced incom¬ 
plete block designs. Ann. Math. Statist., 20 (1949), 619—620 

_, A note on orthogonal arrays (abstract). Ann. Math. 

Statist., 21 (1950), 304-305 

_, Partially balanced incomplete block designs with 

two associate classes involving only two replications. 
Calcutta Statist. Ass. Bull., 3 (1951), 120—125 

-, On the application of finite projective geometry for 

deriving a certain series of balanced Kirkman arrangements. 
Calcutta Math. Soc. Golden Jubilee Commemoration volume 
(1958-1959), 1959, 341-354 

-, Strongly regular graphs, partial geometries and 

partially balanced designs. Pacific J. Math., 13 (1963), 
389-419 


[11] Bose, R.C. and Bush, K.A., Orthogonal arrays of strength 2 and 

3. Ann. Math. Statist., 23 (1952), 508-524 

[12] Bose, R.C. and Clatworthy, W.H., Some classes of partially 

balanced incomplete block designs. Ann. Math. Statist., 

26 (1955), 212-232 



101 


[13] Bose, R.C., Clatworthy, W.H. and Shrikhande, S.S., North 

Carolina Agric. Exp. Station Techn. Bull. no. 107 (1954) 

[14] Bose, R.C. and Connor W.S., Combinatorial properties of group 

divisible incomplete block designs. Ann. Math. Statist., 

23 (1952), 367-383 

[15] Bose, R.C. and Mesner, D.M., On linear association algebras 

corresponding to association schemes of partially balanced 
designs. Ann. Math. Statist., 30 (1959), 21-38 

[16] Bose, R.C. and Nair, K.R., Partially balanced incomplete block 

designs. Sankhya, 4 (1939), 337-372 

[17] - --- On complete sets of Latin squares. 

Sankhya, 5 (1941), 361-382 

[18] - ,-- Resolvable incomplete block designs 

with two replications. Sankhya, A24 (1962), 9-24 

[19] Bose, R.C. and Shimamoto, T., Classification and analysis of 

partially balanced incomplete block designs with two asso¬ 
ciate classes. J. Amer. Statist. Ass., 47 (1952), 151-184 

[20] Bose, R.C. and Shrikhande, S.S., On the falsity of Euler’s con¬ 

jecture about the non-existence of two orthogonal Latin 
squares of order 4*+2. Proc. Nat. Acad. Sci. E.S.A., 45 
(1959), 734-737 

[21] - ,---, On the construction of sets of 

mutually orthogonal Latin squares and the falsity of a con¬ 
jecture of Euler. Trans. Amer. Math. Soc., 95 (I960), 
191-209 

[22] Bose, R.C., Shrikhande, S.S. and Bhattacharya,K.N., On the con¬ 

struction of group divisible incomplete block designs. Ann. 
Math. Statist., 24 (1953), 167-195 

[23] Bose, R.C., Shrikhande, S.S. and Parker, E.T., Further results 

on the construction of mutually orthogonal Latin squares 
and the falsity of Euler’s conjecture. Can. J. Math., 12 
(1960), 189-203 

[24] Brownlee, K.A. and Loraine, P.K., The relationship between 

finite groups and completely orthogonal squares, cubes 
and hypercubes. Biometrika, 35 (1948), 277-282 

[25] Brack, R.H., Finite nets II: Uniqueness and embedding. Pacific 

J. Math., 13 (1963), 421-457 

[26] Brack, F.H. and Ryser, H.J., The non-existence of certain finite 



102 


projective planes. Can. J. Math., 1 (1949), 88—93 

[27] Bush, K.A., Orthogonal arrays of index unity. Ann. Math. Statist., 

23 (1952), 426-434 

[28] _, A generalisation of a theorem due to MacNeish. Ann. 

Math. Statist., 23 (1952), 293-295 

[29] Carmichael, R.D., Introduction to the Theory of Groups of Finite 

Order. Ginn & Co., Boston, 1937 

[30] Chang, L.C., The uniqueness and non-uniqueness of the triangular 

association scheme. Science Record, 3 (1959), 604—613 

[31] -, Association schemes of partially balanced designs 

with parameters v=28, n, = 12, n 2 = 15 and P 1 2 1 =4. Science 
Record, 4 (1960), 12—18 

[32] Chowla, S., Erdos, P. and Straus, E.G., On the maximal number 

of pairwise orthogonal Latin squares of a given order. Can. 

J. Math., 2 (1950), 92-99 

[33] Chowla, S. and Ryser, H.J., Combinatorial problems. Can. J. 

Math., 2 (1950), 92-99 

[34] Clatworthy, W.H., A geometrical configuration which is a par¬ 

tially balanced incomplete block design. Proc. Amer. Math. 
Soc., 5 (1954), 47-55 

[35] -, Contributions on partially balanced block de¬ 

signs with two associate classes. Nat. B. Standards Appl. 
Math. Series no. 47 (1956) 

[36] Connor, W.S., On the structure of balanced incomplete block 

designs. Ann. Math. Statist., 23 (1952), 57—71 

[37] -, Some relations among the blocks of symmetrical 

group divisible designs. Ann. Math. Statist., 23 (1952), 
602-609 

[38] -, The uniqueness of the triangular association 

scheme. Ann. Math. Statist., 29 (1958), 262—266 

[39] Erdos, P. and Kaplansky, I., The asymptotic number of Latin 

rectangles. Amer. J. Math., 68 (1946), 230-236 

[40] Euler, L., Recherches sur une nouvelle espece de carres 

magiques. Verh. Genoot. Wetensch. Vlissingen, 9 (1782), 
85-239 

[41] Evans, T., Embedding incomplete Latin squares. Amer. Math. 

Monthly, 67 ( 1960), 958-961 


103 


/\ 

[42] Fisher, A., A system of confounding for factors with more than 

two alternatives, giving completely orthogonal cubes and 
higher powers. Ann. Eugenics, 12 (1945), 283—290 

[43] Fisher, R.A., An examination of the different possible solutions 

of a problem in incomplete blocks. Ann. Eugenics, 10 (1940), 
52-75 

[44] -, The Design of Experiments. Oliver & Boyd, 

Edinburgh, 1942 (3rd edn) 

[45] Fisher, F.A. and Yates, F., The 6x6 Latin squares. Proc. 

Camb. Phil. Soc., 30 (1934), 492-507 

[46] _,_, Statistical Tables for Biological, 

Agricultural and Medical research. Oliver & Boyd, 

Edinburgh, 1953 (4th edn) 

[47] Hall, Mr, Jr., An existence theorem for Latin squares. Bull. Amer. 

Math. Soc., 51 (1945), 387-388 

[48] -, A survey of combinatorial analysis. Surveys in 

Applied Mathematics, IV: Some aspects of analysis and 
probability. Wiley, 1958 

[49] Hanani, Haim, The existence and construction of balanced in¬ 

complete block designs. Ann. Math. Statist., 32 (1961), 
361-386 

[50] Harshbarger, B., Preliminary report on the rectangular lattices. 

Biometrics Bull., 2 (1946), 115—119 

[51] -, Triple rectangular lattices, Biometrics, 5 

(1949), 1-13 

[52] Hoffman, A.J., On the uniqueness of the triangular association 

scheme. Ann. Math. Statist ., 31 (1960), 492-497 

[53] -, On the exceptional case in a characterization of 

the arcs of a complete graph. IBM Res. and Dev., 4 (1960), 
487-496 

[54] -, On the duals of symmetric partially balanced 

designs. Ann. Math. Statist., 34 (1963), 528-531 

[55] Husain, Q.M., Symmetrical balanced incomplete block designs 

with A =2, k = 7 or 9. Bull. Calcutta Math. Soc., 37 
(1945), 115-123 

^ the totality of the solutions for the symmetrical 

incomplete block design. Sankhya, 7 (1945/6), 204-208 



104 


[57] Husain, Q.M., Impossibility of the symmetrical incomplete block 

design with A =2, k = 7. Sankhya, 7 (1946), 317—322 

[58] -, Structure of some incomplete block designs. 

Sankhya, 8 (1948), 381—383 

[59] -, A note on symmetrical balanced incomplete block 

designs with k = 9, A= 2. Bull. I.S.I., 38 IV (1961), 11-16 

[60] Kishen,-K~, On the construction of Latin and hyper-Graeco-Latin 

cubes and hypercubes. J. Indian Soc. Agric. Stat., 2 (1950), 
20-48 

[61] Levi, F.W., Finite Geometrical Systems. Univ. of Calcutta, 1942 

[62] MacNeish, H.F., Euler squares. Ann. Math. (ser. 2), 23 (1922), 

221-227 

[63] Majumdar, K.N., On some theorems in combinatorics relating to 

incomplete block designs. Ann. Math. Statist., 24 (1953), 
377-389 

[64] Mann, H.B., The construction of orthogonal Latin squares. Ann. 

Math. Statist., 13 (1942), 418-423 

[65] -, On the construction of sets of orthogonal Latin 

squares. Ann. Math. Statist., 14 (1943), 401—414 

[66] -, On orthogonal Latin squares. Bull. Amer. Math. Soc., 

50 (1944), 249-257 

[67] Mendelsohn, N.S., Dulmadge, A.L., Johnson D.M. and Parker, E.T., 

Construction of m.o.l. squares. Abstract 567-4, Notices 
Amer. Math. Soc., 7 (I960), 208 

[68] Moore, E.H., Concerning triple systems. Math. Ann., 43 (1893), 

271-285 

[69] Nair, C.R., On the methods of block section and block inter¬ 

section applied to certain p.b.i.b. designs. Calcutta 
Statist. Ass. Bull., 11 (1962), 49-54 

[70] Nair, K.R., Partially balanced incomplete block designs involv¬ 

ing only two replications. Calcutta Statist. Ass. Bull., 

3 (1950), 83-86 

[71] -, Some 3-replicate partially balanced incomplete block 

designs. Calcutta Statist.Ass. Bull., 4(1951), 39—42 

[72] -- Rectangular lattices and partially balanced incom¬ 

plete block designs. Biometrics, 7 (1951), 145-154 

[73] Nair, K.R. and Rao, C.R., Incomplete block designs for experiments 



105 


[74] 

[75] 

[76] 

[77] 

[78] 

[79] 

[80] 
[81] 
[82] 

[83] 


involving several groups of varieties. Science and Culture, 

7 (1942), 625 

Nandi, H.K., On the relation between certain types of tactical 
configurations. Bull. Calcutta Math. Soc., 37 (1945), 92-94 

---, Enumeration of non-isomorphic solutions of balanced 

incomplete block designs. Sankhya, 7 (1945/6), 305—312 

-, A further note on non-isomorphic solutions of incom¬ 
plete block designs. Sankhya, 7 (1945/6), 313—316 

Norton, H.W., The 7x7 squares. Ann. Eugen., 9 (1939), 269—307 

Cgawa, J., A necessary condition for existence of regular and 
symmetrical experimental designs of triangular type with 
partially balanced incomplete block designs. Ann. Math. 
Statist., 30 (1959), 1063-1071 

-, On a unified method of deriving necessary conditions 

for existence of symmetrical partially balanced block de¬ 
signs of certain types. I.S.I. Bull., 38 IV (1961), 43-57 

Ostrowski, R.T. and van Duren K.D., On a theorem of Mann on 
Latin squares. Math. Comp., 15 (1961), 293-295 

Parker, E.T., Construction of some sets of mutually orthogonal 
Latin squares. Proc. Amer. Math. Soc., 10 (1959), 946-949 

, A computer search for Latin squares orthogonal to 
Latin squares of order 10. Abstract 564-71, Notices Amer. 
Math. Soc., 6 (1959), 798 

- , Orthogonal Latin squares. Proc. Nat. Acad Sci 

E.S.A., 45 (1959), 859-862 


[84] 

[85] 

[ 86 ] 

[87] 

[ 88 ] 


, Computer study of orthogonal Latin squares of 
order 10. Computers and Automation, 11 (1962), 33—35 

-—, Non-extendibility conditions on mutually orthogonal 
Latin squares. Proc. Amer. Math. Soc., 13 (1962), 219-221 
Plackett, R.L. and Burman, J.P., The design of optimum multi¬ 
factorial experiments. Biometrika, 33 (1943-6), 305-325 

Raghavarao, D., A note on the construction of GD designs from 

hyper-Graeco-Latin cubes of the first order. Bull. Calcutta 
Statist.Ass., 9 (1959), 67-70 

- , A generalization of group divisible designs 

Ann. Math. Statist., 31 (1960), 756-765 


, On the block structure of certain p.b.i.b. designs 


[89] 



106 


with two associate classes having triangular and L 2 asso¬ 
ciate schemes. Ann. Math. Statist., 31 (1960), 787—791 

[90] Rao, C.R., Hypercubes of strength “d” leading to confounded 

designs in factorial experiments. Bull. Calcutta Math. Soc., 
38 (1946), 67-78 

[91] _ , A general class of quasi-factorial and related designs. 

Sankhya, 17 (1956-7), 165-174 

[92] _, Combinatorial arrangements analogous to orthogonal 

arrays. Sankhya, 23A (1961), 283—286 

[93] Rao, P.V., On the construction of some partially balanced in¬ 

complete block designs with more than three associate 
classes. Bull. Calcutta Statist. Ass. 9 (1959), 87—92 

[94] _ # The dual of a balanced incomplete block design. Ann. 

Math. Statist., 31 (I960), 779—785 

[95] Ray-Chaudhuri, D.K., Application of the geometry of quadrics 

for constructing p.b.i.b. designs. Ann. Math. Statist., 33 
(1962), 1175-1186 

[96] Reiss, M., Uber eine Steinersche kombinatorische Aufgabe, 

welche im 45. Band dieses Journals, Seite 181, gestellt 
worden ist. Crelle’s J. reine und angew. Math., 56 (1859), 
326-344 

[97] Rojas, B. and White, R.F., The modified Latin square. J. Royal 

Stat. Soc. (B), 19 (1957), 305-317 

[98] Rouse-Ball, W.W., Mathematical recreations and essays (revised 

by H.S.M. Coxeter). MacMillan, London, 1942 

[99] Roy, J. and Laha, R.G., Two associate partially balanced de¬ 

signs involving three replications. SankhyS, 17 (1956-7), 
175-184 

[100] Roy, P.M., Hierarchical group divisible incomplete block designs 

with m associate classes. Science and Culture, 19(1953-4), 
210-211 

[101] _, On the method of inversion in the construction of 

partially balanced incomplete block designs from the cor¬ 
responding balanced incomplete block designs. Sankhya, 

14 (1954), 39-52 

[102] -, On the properties and construction of HGD designs 

with m associate classes. Calcutta Statist. Ass. Bull., 11 
(1962), 10-38 



107 


[103] Ryser, H.J., A combinatorial theorem with application to Latin 

rectangles. Proc. Amer. Math. Soc., 2 (1951), 550-552 

[104] Sade, A., An omission in Norton’s 7 x 7 squares. Ann. Math. 

Statist., 22 (1951), 306-307 

[105] Schiitzenberger, M.P., A non-existence theorem for an infinite 

family of symmetric block designs. Ann. Eugen., 14 (1949), 
286-287 

[106] Seiden, E., A remark on the geometrical method of construction 

of an orthogonal array. Ann. Math. Statist., 25 (1954), 
177-178 


[107] 

[108] 


-, On the maximum number of constraints of an 

orthogonal array. Ann. Math. Statist., 26 (1955), 132-135 

-, Further remark on the maximum number of constraints 

of an orthogonal array. Ann. Math. Statist., 26 (1955), 
759-763 


[109] 


-, On a geometrical method of construction of partially 

balanced designs with two associate classes. Aiin. Math. 
Statist., 32 (1961), 1177-1180 


L110J --, On necessary conditions for the existence of some 

symmetrical and unsymmetrical triangular p.b.i.b. designs 
and b.i.b. designs Ann. Math. Statist., 34 (1963), 348—351 

[111] Shah, B.V., On a generalization of the Kronecker product design. 

Ann. Math. Statist., 30 (1959), 48-54 

[112] Shnkhande, S.S., The impossibility of certain symmetrical 

balanced incomplete block designs. Ann. Math. Statist 
21 (1950), 106-111 


[113] 

[114] 


-, On the dual of some incomplete block designs 

Biometrics, 8 (1952), 66-72 

—-- , The non-existence of certain affine resolvable 

balanced incomplete block designs. Can. J. Math 5 (1953) 

A 13 /no ‘ J 


[115] 

[116] 

[117] 


-, Cyclic solutions of symmetrical GD-designs. 

Calcutta Statist. 4ss. Bull., 5 (1953), 36-39 

, Affine resolvable balanced incomplete block 
designs and non-singular GD-designs. Calcutta Statist. Ass. 
Bull., 5 (1954), 139-141 

~ 7 ' a c haracterization of the triangular asso¬ 

ciation scheme. Ann. Math. Statist., 30 (1959), 39-47 



108 


[118] Shrikhande, S.S., The uniqueness of the L 2 association scheme. 

Ann. Math. Statist., 30 (1959), 781—798 

[119] _, Relations between certain incomplete block de¬ 

signs. Essays in honor of Harold Hotelling. Stanford U.P., 1960 

[120] -, A note on mutually orthogonal Latin squares. 

Sankhyd, A, 23 (1961), 115-116 

[121] Shrikhande, S.S. and Singh, N.K. On a method of constructing 

symmetrical balanced incomplete block designs. Sankhya 
A, 24 (1962), 25-32 

[122] Smith, H.J.S., On systems of linear indeterminate equations and 

congruences. Phil. Trans. Royal Soc. London, 151(1961—2), 
293-326 

[123] Sprott, D.A., Some series of partially balanced incomplete block 

designs. Can. J. Math., 7 (1955), 369—381 

[124] _, A series of symmetrical group divisible incomplete 

block designs. Ann. Math. Statist., 30 (1959), 249—250 

[125] Stevens, W.L., The completely orthogonalized Latin square. Ann. 

Eugen., 9 (1939), 82-93 

[126] Tarry, G., Le probleme des 36 officiers. C.R. Acad. Franc, pour 

VAvancement de Science Naturel, 1 (1900), 122—123; 2 
(1901), 170-203 

[127] Thompson, W.A. Jr., A note on p.b.i.b. design matrices. Ann. 

Math. Statist., 29 (1958), 919-922 

[128] Todd, J.A., A combinatorial problem. J. Math. Phys., 12 (1932—3), 

321-333 

[129] Vajda, S., Patterns and Configurations in Finite Spaces. Griffin, 

London, 1967 

[130] Yates, F., Complex experiments. J. Royal Stat. Soc. Suppl., 2 

(1935), 181-247 

[131] -, Incomplete randomized blocks. Ann. Eugen., 7(1936), 

121-140 



109 


INDEX 


Abelian Group, 1 
affine resolvable, 7 
assembly, 40 
association matrix, 62 
— scheme, 54 
automorphism, 2 


Balanced design, 46 
b.i.b.d., 7 

block intersection, 11 
— section, 11 


hyperbola, 9 

Hyper-Graeco-Latin square, 25 
Hyper-Latin cube, 33 
hyperplane, 4 


Incidence matrix, 12 
incomplete Latin square, 24 
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